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I.  INTRODUCTION 


Since  the  prediction  of  the  pion  by  Yukawa  [Yu-35 1  in  1935  and 
its  production  in  the  University  of  California  Radiation  Laboratory  by 
Gardner  and  Lattes  in  1948,  there  has  been  extensive  theoretical  and 
experimental  work  with  the  pion  as  a  probe  of  nuclear  structure.  Such 
extensive  work  was  possible  and  is  continuing  due  to  the  construction  and 
subsequent  operation  of  three  meson  production  facilities:  the  Los 
Alamos  Meson  Physics  Facility  (LAMPF)  in  the  United  States,  the  Swiss 
Institute  for  Nuclear  Research  (SIN)  in  Switzerland,  and  the  Tri- 
University  Meson  Facility  (TRIUMF)  in  Canada.  LAMPF  is  the  highest 
intensity  meson  production  facility.  However,  all  three  facilities  were 
constructed  to  produce  high  energy  resolution  secondary  beams  of  pions 
required  in  the  study  of  closely  spaced  (<  l  MeV  separation)  excited 
levels  of  nuclei  with  sometimes  extremely  small  pion-induced  inelastic 
scattering  cross  sections. 

The  pion  is  a  pseudoscalar,  J^-CT,  isovector,  T-I,  meson  and  is 
the  lightest  known  strongly  interacting  particle  (mn  -  139.6  MeV).  As  a 
field  quantum  of  the  nuclear  force,  the  pion  is  responsible  for  the  long- 
range  part  of  the  hadronic  potential.  Many  of  the  pion's  properties  make 
it  an  extremely  useful  hadronic  probe  due  to  resulting  theoretical 
simplifications  in  the  analysis  of  pion-nucleus  scattering  data,  and  the 
ability  to  perform  experiments  with  the  pion  which  complement  those  using 
other  hadronic  and  electromagnetic  probes  or  which  cannot  be  done  with 
any  other  probe:  (1)  There  are  three  charge  states  of  the  pion,  ir+,  w°, 


and  *  ;  thus,  one  can  use  the  pion  for  elastic  scattering,  single-charge- 
exchange  (SCX)  scattering,  and  double-charge-exchange  (DCX)  scattering. 
The  latter  process  can  lead  to  nuclear  levels  with  T,  Tz  values  not 
obtainable  with  a  T-l/2  hadronic  probe  such  as  the  proton.  (2)  The  pion 
possesses  spin  zero.  As  a  boson,  it  can  be  absorbed  by  clusters  of 
nucleons  in  the  nucleus,  and  pion  absorption  experiments  may  therefore 
provide  useful  information  about  nucleon  correlations  within  the  nucleus. 
Furthermore,  the  pion's  zero  spin  results  in  a  considerably  simpler  two- 
body  scattering  amplitude  for  pion-nucleon  reactions  than  for  nucleon- 
nucleon  reactions.  Also,  along  with  the  three  charge  states,  the  pion's 
spinlessness  allows  an  easier  parameterization  of  the  scattering 
amplitude  and  fits  to  the  experimental  i»+p-*-x+p,  ir~p-*-ir-p,  and  ir~p-M»°n 
data.  Thus,  the  most  important  pion-nucleon  phase  shifts  have  been 
determined  accurately  over  a  large  incident  pion  energy  region 
(T^dab)  <  500  MeV).  (3)  The  pion's  light  mass  (M/7  of  the  mass  of  a 
nucleon)  and  its  distinguishability  from  nucleons  means  that  recoil 
effects,  projectile  structure,  and  particle  exchange  between  the 
projectile  and  target  are  not  as  important  to  consider  or  do  not  need  to 
be  considered  in  a  theoretical  treatment  of  pion-nucleus  scattering  as  is 
the  case  for  other  hadronic  probes  such  as  protons  and  alpha  particles. 
(4)  Lastly,  the  total  x+p  and  Tr~p  cross  sections  (see  Fig.  1-1)  exhibit 
many  resonances  resulting  in  strong  energy-  and  charge-dependent  cross 
sections.  The  total  cross  section  for  p  +  p  shown  in  Fig.  1-1  shows  no 


resonant  structure.  At  low  to  medium  incident  pion  energies  of 
100  <  T  <  300  MeV,  the  x+p  and  n”p  cross  sections  are  dominated  by  the 


T,  (MoV) 


Fig.  1-1:  Total  cross  sections  for  w+p  (solid  line)  and  *  p  (dot-dashed 
line)  scattering  as  a  function  of  lab  kinetic  energy  and  the  total  pp 
cross  section  for  comparison  (Ei-80], 
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4(1232)  resonance  (T  **  110-122  MeV)  which  occurs  when  the  phase  shift  of 
the  1-1 ,  J-3/2,  T-3/2  partial  wave  rises  through  90°  at  »  195  MeV. 

The  4(1232)  resonance  Is  extremely  Important  in  pion-nucleus 
scattering.  The  most  prevalent  theoretical  description  of  pion-nucleus 
inelastic  scattering  uses  the  distorted-wave  Impulse  approximation 
(DWIA) ,  treating  nuclear  transitions  in  the  Born  approximation.  Such  a 
description  involves  the  basic  pion-nucleon  interaction;  therefore,  the 
energy,  spin,  and  isospin  dependences  of  the  pion-nucleus  Interaction 
result  from  the  energy,  spin,  and  lsospin  dependences  of  the  pion-nucleon 
interaction,  which  is  dominated  by  formation  of  the  4(1232)  resonance. 
Thus,  to  understand  the  utility  of  the  pion  as  a  probe  of  nuclear 
structure,  one  needs  to  look  at  pion-nucleon  scattering  in  the  energy 
region  of  the  4(1232)  resonance. 

The  center-of-mass  differential  cross  section  for  the  elastic 
scattering  of  two  particles  with  spin  Sj  and  s2  is  [Ta-72] 

do/dftc.m/P'.X'^P.X)  -  |  <X" l FCp'^p)  |  X> |  2  ,  (1-1) 

where  p,  p'  are  the  initial  and  final  relative  momenta,  |x>.  I x'>  ate  the 

general  initial  and  final  spin  states  of  the  two  particle  system 

( |  x>  ■  JxclP).  ai*d  F(p'«-p)  is  the  scattering  amplitude  matrix  consisting 

of  the  I(2$1  +  l)(2s2  +  l)]2  individual  scattering  amplitudes  f^'^(p'+p). 

|  C>  labels  the  basis  vector  )m  m2>  or  |s,m>-  [  |ra,  ,m2><s1  s2m,m  |sm>. 

,mo  ^ 

1  *  ♦  ♦ 

For  the  case  of  pion-nucleon  scattering  (s^O,  s2-l/2)#  F(p'«*p)  is  a 
(2  «  2)  matrix  and  its  most  general  form  which  is  compatible  with 
invariance  under  rotations  and  parity  is  [Ta-72] 
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F(p'<-p)  -  a(Ep,0)I  +  ib(Ep,6)n»o  ,  (1-2) 

where  9  is  the  scattering  angle,  Ep  is  the  center-of-mass  energy,  I  is 
the  unit  (2  *  2)  matrix,  a  denotes  the  Pauli  spin  matrices,  and  n  is  the 
unit  normal  to  the  scattering  plane  (n  -  (p  x  p')/(|p  x  p'|)).  With  this 
form  for  F(p'+p),  the  unpolarized  center-of-mass  differential  cross 
section  for  elastic  pion-nucleon  scattering  is 

d°/dnc.m.<P'*P)  “  |a(Ep,0)l2  +  I b(Ep,9) | 2  .  (1-3) 


Both  orbital  and  total  angular  momentum  are  conserved  as  well  as 
the  total  isospin  in  pion-nucleon  scattering.  One  can  thus  expand 
F(p'*p)  in  partial  waves  according  to 

(21  +  l)yP*,j«2T,2j(P>pi(cos9>  *  (!-*) 


F(p'*p)  -  |Qt^ 


where  QT  and  Pjj  project  onto  states  of  definite  total  isospin  T  and 

total  angular  momentum  J,  respectively,  and  the  partial-wave  amplitudes 
0 

ogj  2J  are  expressed  in  terms  of  the  phase  shifts  according  to 


ei5tT,2J<P)sin{S£T>2J(p)l 
~~~  P 


(1-5) 


The  total  isospin  is  T-l/2  or  3/2  and  the  total  angular  momentum  is 
J-i+1/2  or  t-l/2.  The  projection  operators  are 


1/3(1  -  I*T) 


Q3/2  “  1/3(2  +  I*t)  , 


( I-6a) 


Ql/2  “ 


P4,J*4-l/2 


♦  ♦ 
4  -  a*L 

24  +  1 


♦  ♦ 

„  4  +  1  +  a*L 

P4,J-4+l/2 - 2t~n - 


(I-6b) 


♦  + 
with  I  Che  pion  Isospln  operator  and  t  che  nucleon  Pauli  Isospin 

operator.  Substituting  these  expressions  Into  Eq.  (1-4),  keeping  only 

the  4*0  (s  wave)  and  4-1  (p  wave)  terms  (since  these  are  the  most 

Important  partial-wave  amplitudes  in  the  energy  range  0  <  <  345  MeV 

[Ei-80 ] ) ,  and  comparing  to  Eq.  (1-2), 

F(p'+p)  -  (aQ  +  a^'t)!  +  i(bQ  +  bt !•  r )n-a  ,  (I-7a) 


with 

aQ  -  1/3(0?^  +  20^)  +  l/3(m{fl  +  20^!  +  2a{  3  +  4a^  3)cos0  ,  (I-7b) 

-  l/3(-atj  ^  i  +  a3,i)  +  l/3(-ai>i  +  <*3>i  “  2c,|t3  +  2aj>3)cos9  ,  (I-7c) 

bQ  -  l/3(a|>3  “  <*1^  +  2a3)3  "  2a^>1)sin8  ,  (I-7d) 

bt  =  l/3(~ai>3  +  «}fi  +  ^^3  "  a3ti)sin9  .  (I-7e) 

aQ  and  bQ  are  the  isoscalar  terms,  a^  and  bj  the  isovector  terms,  and  b0 
and  bt  the  spin-flip  terms. 

As  mentioned  previously,  for  100  <  Tn  <  300  MeV,  pion-nucleon 
scattering  is  dominated  by  the  4(1232)  resonance.  Therefore,  the  largest 
phase  shift  in  this  energy  range  is  Si  With  Eqs.  (1-3)  and  (1-7)  and 
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For  example,  keeping  only  the  parcial-wave  amplitude,  Eq. 

(1-7)  la 

F(p'«-p)  -  1/30^3(2  +  I*t)[2cos0  +  in*osin0]  .  (1-11) 

Slclllano  and  Walker  [Sl-81],  following  the  approach  of  Gupta  and  Walker 
[Gu-76]  and  using  a  single-step,  impulse-approximation  mechansim,  showed 
that  the  differential  cross  section  for  plon-nucleus  inelastic  scattering 
retains  the  same  basic  energy,  isospin,  and  spin  dependence.  For 
unnatural-parity  transitions  (Ax*(-l)*J+l ,  where  Ax  is  the  parity  change 
and  J  is  the  total  angular  momentum  transfer)  the  spin  transfer  is  AS-1, 
and  such  transitions  must  proceed  through  the  n*o  operator  of  Eq.  ( I— 11). 
As  a  result,  the  constant-momentum-transfer  (q)  excitation  functions 
should  decrease  with  increasing  incident  pion  energy  according  to  sin20. 
Such  an  energy  dependence  was  first  observed  [Co-79]  for  unnatural-parity 
transitions  to  the  4”,  19.25-MeV  state  and  2”  structure  at  *18.4  MeV  in 
12C  and  later  [Se-81]  to  the  9/2+,  9.5-MeV  state  in  i3C.  Microscopic 
DWIA  calculations,  using  harmonic  oscillator  forms  for  the  spin 
transition  densities,  reproduce  well  both  the  excitation  functions  and 
the  shapes  of  the  angular  distributions  for  the  4"  state  of  i2c  [Co-84] 
and  the  9/2+  state  of  i3C  [Se-81a,Se-82] . 

Transitions  to  states  of  stretched  configuration,  (jpjh~l)Js  with 
j  ■  ip  +  1/2,  -  ih  +  1/2,  and  Jg  -  jp  +  jh,  are  an  important  subclass 

of  unnatural-parity  transitions.  For  these  transitions  only  the 
L  ■  J9  -  1  (L  is  the  orbital  angular  momentum  transfer)  spin  transition 
density  is  involved  in  the  (e,e'),  (p,p'),  and  (x,x')  reactions  [Mo-74], 


Using  harmonic  oscillator  forms  for  the  spin  transition  densities  with 
the  size  and  strength  parameters  fixed  from  (e,e'>  and  (p»p')  data,  Carr, 
et  al.  [Ca-83]  obtained  good  agreement  between  DWIA  calculations  and  the 
(w,w')  data  for  the  4“,  17. 79-,  18. 98-,  and  19.80-MeV  states  of  160  and 
the  6~,  T-0  and  T-l  states  at  11.58  and  14.36  MeV  in  28Si. 

The  above  examples  illustrate  that,  in  the  energy  range  of  the 
A(1232)  resonance,  calculations  using  the  DWIA  and  well-known  transition 
densities  are  adequate  for  describing  the  pion-induced  excitation  of 
nuclear  states  that  are  strongly  excited.  However,  for  pion-induced 
excitation  of  nuclear  states  which  are  weakly  excited  by  a  one-step, 
impulse-approximation  mechanism,  nuclear  medium  effects  and  multistep 
processes  may  be  important  and  there  may  not  even  be  qualitative 
agreement  between  pion  inelastic  scattering  data  and  DWIA  calculations. 
This  is  Indeed  the  case  for  the  T-l  member  of  the  weakly  excited  1  + 
doublet,  12.71  MeV  (T-0)  and  15.11  (T-l),  of  l2C  [Mo-82],  As  can  be  seen 
from  Eq.  ( I— 11),  the  ratio  of  cross  sections  for  AS- 1  transitions  should 
be  a(AT-0)/a(AT-l)  -  4/1,  provided  that  the  form  factors  for  the  two 
transitions  are  the  same  [Gu-76].  Furthermore,  this  ratio  is 
approximately  correct  in  DWIA  calculations  which  include  all  s-  and  p- 
wave  partial  amplitudes.  Cohen-Kurath  p-shell  wave  functions  [Co-65] 
describe  the  12.71-  and  15.11-MeV  states  as  near  analogs  (i.e.,  the 
space-spin  pieces  of  the  wave  functions  are  approximately  the  same),  and 
this  description  for  the  spin  densities  is  supported  by  (e,e')  data 
[Co-84].  Nonetheless,  although  both  the  experimental  excitation  function 
and  the  angular  distributions  for  the  12.71-MeV  state  are  well  reproduced 
by  DWIA  calculations,  the  measured  ratio  of  four  times  the  averaged  *  + 


and  »"  differential  cross  sections  for  the  15.11-MeV  state  to  the 
averaged  w+  and  w~  differential  cross  sections  for  the  12,71-MeV  state 
(see  Fig.  1-2)  deviates  significantly  from  one,  especially  at  energies 
near  180  MeV  (A(1232)  dominance),  and  displays  a  rapidly  varying  energy 
dependence.  (Averaging  the  n+  and  w"  differential  cross  sections  removes 
the  effect  of  isospln  mixing  between  the  two  states  on  the  ratio  to 
better  than  It  [Mo-82].)  Also,  DWIA  calculations  do  not  agree  with  the 
15.11-MeV  angular  distributions  at  energies  near  the  A(1232)  resonance 
[Co-84] . 

Uncertainties  in  the  spin  transition  density  or  the  spin- 
dependent  piece  of  the  effective  pion-nucleus  interaction  do  not  seem  to 
be  a  plausible  explanation  [Mo-82]  for  the  anomalous  excitation  function 
for  the  AS«AT»1  transition  to  the  15.11-MeV  level  of  12C.  Rather,  a  more 
likely  explanation  considers  another  process,  in  addition  to  a  one-step, 
impulse-approximation  mechanism,  contributing  to  the  isovector 
transition.  A  possible  mechanism  is  the  previously  proposed  [Mo-82] 
direct  excitation  of  A  particle-nucleon  hole  (A-h)  components  in  the 
final  state  wave  function.  Thus,  to  further  investigate  pion-induced 
excitation  of  nuclear  states  which  are  weakly  excited,  where  the  DWIA  may 
be  inadequate  and  the  direct  excitation  of  A-h  components  may  be  an 
Important  piece  of  the  scattering  amplitude,  this  dissertation  is 
concerned  with  the  AS“AT*i  transition  from  the  1+,  T»0,  ground  state  of 
6Li  to  the  0+,  T”l,  3.563-MeV  level.  As  for  the  15.11-MeV  level  of  12C, 
the  DWIA  analysis  of  this  unnatural-parity  transition  is  particularly 
simple,  involving  only  the  spin-dependent  piece  of  the  effective  pion- 
nucleus  Interaction  and  the  spin  transition  density. 


Fig.  I~2^  RatiOjOf  four  times  the  differential  cross  sections  (averaged 
it+  and  it”)  for  it  inelastic  scattering  to  the  1+,  T-i,  15.11-MeV  state  o£ 
12C  and  the  differential  cross  sections  (averaged  it  and  it  )  for  n“ 
inelastic  scattering  to  the  1+,  T«0,  12.71-MeV  state  of  i2c  [Mo-821. 


Chapter  II  of  this  dissertation  describes  the  experimental 

facility  used  to  acquire  the  data.  Chapter  III  describes  the  extraction 

of  differential  cross  sections  from  the  various  (*+,Tr+')  Q-value  spectra. 

Also  in  this  chapter  is  the  presentation  of  the  data,  which  consist 

of:  (1)  the  first  detailed  measurement  of  a  constant  q,  q  “  109  MeV/c, 

excitation  function  for  it+  inelastic  scattering  to  the  weakly  excited 

3.563-MeV  state  of  6Li  at  incident  pion  energies  from  100  to  260  MeV,  (2) 

partial  angular  distributions  at  T_  -  120  and  180  MeV  for  t+  Inelastic 

w 

scattering  to  the  3.563-MeV  level,  and  (3)  differential  cross  sections 
for  x+  elastic  scattering  and  for  i»+  inelastic  scattering  to  the 
2.185-MeV  (3+,  T*0)  and  4.25-MeV  (2+,  T-0)  states.  A  discussion  of  non- 
relativlstic  multiple-scattering  theory  and  its  application  to  pion- 
nucleus  scattering  is  given  in  Chapter  IV.  Chapter  V  presents  the  DWIA 
calculations  performed  for  the  three  low-lying  excited  states  of  6Li  and 
Includes  discussion  of  the  data  and  calculations.  Appendix  C  is  a  copy 
of  the  paper  resulting  from  this  dissertation  experiment,  which  was 
received  by  Physical  Review  C  in  August  1984. 


II.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


The  data  presented  in  this  dissertation  were  collected  using  the 
Energetic  Pion  Channel  and  Spectrometer  (EPICS)  system  at  the  Clinton  P. 
Anderson  Meson  Physics  Facility  in  Los  Alamos,  New  Mexico  (LAMPF).  LAMPF 
is  described  in  various  Los  Alamos  Scientific  Laboratory  Reports 
[Li-72 ,Al-77 , Li-77] .  Briefly,  LAMPF  consists  of  an  800-MeV  proton  linear 
accelerator,  capable  of  accelerating  simultaneously  H+  ions  and  H~  ions 
(or  polarized  H~  ions),  and  several  experimental  areas  for  nuclear 
structure  and  chemistry  studies.  LAMPF  was  designed  to  produce  a  proton 
beam  with  an  average  current  of  1  mA  at  a  duty  factor  of  6-122.  However, 
during  the  experiment  (June  and  November  1982)  which  provided  the  data 
for  this  dissertation,  the  average  current  was  *  700-750  pA  at  a  7.52 
duty  factor.  The  H+  beam  enters  Experimental  Area  A  and  impinges  on  the 
first  graphite  target  (the  A-l  target),  producing  pions  (and  other 
secondary  particles)  which  the  EPICS  channel  accepts  at  an  angle  of  35° 
from  the  H4’  beam  direction.  Experimental  Area  A  is  shown  in  Fig.  II-l. 

A.  EPICS  system 

The  EPICS  system  [Th-70]  was  designed  to  provide  good  energy 
resolution  and  good  angular  resolution  studies  of  pion-induced  excitation 
of  nuclear  levels,  with  cross  sections  as  low  as  several  nb/sr,  from  low 
incident  pion  energies  of  -  70  MeV  to  energies  above  the  A(1232)  reso¬ 
nance,  -  300  MeV.  EPICS  consists  of  a  high  intensity,  high  resolution 
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Fig.  II— i :  Experimental  Area  A  of  the  LAMPF  facility  [ LI— 7 7 J 


pion  bean  channel,  a  scattering  chamber,  and  a  high  resolution  spec¬ 
trometer.  A  diagram  of  the  EPICS  channel,  scattering  chamber,  and 
spectrometer  is  shown  in  Fig.  II-2.  Table  II— 1  lists  the  specifications 
of  the  EPICS  system. 

The  EPICS  channel,  shown  in  Fig.  I 1-3,  consists  of  four  dipole 
magnets  (BM01-BM04),  three  multipole  focusing  magnets  (FM01-FM03),  and 
four  adjustable  collimating  jaws  (FJ01-FJ04).  The  channel's  length  of 
15.24  m  was  chosen  in  order  to  momentum  select  pions  of  maximum  kinetic 
energy  of  300  MeV,  but  be  short  enough  to  maintain  a  reasonable  flux  of 
70  MeV  pions  even  after  pion  decay  through  the  channel.  The  four  bending 
magnets  provide  charge  and  momentum  selection  and  momentum  dispersion  of 
the  channel  beam.  (In  such  a  dispersed  beam,  the  momentum  of  a  particle 
is  correlated  with  its  position  in  the  beam.)  The  optical  mode  is  point 
to  point  to  point  to  point  in  the  vertical  direction  and  point  to 
parallel  to  point  to  parallel  in  the  horizontal  direction.  FM01-FM03 
each  have  three  windings,  quadrupole,  half-sextupole  (top),  and  half- 
sextupole  (bottom),  and  are  used  for  removing  higher-order  aberrations  of 
the  channel  optics.  The  four  adjustable  collimating  jaws  determine  the 
phase  space  and  flux  of  the  channel  beam.  FJ01  has  only  vertical  jaws 
(disabled),  and  along  with  a  fixed  collimator,  defines  the  channel 
acceptance.  When  FJ01  is  operative,  regulation  of  the  channel  beam 
intensity  and  vertical  divergence  is  accomplished  with  FJ01.  FJ02-FJ04 
each  have  two  pairs  of  jaws,  one  horizontal  and  one  vertical.  The 
vertical  jaws  of  FJ04  define  the  channel  momentum  acceptance  and,  hence, 
the  vertical  size  of  the  channel  beam  at  the  scattering  target,  while  the 
horizontal  jaws  define  the  horizontal  divergence  of  the  beam.  With  all 
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TABLE  II-l:  EPICS  system  specif ications.a 


CHANNEL 


SPECTROMETER 


Solid  angle  3.4  msr 

Ap/p  2% 

Beam  size  (horizontal)  8  cm 

Beam  size  (vertical)  20  cm 

Beam  divergence  (horizontal)  <10  mrad 

Beam  divergence  (vertical)  100  mrad 

Energy  range  70-300  MeV 


Solid  angle  *10  msr 

Ap/p  14% 

Momentum  range  100-750  MeV/c 
Flight  path  >125  cm 

Dispersion  4  cm/X 


*Ref.  [LA-80]. 


TABLE  II-2:  EPICS  channel  beam  pion  flux  and  composition. 


Plon  Fluxa 


Beam  Content 


jaws  fully  open,  Che  channel  beam  at  Che  scaccerlng  cargec  has  a  vercical 
dispersion  of  10  cm/Z,  *  momentum  bice  of  £  1Z,  and  a  divergence  of  100 
mrad.  In  che  horizoncal  plane,  Che  divergence  is  <  10  mrad.  The  resul- 
canc  beam  spoc  is  20  ca  (vercical)  by  8  cm  (horizoncal)  [LA-80],  and  che 
channel  beaa  pion  flux  and  coaposicion  are  presenced  in  Table  II-2. 

Ac  che  cencer  of  che  scaccerlng  chaaber,  locaced  ac  che  focal 
plane  of  che  channel,  is  che  cargec  ladder  concaining  up  Co  four  full- 
size  cargecs  (22.9  cm  x  15.2  ca).  Differenc  targets  are  selecCed  by 
moving  Che  cargec  ladder  vertically.  Also,  Che  cargec  ladder  can  be 
rotated  relative  Co  che  channel  beam.  The  normal  procedure  is  Co  bisect 
the  scaccerlng  angle,  thus  minimizing  the  path  length  traversed  Chrough 
the  cargec.  Mimlmum  beam  travel  through  Che  target  minimizes  straggling 
which,  in  turn,  helps  Co  maximize  che  resolution.  To  balance  acceptable 
resolution  against  count  rates,  target  areal  densities  for  pion  elastic 
and  inelastic  scattering  experiments  range  from  100  to  300  mg/cm2. 

The  EPICS  spectrometer  (see  Fig.  II-4)  consists  of  a  quadrupole 
triplet  (QM01-QM03),  two  dipole  magnets  (BM05-BM06),  and  front  and  rear 
focal  plane  deteccor  systems.  The  spectrometer  rotates  about  its  pivot, 
situated  at  the  focal  plane  of  the  channel,  through  an  angular  range  of 
-10°  to  120°.  The  optical  mode  from  the  scattering  target,  through  the 
quadrupole  triplet,  to  the  front  focal  plane  (located  just  before  the  two 
dipole  magnets)  is  point  to  point  in  the  vertical  direction  and  parallel 
to  point  in  the  horizontal  direction.  A  magnification  of  -1  is  provided 
by  che  quadrupole  triplet,  forming  an  inverted  image  of  the  scattering 
target  in  which  x  is  proportional  to  x  at  the  scattering  target  and  y  is 
proportional  to  the  scattering  angle  at  the  scattering  target.  The  two 
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dipole  magnets  vertically  disperse  the  beam,  4  cm/%,  and  the  useful 
momentum  bite  of  the  spectrometer  is  *  6%.  Their  optical  mode  is  point 
to  point  in  the  vertical  direction  and  point  to  parallel  in  the  hori¬ 
zontal  direction.  Fig.  II-5  illustrates  the  optics  of  the  spectrometer 
and  the  coordinate  system  used  in  defining  particle  trajectories.  Also 
shown  are  the  front  and  rear  focal  plane  detector  systems.  The  front 
focal  plane  detector  system  is  a  set  of  four  multiwire  proportional  drift 
chambers  (F1-F4)  with  F4  located  at  the  focal  plane  of  the  quadrupole 
triplet.  The  rear  focal  plane  detector  system  consists  of  four  multiwire 
proportional  drift  chambers  (R5-R6,  R9-R10),  with  R5  located  at  the  rear 
focal  plane,  a  scintillator  (S2),  a  slab  of  Lucite  used  to  range  out 
protons,  a  scintillator  (S3),  and  a  muon  rejector. 

B.  Detector  system  and  beam  monitoring 

Since  the  spectrometer  Is  not  dispersion  matched  to  the  channel, 
the  detector  systems  at  the  front  and  rear  focal  planes  must  measure  both 
particle  positions  and  angles  in  order  to  determine  the  incident 
particles'  momenta  and  the  scattered  particles'  momenta.  These  two 
quantities  allow  determination  of  the  reaction  kinematics.  Thus,  in 
order  to  obtain  good  momentum  resolution,  the  multiwire  proportional 
drift  chambers  must  provide  good  position  resolution.  The  multiwire  pro¬ 
portional  drift  chambers  used  in  the  EPICS  spectrometer  [ At— 8 1 ,Mo-82a] 
provide  position  resolutions  of  125  urn  (FWHM)  and  may  operate  at  count 
rates  up  to  106  Hz.  The  set  of  four  multiwire  proportional  drift 
chambers  (F1-F4)  at  the  front  focal  plane  are  packaged  in  one  chamber 


Fig.  II-5:  The  EPICS  spectrometer  optics  and  the  detector  system 


assembly.  This  chamber  assembly  consists  of  eight  signal  planes,  two 
orthogonal  sets  of  four  planes  each  in  order  to  obtain  vertical  and  hori¬ 
zontal  position  (x,y)  and  angle  (9,$)  information,  with  each  signal  plane 
separated  by  grounded  foil  planes.  Each  set  of  four  signal  planes  with 
the  signal  wires  in  the  same  direction  are  arranged  as  two  sets  of  two 
signal  planes.  The  two  sets  of  signal  planes  are  separated  by  10  cm. 
Within  each  set  of  two  signal  planes,  the  two  signal  planes  are  separated 
by  1  cm  and  offset  with  respect  to  each  other  by  one-half  of  a  wire 
spacing.  The  set  of  four  multiwire  proportional  drift  chambers  (R5-R6, 
R9-R10)  at  the  rear  focal  plane  are  packaged  in  two  chamber  assemblies, 
each  consisting  of  four  signal  planes,  two  orthogonal  sets  of  two  planes 
each.  Their  construction  is  the  same  as  the  front  chamber  assembly.  All 
signal  planes  are  of  an  alternating  gradient  design,  with  the  anode  wires 
at  positive  high  voltage  (typically  2150  V)  and  a  wire  spacing  of  8  mm, 
and  cathode  wires  at  -200  V  centered  between  the  anode  wires.  The  anode 
wires  are  directly  coupled  into  a  fast  (2.5  ns/cm)  delay  line,  thus 
giving  two  anode  signals  per  signal  plane,  while  alternate  cathode  wires 
are  bussed  together,  thus  giving  two  cathode  signals  per  signal  plane. 
Presently,  only  the  two  anode  signals  per  plane  are  used  to  calculate 
positions  and  angles.  Calculation  of  position  is  as  follows.  A  time 
difference  and  time  sum  from  the  two  anode  signals  are  formed.  The  time 
difference  is  directly  related  to  the  position  of  the  wire  closest  to  the 
ionizing  event,  determining  the  position  of  an  event  to  ±  4  mm.  Further 
position  resolution  is  obtained  from  the  time  sum,  which  is  equal  to 
twice  the  drift  time  to  the  event  wire  plus  a  constant  offset.  However, 


one  must  decide  whether  to  add  to  or  subtract  from  the  drift  distance  the 
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anode  wire  position.  This  left-right  ambiguity  can  be  solved  with  four 
signal  planes,  hence,  the  need  for  eight  signal  planes  at  the  two  focal 
planes,  and  is  described  in  [At-81 ,Mo-82a,Iv-79] .  Angles  are  determined 
from  the  spacing  ( 10  cm)  of  the  two  sets  of  two  signal  planes  which 
comprise  a  set  of  four  signal  planes  whose  anode  and  cathode  wires  are  in 
the  same  direction  in  conjunction  with  the  positions  measured  in  each  set 
of  two  signal  planes.  Thus,  the  information  provided  by  the  multiwire 
proportional  drift  chambers  located  at  the  front  and  rear  focal  planes  of 
the  EPICS  spectrometer  consists  of  the  eight  quantities  x^,  yf ,  9f 
(dxf/dz),  $f  (dyf/dz),  xr,  yr,  8r,  $r. 

Two  scintillators  (S2,  S3)  were  the  only  scintillators  used  as 
part  of  the  EPICS  detector  system  during  this  experiment.  They  are 
located  just  after  the  rear  set  of  multiwire  proportional  drift  chambers 
and  are  separated  by  a  slab  of  Luclte.  Photomultiplier  tubes  are 
optically  connected  to  each  end  of  the  scintillators.  The  purpose  of  the 
scintillators  was  twofold.  They  were  used  as  part  of  the  logic  defining 
a  hardware  trigger  signal  for  events  and  for  particle  identification. 
For  particle  identification,  the  signals  from  S2  and  S3  are  used  to 
measure  a  time  of  flight  between  S2  and  S3  and  pulse  heights,  which  are 
linearly  proportional  (Me-66]  to  the  energy  loss  in  the  scintillators. 
Since  time  of  flight  is  proportional  to  E/pc2  and  energy  loss  is  propor¬ 
tional  to  (ZeE/pc2)2  [Me-66]  (E  is  the  total  energy  of  the  particle,  p  is 
the  momentum,  and  Ze  is  the  charge),  these  two  measurements  aid  in 
distinguishing  particles  of  different  masses  and  charges. 


Following  scintillators  S2  and  S3  is  a  series  of  scintillators 
(S4-S9)  separated  by  carbon  wedges  of  varying  thicknesses.  (The  wedge 
shape  is  needed  to  account  for  the  variation  of  momentum  across  the  rear 
focal  plane.)  This  arrangement  of  six  scintillators  and  six  carbon  wedges 
is  known  as  the  muon  rejector  [Mo-84].  The  muon  rejector  is  designed  to 
eliminate  the  muon  background  from  in-flight  decay  of  pions  before  the 
front  set  of  multiwire  proportional  drift  chambers  and  elastic  scattering 
of  channel  beam  muons.  The  operation  of  the  muon  rejector  is  as  follows. 
Since  muons  and  pions  both  have  a  single  unit  of  charge,  muons  and  pions 
of  the  same  momentum  are  transported  to  the  rear  focal  plane.  For  the 
same  momentum,  the  total  energy  of  a  pion  is  greater  than  a  muon  and  thus 
the  plon's  range  In  matter  is  less  than  that  of  the  muon 
(-dE/dx  «  (ZeE/pc2)2nfi ,  where  ng  is  the  electron  density  of  the  ranging 
material  [Ja-75]).  Furthermore,  pions  also  lose  energy  due  to  nuclear 
collisions.  The  first  carbon  wedge  is  of  the  correct  thickness  to  remove 
pions  of  100  MeV  kinetic  energy  but  not  muons.  The  second  carbon  wedge 
removes  pions  of  140  MeV  kinetic  energy  but  not  muons,  and  so  on.  To 
remove  pions  of  intermediate  energies,  there  is  space  between  S3  and  the 
muon  rejector  to  place  varying  thicknesses  of  aluminum.  The  signals  from 
the  scintillator  following  the  carbon  wedge  of  the  desired  pion  energy 
are  then  used  to  veto  muons  (usually  in  the  software).  During  this 
experiment,  the  muon  rejector  rejected  most  of  the  muon  background 
(90-952),  while  rejecting  only  2-5%  of  the  pions. 

To  obtain  the  normalization  required  in  the  calculation  of  pion 
elastic  and  inelastic  differential  scattering  cross  sections,  three  beam 
monitors  are  used.  Several  beam  monitors  provide  cross  checks  on  the 


stability  of  each  monitor.  An  ionization  chamber  (IC1)  located  down¬ 
stream  of  the  scattering  target  monitors  the  EPICS  channel  beam  current. 
The  primary  beam  current  is  monitored  by  an  ionization  chamber  (BOT) 
located  within  the  plon  production  target  cell  and  a  charge  integrating 
toroidal  coil  (1ACM02)  located  upstream  of  the  plon  production  target. 
Both  the  BOT  and  1ACM02  signals  are  gated  by  a  RUN  gate  so  that  the 
primary  beam  current  is  monitored  only  during  data  acquisition.  BOT  and 
1ACM02  are  reliable  indirect  measures  of  the  EPICS  channel  beam  current 
provided  the  proton  beam  is  always  properly  steered  and  strikes  the  A-l 
target.  For  portions  of  this  experiment,  both  BOT  and  IACM02  were  used 
for  normalization,  because  ICI  partially  blocks  the  spectrometer  entrance 
for  9^3  <  25°  and  is  not  used,  and  BOT  was  not  operational  near  the  end 
of  the  experiment.  A  careful  check  of  the  ratios  IC1/B0T  and  ICI/1ACM02 
( <  2%  fluctuations)  for  >  25°  proved  the  reliability  of  BOT  and 
1ACM02 . 

C.  Data  acquisition  and  analysis  system 

The  data  acquisition  and  analysis  system  for  EPICS  consists  of  an 
on-line  PDP-ll/45  computer  which  runs  under  DEC'S  RSX11-D  operating 
system,  several  CAMAC  crates  containing  CAMAC  electronic  modules,  a 
Microprogrammable  Branch  Driver  (MBD)  interfacing  the  PDP-11/45  and  CAMAC 
crates,  computer  peripherals  (tape  drive,  two  graphics  terminals,  disk 
drive  system,  and  a  printer/plotter),  and  the  LAMPF  standard  data 
acquisition  software  package  "Q"  for  managing  the  data  acquisition  and 
display  of  the  experimental  data  [Am-79J.  The  scintillator  signals  and 


chamber  signals  are  passed  via  coaxial  cables  from  Che  focal  planes  to 
che  councing  house  and  fed  inco  NIM  modules,  discriminaCors ,  Meancimers, 
ADC's,  and  TDC's  for  Che  scinCillaCor  signals  and  discriminaCors  and 
TDC's  for  Che  chamber  signals  (see  Fig.  II-6).  The  logical  "and”  of  che 
meancimed  signals  from  S2  and  S3  and  ac  lease  evo  signals  from  boch  che 
four  x  planes  and  four  y  planes  of  che  fronc  multiwire  propordonal  drift 
chambers  define  a  trigger  signal  (trigger  signal  »  (S2»S3)*(F1  or  F2)*(F3 
or  F4) ) .  This  signal  along  with  a  logical  "and"  of  the  not-busy  signal 
issued  by  the  PDP-11/45  computer  when  it  is  ready  to  accept  data 
constitutes  the  hardware  trigger  (hardware  trigger  ■  (trigger 
signal) *(not-busy  signal)).  The  hardware  trigger  signal  initiates  the 
reading  of  the  data  words  (20-100)  for  that  event.  The  "Q"  software 
[Ke-78]  controls  this  reading,  the  writing  of  the  data  words  for  each 
event  to  magnetic  tape,  and  processing  of  the  events  if  time  is  available 
(MAY  PROCESS  mode  of  on-line  data  acquisition).  Both  the  trigger  signals 
and  the  hardware  trigger  signals  are  scaled  along  with  other  quantities 
such  as  S2 *S3 ,  (FI  or  F2)*(F3  or  F4),  ICl,  BOT,  1ACM02,  etc.  The  "Q" 
software  reads  che  scalers  at  fixed  time  intervals  and  stores  all 
quantities  in  a  scaler  file.  The  ratio  of  the  number  of  hardware  trigger 
signals  (EVENT*BUSY)  to  the  number  of  trigger  signals  (EVENT)  is  a 
measure  of  the  computer's  rate  of  data  taking,  the  computer  live  time 
(CLT  -  EVENT  •‘BUSY/ EVENT). 

An  important  part  of  the  EPICS  electronics  not  shown  in  Fig. 
II-6,  but  used  during  the  experiment,  is  a  fast  clear  circuit  used  to 
eliminate  the  majority  of  the  hardware  trigger  signals  corresponding  to 
elastic  scattering  events  [Se-81a].  The  fast  clear  circuit  aids  in  the 


Fig.  II— 6 :  Block  diagram  of  the  electronics  for  the  EPICS  data 
acquisition  system. 


measurement  of  small  Inelastic  scattering  cross  sections  relative  to 
elastic  scattering.  This  circuit  produces  a  signal  which  is  used  to  veto 
9  out  of  10  or  99  out  of  100  hardware  trigger  signals;  thus,  the  data 
words  for  only  one  of  10  or  one  of  100  elastic  scattering  events  are  read 
into  the  computer,  written  to  magnetic  tape,  and  processed  (if  time  is 
available).  The  number  of  vetoed  hardware  trigger  signals  (FAST  CLEAR) 
are  scaled. 

Once  the  data  words  for  an  event  are  read  into  the  computer,  the 
processing  of  the  data  is  done  on  an  event  by  event  basis  with  the 
standard  analyzer  task  for  EPICS  data-taking  (QEPX)  in  conjunction  with  a 
test  file  (experiment-specific)  and  a  display  package  (experiment- 
specific).  The  goal  of  the  data  processing  is  to  determine  which  events 
are  good  pion  events  and  calculate  the  excitation  energy  given  to  the 
target  nucleus  by  these  pions.  This  goal  is  accomplished  through  the 
main  program  of  the  analyzer  (PROC06).  PROC06  examines  each  event  and 
calls  subroutines  to  calculate  all  chamber  and  scintillator  quantities 
from  the  outputs  of  the  TDC's  and  ADC's,  calculate  trajectory  infor¬ 
mation,  use  the  test  file  to  test  the  various  quantities,  compute  the 
missing  mass  associated  with  each  event,  and  arrange  the  data  in  bins  for 
storage  in  histograms  determined  by  the  display  package. 

The  test  file  [Am-78]  consists  of  the  user-specified  tests  to  be 
performed  on  each  event.  The  test  file  is  structured  into  loops,  and  the 
loops  are  called  separately  and  in  order  by  the  analyzer.  Each  loop  may 
contain  two  types  of  tests,  microtests  and  macrotests.  The  microtests 
are  tests  on  the  data  words,  and  the  macrotests  are  tests  on  low' cal 
combinations  of  previous  microtests  or  macrotests.  Microtests  can  be  a 


BIT  test  which  Is  true  if  a  specified  bit  is  sec  In  a  specified  daca 
word,  a  GATE  (DIRECT  or  INDIRECT)  test  which  is  true  If  the  specified 
integerlzed  data  word  is  greater  than  a  given  lower  limit  and  less  than 
or  equal  to  a  given  upper  limit,  and  a  BOX  test  which  is  a  combination  of 
two  INDIRECT  GATE  tests.  Macrotests  are  logical  "and's",  "inclusive 
or's",  or  "exclusive  or's"  on  specified  microtests  or  macrotests. 

The  display  package  (DSP)  is  used  to  construct  and  display  both 
one-dimensional  histograms  and  two-dimensional  dot  plots.  Histograms  and 
dot  plots  are  defined  by  the  data  words  to  be  hlstogrammed  or  plotted, 
lower  and  upper  limits  of  display,  a  binning  factor,  and  the  tests 
defined  in  the  test  file  which  the  event  must  pass.  Both  histograms  and 
dot  plots  can  be  viewed  on  a  graphics  terminal  during  on-line  data 
acquisition.  However,  only  histograms  can  be  saved  for  later  data 
analysis.  In  addition  to  saving  histograms,  scaled  quantities,  the 
results  of  the  cest  file,  and  the  relevant  kinematic  information  for  each 
experimental  run  (period  of  data  acquisition,  ~  155,000  taped  events)  are 
saved  for  off-line  data  analysis. 

D.  Event  analysis 

The  following  paragraphs  describe  in  more  detail  the  various 
steps  in  the  data  analysis,  leading  to  a  missing  mass  histogram  of  good 
pion  events  which  can  be  used  to  calculate  scattering  cross  sections. 
Fig.  II-7  is  a  test  file  used  during  the  experiment  for  this 
dissertation.  The  first  calculation  for  an  event  performed  by  the 
analyzer  is  the  particle  identification  test  (PID)  (see  test  #11  in  Fig. 
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sALLBUT  IBA 

sALLBUT  IC 

sALLBUT  2C 

sALLBUT  3C 

sALLBUT  4C 

sALLBUT  5C 

sALLBUT  SC 

sALLBUT  9C 

sALLBUT  1BC 


Fig.  XI— 7 :  Sample  test  file 
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II— 7 ) •  The  PID  box  consists  of  two  indirect  gates,  one  on  the  mean 
energy  loss  in  scintillators  S2  and  S3,  and  one  on  the  time  of  flight 
between  S2  and  S3.  The  PID  test  eliminates  protons  but  not  muons.  The 
PID  test  constitutes  the  first  loop  of  the  test  file. 

If  an  event  passes  the  PID  test,  the  analyzer  proceeds  to 
calculate  the  positions  in  the  front  and  rear  multiwire  proportional 
drift  chambers  from  the  time  difference  and  time  sum  of  the  chamber 
signals.  At  this  point,  a  number  of  chamber  tests  are  performed  (second 
loop  of  the  test  file).  The  tests  consist  of  drift-time  tests  and  drift- 
difference  tests.  The  drift-time  tests  are  gates  (100  ns)  on  the  drift 
times  from  the  eight  front  and  eight  rear  chamber  planes  (see  test  #'s 
17-32  in  Fig.  II-7).  The  drift-time  tests  eliminate  events  which  did  not 
produce  signals  in  one  or  more  chamber  planes.  The  drift-difference 
tests  are  applied  to  groups  of  four  chamber  planes,  the  front  x  planes, 
the  front  y  planes,  the  rear  x  planes,  and  the  rear  y  planes.  These 
tests  compare  the  positions  calculated  for  the  four  chamber  planes  to  a 
straight  line  trajectory.  From  this  comparison,  the  drift-difference 
tests  eliminate  events  for  which  two  particles  passed  through  the 
chambers  simultaneously  and  pion  events  which  created  knock-on  electrons 
with  the  group  of  four  chamber  planes.  At  this  stage,  an  event  has 
passed  the  PID  test  and  chamber  tests,  and  the  analyzer  has  computed  x^, 
yf ,  9j,  $£,  xr ,  yf ,  9^, 

The  next  step  is  the  calculation  of  the  quantities  xt  ,  ytgt, 
9tgt»  ♦tgt ,  the  incident  momentum,  the  scattered  momentum,  and  the 
scattering  angle.  The  quantities  at  the  scattering  target  are  calculated 
using  the  measured  positions  and  angles  at  the  front  chambers  in 


conjunction  with  the  transport  matrix  of  the  spectrometer  quadrupoles 
To  first  order 


xtgt  "  “xf  • 


(II-l) 


However,  due  to  abberations  in  the  quadrupole  fields,  within  the  analyzer 
xtgt*  ytgt’  ®tgt*  and  $tgt  are  calculated  as  polynomials  of  up  to  fourth 
order  in  the  front  chamber  positions  and  angles  [Bo-79].  From  xtgt,  one 
can  obtain  the  incident  momentum. 


jr  m  D  •  6 

Xtgt  UC  C 


(I 1-2) 


where  Dc  is  the  dispersion  of  the  channel  (10  cm/Z)  and  5C  is  the  percent 
difference  between  the  Incident  momentum  for  the  event  (Prevent))  and 
the  central  momentum  of  the  channel  (Pc,  determined  from  the  channel 
magnetic  field  settings). 


Prevent)  -  P£ 


(II-3) 


The  scattered  momentum  is  determined,  to  first  order,  from  the  relation 


D  *5 
sp  sp 


( II-4) 


where  D3p  is  the  dispersion  of  the  spectrometer  (4  cm/Z)  and  5sp  is  the 
percent  difference  between  the  scattered  momentum  of  the  event 


(Pf (event))  end  the  centrel  momentum  of  the  spectrometer  (P8p,  determined 
from  the  spectrometer  magnetic  field  settings). 


sp 


Pj (event)  -  Pgp 


sp 


(II-5) 


After  calculating  the  above  quantities,  the  analyzer  calls 
CALKIN,  a  subroutine  to  compute  the  missing  mass  of  the  event.  Missing 
mass  is  the  difference  between  the  mass  of  the  recoil  nucleus  and  the 
mass  of  the  nucleus  before  the  collision  (i.e.,  the  excitation  energy 
given  to  the  nucleus).  Thus, 

Q  -  [(E1(event)  -  Ef( event)  +  M)2  -  P2]1^2  -  M  ,  (II-6) 


where  Q  is  the  missing  mass,  M  is  the  ground  state  mass  of  the  target 
nucleus,  P  is  the  recoil  momentum  of  the  target  nucleus,  and  E^(event) 
(Ej(event))  is  the  total  energy  of  the  incident  (outgoing)  event. 

The  final  step  of  the  analyzer  is  background  rejection  tests  and 
passing  the  binned  data  to  DSP  for  use  in  background  rejection  histograms 
and  missing  mass  histograms.  Four  of  the  background  rejection  tests  are 
gates  on  the  calculated  quantities  xtgt»  ytgt,  9tgt»  and  *tgt* 
Histograms  of  these  quantities  aid  in  determining  the  lower  and  upper 
limits  of  the  gates.  Other  background  rejection  tests  include  gates  on 
the  calculated  quantities  9check  and  *check.  9check  and  *check  are  the 
difference  between  the  angles  measured  at  the  front  chambers  and  the 
angles  calculated  at  the  front  focal  plane  from  the  information  measured 


at  tha  rear  chambers.  Once  again,  histograms  of  0check  and  ♦check.  helP 
to  set  the  gate  limits.  The  9check  and  *check  testa  eliminate  -  99%  of 
the  muons  resulting  from  the  decay  of  pions  after  the  front  chambers. 
However,  these  tests  will  not  eliminate  muons  due  to  pion  decay  before 
the  front  chambers  or  elastically  scattered  channel  beam  muons.  The  muon 
rejector  (muon  rejection  test  #100  in  Fig.  II-7)  is  used  to  eliminate 
3  96%  of  these  muons.  The  missing  mass  histograms  consist  of  the  good 
pion  events  (test  #108),  the  good  events  (test  #93),  and  the  one  out  of 
10  or  100  fast  cleared  good  pion  events  (test  #112).  The  fast  cleared 
missing  mass  histogram  is  used  to  monitor  the  fast  clear  system  in  order 
to  prevent  rejection  of  events  from  excited  states. 


III.  DATA  REDUCTION  AND  RESULTS 


A.  Missing  mass  histograms  and  targets 

Pion  elastic  and  inelastic  scattering  cross  sections  are 
calculated  from  the  peak  areas  (counts)  extracted  from  the  missing  mass 
histograms  of  good  pion  events.  These  histograms  contain  only  those 
events  which  have  passed  test  #108  of  the  test  file  (see  Fig.  II-7). 
Thus,  the  good  pion  events  have  passed  the  particle  identification  tests, 
the  chamber  tests,  and  the  background  rejection  tests.  Each  missing  mass 
histogram  contains  good  pion  events  for  the  full  angular  acceptance  of 
the  spectrometer,  9gpect  ±  ~1.5°,  where  0spect  is  the  scattering  angle 
for  which  the  spectrometer  has  been  set.  Fig.  III-l  shows  a  missing  mass 
histogram  for  6Li(*+,x+')  using  a  205  mg/cm2,  20  cm  x  10  cm  sheet  of 
enriched  lithium  metal  (>  98%  6Li)  for  T^  -  140  MeV  and  9ub  -  26°. 
Although  the  experimental  resolution  was  not  the  best  for  the  205  rag/cm2 
target  (240  keV  full  width  at  half  maximum  (FWHM)),  the  0+,  T*l, 
3.563-MeV  level  is  clearly  resolved  from  the  broad  2+,  T»0,  4.25-MeV 
level  and  the  large  background  due  to  the  continuum  of  final  states  in 
the  alpha-deuteron  breakup  of  6Li.  Also  clearly  evident  in  the  Q-value 
spectrum  is  the  3+,  T*0,  2.185-MeV  level.  There  is  no  indication  of  the 
2+ ,  T-l,  5. 37 -MeV  or  1+,  T-0,  5.65-MeV  levels.  (The  5.65-MeV  level  has 
never  been  observed  in  electron  scattering  [ Be— 82 ] . )  The  discontinuity  in 
the  Q-value  spectrum  between  the  ground  state  and  2.185-MeV  state  is  an 
artifact  of  the  fast  clear  system  used  to  veto  99  out  of  100  events  with 
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a  Q  value  below  «  1.2  MeV.  The  Insert  of  Fig.  III-l  shows  the  same 
missing  mass  histogram  acquired  without  the  use  of  the  fast  clear  system 
and  for  a  shorter  period  of  data  acquisition. 

Four  separate  targets  consisting  of  sheets  of  enriched  lithium 
metal  (>  98Z  6 LI)  fabricated  by  the  Oak  Ridge  National  Laboratory  were 
used  during  the  experiment.  Two  of  the  targets,  with  dimensions  of  22.9 
cm  x  15.2  cm  and  areal  densities  of  202  mg/cm2  and  100  mg/cm2 ,  were  used 
for  approximately  one-half  of  the  total  data  acquisition  and  contained  no 
discernible  contamination.  The  remaining  two  targets,  with  dimensions  of 
20  cm  x  10  cm  and  areal  densities  of  205  mg/cm2  and  95  mg/cm2 ,  were 
contaminated  by  exposure  to  the  air  during  shipping.  (Further 
contamination  was  prevented  by  handling  the  targets  In  a  sealed  glove  bag 
filled  with  argon.  Also,  when  not  In  use,  the  targets  were  stored  In  an 
evacuated  canister.)  Since  lithium  metal  exposed  to  the  air  will  form  the 
compounds  lithium  hydroxide  (L10H)  and  lithium  nitride  (L13N),  possible 
contaminants  were  oxygen,  nitrogen,  and  hydrogen.  The  areal  densities 
for  the  20  cm  x  10  cm  targets  Include  corrections  for  these  contaminants, 
as  described  In  Sec.  III-B. 

Full-size  targets  for  the  target  ladder  are  22.9  cm  x  15.2  cm. 
Two  of  the  targets  (20  cm  x  10  cm)  were  not  full-size  targets  and  care 
was  taken  to  ensure  proper  normalization  of  the  cross  sections.  The  CH2 
targets  used  for  absolute  normalization  (see  Sec.  III-C)  were  of  the  same 
dimensions  as  the  6 LI  targets,  and  both  the  6 Li  and  CHg  targets  were 
mounted  In  the  target  ladder  In  the  same  configuration.  Furthermore,  due 
to  possible  misalignment  of  the  targets  in  the  target  ladder  and  the 
variation  of  the  target  area  Intercepted  by  the  channel  beam  as 


co8^0target^  ^target  18  the  an8^e  between  a  perpendicular  plane  to  the 
bean  and  the  plane  of  the  target ),  @CargeC  was  set  to  bisect  the  central 
angle  of  the  angular  distributions  acquired  with  the  20  cm  x  10  cm 
targets.  For  example,  at  Tff  -  120  MeV  the  thick  and  thin  20  cm  x  10  cm 
targets  were  used  for  9^ab  ■  21°  to  42°.  Thus,  9target  waa  set  t0  16° 
(bisects  the  central  angle  of  32°)  for  each  8lab.  Such  a  procedure 
eliminated  possible  errors  which  could  result  from  differing  target  areas 
intercepted  by  the  beam  as  the  target  angle  varied.  Also,  the  central 
angle  of  the  angular  distribution  was  chosen  to  minimize  straggling  for 
the  set  of  angles.  For  the  two  full-size  targets,  the  target  area  was 
much  greater  than  the  channel  beam  area  and  the  above  precautions  were 
not  needed. 

B.  Extraction  of  peak  areas 

Peak  areas  for  the  elastic  peak  of  6Li,  the  first  three  excited 
states  of  6  Li,  and  the  elastic  peaks  of  the  contaminants  were  extracted 
from  the  6Li(n+,it+' )  Q-value  spectra  using  a  lineshape  oriented  fitting 
program  LOAF  [Sm-78].  LOAF  allows  simultaneous  fitting  of  a  maximum  of 
ten  peaks,  each  with  a  separate  lineshape.  The  peaks  can  be  fit  as 
"free"  peaks  (centroid  of  each  peak  is  allowed  to  move  Independently)  or 
as  "delta"  peaks  (centroids  of  a  group  of  peaks  are  constrained  to  have 
fixed  separations  but  the  position  of  the  group  of  peaks  may  vary). 
Llneshapes  are  taken  from  a  user-specified  peak  in  the  spectrum  to  be 
fitted.  LOAF  automatically  determines  the  background  as  a  polynomial  of 
order  zero  to  ten  (the  polynomial  order  is  specified  by  the  user)  from 


Che  displayed  region  of  Che  histogram,  or  Che  user  may  fix  Che  background 
by  specifying  a  number  of  background  poinCs  which  LOAF  Chen  fiCs  wich  a 
spline  curve.  For  a  displayed  region  of  a  Q-value  spectrum,  LOAF 
deCermines  Che  besC  fit  background  and  then  subtracts  this  best  fit  back¬ 
ground  from  Che  specCrum.  The  remaining  spectrum  is  assumed  Co  be  peaks 
and  is  fitted  with  che  user-specified  lineshapes,  where  only  the  area  and 
centroid  of  each  peak  are  variable  parameters.  A  disadvantage  of  this 
peak  fitting  procedure  is  that  uncertainties  for  peak  areas  extracted 
from  peaks  strongly  correlated  with  the  background  may  be  underestimated. 
Due  to  the  relative  ease  of  use  and  constant  user  Interaction  with  the 
program  LOAF,  LOAF  was  used  for  the  extraction  of  all  peak  areas,  even 
though  the  weakly  excited  3.563-MeV  state  lies  on  top  of  a  large 
continuum  background.  However,  to  ensure  that  LOAF  provided  realistic 
uncertainties  for  the  peak  areas,  one  spectrum  at  each  incident  pion 
energy  was  fitted  with  another  peak  fitting  routine  [Bl-84].  This 
alternate  fitting  program  optimizes  the  background  polynomial  and  peak 
parameters  simultaneously;  thus,  the  peak  area  errors  contain  any 
correlations  between  the  background  and  peak  parameters.  Both  fitting 
programs  returned  consistent  peak  areas  and  errors. 

Due  to  the  differences  in  recoil  energies  of  different  nuclei, 
contaminants  in  the  6Li  targets  appear  as  separate  peaks  in  the 
6Li(»+,7r+')  Q-value  spectra.  The  centroids  of  these  peaks  allow 
determination  of  the  types  of  contaminants,  whereas  the  areas  determine 
the  quantity  and  thus  the  correct  6Li  target  areal  density.  As  a  result 
of  its  light  mass,  the  hydrogen  contamination  was  easy  to  calculate.  For 
both  the  thick  and  thin  20  cm  x  10  cm  targets,  the  hydrogen  peak  was 


fitted  for  ell  spectra  where  It  was  clearly  resolvable  from  the  6 Li 
inelastic  peaks*  These  hydrogen  peak  areas  were  used  in  conjunction  with 
the  tt+  +  p  cross  sections  generated  by  the  computer  program  CROSS 
[Bu-n.d.]  to  calculate  an  average  target  areal  density  due  to  hydrogen. 
The  result  was  <  1%  (by  areal  density)  hydrogen  contamination  for  the 
thick  20  cm  x  10  cm  target  and  *  5 Z  for  the  thin  20  cm  x  10  cm  target. 
The  elastic  peaks  of  16  0  and  14  N  were  not  separable  within  the  elastic 
impurity  peak  appearing  at  the  low  excitation  energy  side  of  the  6 Li 
elastic  peak  (see  Fig.  III-l).  However,  since  in  some  Q-value  spectra 
the  6.13-  and  6.92-MeV  states  of  160  but  no  inelastic  peaks  of  14N  were 
observed,  the  heavy  elastic  impurity  peak  was  assumed  to  be  16 0.  With 
this  assumption,  the  extracted  yields  for  the  heavy  impurity  elastic  peak 
resulted  in  smooth  120-  and  180-MeV  angular  distributions  whose  shapes 
were  in  rough  agreement  with  predicted  x+  +  160  angular  distributions. 
As  for  hydrogen,  the  16 0  elastic  peak  areas  were  used  in  conjunction  with 
theoretical  tt+  +  160  cross  sections  (tested  against  the  data  of  [Ho-80]) 
to  calculate  an  average  target  areal  density  due  to  oxygen.  The  160 
target  areal  density  and  the  hydrogen  target  areal  density  resulted  in 
total  impurities  of  2.3%  and  18%  for  the  thick  and  thin  20  cm  x  10  cm 
targets.  Similar  calculations  assuming  the  heavy  elastic  impurity  peak 
to  be  all  nitrogen  resulted  in  the  same  amount  of  total  impurities  to 
within  t  10%.  In  addition,  a  small  amount  of  nitrogen  contamination 
would  negligibly  affect  the  ®Li  inelastic  cross  sections  (<  1%)  since  all 
the  (x+,tr+')  cross  sections  of  l4N  are  a  factor  of  100  to  1000  less  than 
the  elastic  cross  sections  [Ge-83].  As  a  further  check  on  the 
contaminant  correction  to  the  6Li  target  areal  density,  the  differential 


cross  sections  for  the  elastic  scattering  from  6L1  extracted  from  the 
contaminated  and  uncontaminated  targets  a>.  the  same  Incident  plon  energy 
and  scattering  angle  were  compared  and  found  to  agree  within  2 Z.  Fig. 
III-2  shows  the  fits  from  the  program  LOAF  for  two  spectra  acquired  at 
■  120  MeV  and  0^^  -  27°,  one  with  the  contaminated  thin  20  cm  *  10  cm 
target  and  the  other  with  the  thin  22.9  cm  x  15.2  cm  target.  Due  to  the 
large  amount  of  contamination  in  the  thin  20  cm  x  10  cm  target 
(particularly,  the  broad  hydrogen  peak),  this  target  was  used  only  for 
the  elastic  and  2.185-MeV  states  of  6L1  at  Tff  -  120  MeV  and  0^^  ■  21°, 
27°,  29°,  and  32°,  where  the  oxygen  elastic  was  clearly  resolved  from  the 
6L1  elastic  and  the  oxygen  inelastics  and  hydrogen  elastic  appeared  at 
2.77  MeV  and  above  In  the  Q-value  spectra.  No  7  LI  Impurity  was  observed 
in  the  spectra  for  any  of  the  four  targets. 

To  determine  the  areas  for  the  elastic  peak  of  6 LI  and  the  3+ 
peak,  the  Q-value  spectra  acquired  without  the  fast  clear  system  were 
fitted  in  the  region  between  -1.5  MeV  and  3.0  MeV  with  a  first-order 
polynomial  background  and  two  peaks  constrained  to  be  at  0.0  and  2.185 
MeV  [Aj-84]  (see  Fig.  III-2).  Both  the  elastic  and  3+  peaks  were  fitted 
at  each  incident  plon  energy  with  their  own  llneshapes.  To  obtain  areas 
for  the  0+  and  2+  peaks,  the  Q-value  spectra  acquired  with  the  fast  clear 
system  were  fitted  in  the  region  between  1.4  MeV  and  6.0  MeV  with  a 
third-order  polynomial  background  and  three  peaks  constrained  to  be  at 
2.185  MeV,  3.563  MeV  [Aj-84],  and  4.25  MeV.  The  llneshapes  used  for  the 
2.185-  and  3.563-MeV  peaks  were  extracted  from  the  3+  in  spectra  acquired 
without  the  use  of  the  fast  clear  system.  Llneshapes  were  determined  at 
each  incident  plon  energy  because  the  resolution  varied  with  energy  from 
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Fig.  III-2:  ir+  energy-loss  spectra  taken  at  *  120  MeV  and  8^  -  27°. 
The  top  spectrum  Is  from  the  contaminated  95  mg/cm2 ,  20  cm  >  10  cm 
target.  The  bottom  spectrum  Is  from  the  100  rag/cra2 ,  22.9  cm  x  15.2  cm 
target.  LOAF  was  used  to  fit  the  spectra. 
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200  keV  (220  keV)  FWHM  at  T  -  100  MeV  to  260  keV  (260  keV)  FWHM  at 

n 

Tw  ■  260  MeV  for  the  thin  (thick)  targets.  The  llneshape  used  for  the  2+ 
state  was  constructed  by  folding  the  instrumental  llneshape  taken  from 
the  3*  with  a  Lorentzian  of  680  keV  FWHM.  Both  the  position  and  natural 
width  of  the  2*  level.  4.25  ±  0.02  MeV  and  680  ±  20  keV,  were  determined 
from  fits  to  several  spectra  where  this  state  was  predominant.  Previous 
determinations  from  other  experiments  are  4.27  ±  0.04  MeV  and  690  ±  120 

keV  (6Li(e,e'))  [El-69],  4.29  i  0.02  MeV  and  850  ±  50  keV  (9Be(p,a)) 

[De-83],  4.30+0.1  MeV  and  480+80  keV  (9Be(p,a))  [De-83],  and 

4.3  ±  0.1  MeV  and  600  ±  100  keV  (7Li(3He,a))  [Sc-75].  Fig.  1II-3  shows 
the  fits  to  several  spectra  taken  at  different  incident  pion  energies  and 
scattering  angles. 

The  consistency  of  the  fits  to  the  various  Q-value  spectra  was 
checked  by  extracting  areas  for  100  keV  wide  segments  of  background 
(automatically  determined  by  the  program  LOAF)  centered  about  3,563  and 
4.25  MeV.  Plots  of  the  background  yields  for  both  the  3.563-  and 

4.25-MeV  segments  result  in  smooth  and  continuous  angular  distributions 
at  Tn  -  120  and  180  MeV  (see  Figs.  III-4  and  III-5).  Also,  the  back¬ 
ground  yields  versus  incident  pion  energy  are  smooth  and  continually 


increasing  (see  Figs.  III-6  and  III-7).  Therefore,  there  is  indication 
that  the  extraction  of  the  peak  areas  from  the  large  continuum  background 
is  consistent  for  different  incident  pion  energies  and  scattering  angles. 


However,  this  background  check  does  not  eliminate  the  possibility  of  an 
overall  underestimation  or  overestlmatlon  of  the  background  underneath 
the  3.563-  and  4,25-MeV  peaks. 
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Fig.  III-3:  ir  +  energy-loss  spectra  and  fits  using  the  computer  program 
LOAF. 
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Fig.  III-6:  Background  yields  versus  Incident  plon  energy  for  100  keV 
wide  segments  of  background  centered  about  3.563  MeV.  The  data  points 
correspond  to  a  constant  momentum  transfer  q  »  109  MeV/c. 
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Fig.  III-7:  Background  yields  versus  Incident  pion  energy  for  100  keV 
wide  segments  of  background  centered  about  4.25  MeV.  The  data  points 
correspond  to  a  constant  momentum  transfer  q  =  109  MeV/c. 


C.  Cross  section  computation  and  errors 

The  differential  scattering  cross  section  In  the  center-of-mass 
system  Is  given  by 


do  _  dndt 

dfl  “  FdtNdfl 


(III-l) 


where  dndt  Is  the  number  of  pions  scattered  at  an  angle  8  into  the  solid 
angle  dft  In  the  time  Interval  dt,  Fdt  Is  the  number  of  pions  crossing  a 
unit  area  perpendicular  to  the  channel  beam  In  the  time  interval  dt,  N  is 
the  number  of  Independent  scattering  centers  In  the  target  which  are 
Intercepted  by  the  beam,  and  G  Is  the  Jacobian  of  the  transformation  from 
the  laboratory  to  the  center-of-mass  system.  dndt  Is  given  by  the 
extracted  peak  area  from  a  Q-value  spectrum  corrected  for  the 
lnef f lclences  of  the  EPICS  data  acquisition  system. 

dndt  -  (Peak  Area)*CF  ,  (III-2a) 


where 


CF 


l 


1 


CLT  DREFF 


1  1  „  I  1 

CHEFF  LE  *  SF  SACORR 


(III-2b) 


CLT  is  the  computer  live  time  as  discussed  In  Sec.  II-C.  DREFF  is  the 
drift  efficiency  of  the  raultiwire  proportional  drift  chambers  and 
corrects  for  those  events  producing  signals  in  all  chamber  planes  but  not 


passing  the  drift-difference  tests.  CHEFF  is  the  chamber  efficiency  and 
Is  the  product  of  the  sixteen  Individual  chamber  efficiencies,  each  given 
by 


CHEFFi 


All  Chambers  OK 
All  Chambers  but  1  OK 


(III-3) 


LE  is  the  fraction  of  analyzed  events.  For  on-line  data  acquisition  In 
the  MAT  PROCESS  mode,  LE  <  1.0.  However,  all  data  tapes  are  replayed  In 
the  MUST  PROCESS  mode  so  that  LE  -  1.0  provided  there  are  no  Input/output 
tape  errors.  SF  corrects  for  the  plons  that  decay  after  the  scattering 
target  and  before  scintillator  S3. 

SF  -  exp(-t/r ' )  ,  (III-4a) 


where 


t 


(III-Ab) 


(III-4c) 


is  the  relativistic  momentum  of  the  plon,  mw  is  the  rest  mass  of  the 
plon,  tw  Is  the  mean  lifetime  of  the  plon,  y  ■  l/(l-v2/c2  and  L  Is 

the  path  length  from  the  scattering  target  to  S3  which  Is  corrected  to 
first  order  In  5gp  for  varying  path  lengths,  L  -  (12.28  -  0.035»5gp) 
meters.  SACORR  is  a  correction  factor  for  the  variation  of  the 
spectrometer's  solid  angle  as  a  function  of  This  correction  factor 


53 


Is  determined  by  measuring  the  yield  of  a  state  (usually  an  elastic  peak) 
for  -6?  <  6gp  <  6X.  5gp  is  varied  by  varying  the  spectrometer  field 
settings.  For  all  of  the  experimental  runs,  the  spectrometer  quadrupole 
and  dipole  magnetic  fields  were  set  to  values  such  chat  pions  exciting 
the  3.563-MeV  level  had  the  same  6gp,  where  6gp  was  chosen  to  correspond 
to  the  flat  region  of  the  spectrometer's  solid  angle  variation  versus  6gp 
curve  (acceptance  curve).  This  procedure  minimized  solid  angle 
correction  errors.  Fdt  is  given  by  the  relation 


Monitor  *cos( 0t  t) 
Fdt  «  _ car  gee 


(1II-5) 


where  Monitor  indicates  the  number  of  scaled  XC1,  BOT,  or  1ACM02  pulses, 
cos(  ®target^  corrects  for  the  scattering  target  not  being  perpendicular 
to  the  channel  beam,  a  is  the  target  area  intercepted  by  the  beam,  and  £ 
is  a  proportionality  factor  since  the  beam  monitors  are  only  relative 
monitors  of  the  pion  current.  Lastly, 

N  -  a*n*d  ,  (IIl-6a) 


with  a  Che  target  area  intercepted  by  the  beam,  d  the  target  thickness, 
and  n  the  number  of  scattering  centers  per  unit  volume, 

No' 

n  ■  _ 

A 


(III-6b) 
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NQ  Is  Avagadro's  number,  p  is  Che  target  density,  and  A  is  the  mass  of 
the  target  in  graas/mole.  Combining  Eqs.  C III— 1  to  III-6), 

do  _  C«(Peak  Area) *NORM*A  > 

dfl  N0*areal  density*€*dft  *  ' 


where  NORM*  CF/(Monitor *cos( ®targec))  anc*  ctie  area^  density  is  d*p. 
Since  neither  Z  nor  dfi  are  known,  the  normal  procedure  is  to  determine 
Che  absolute  magnitude  of  the  differential  scattering  cross  section  by 
measuring  a  known  differential  scattering  cross  section.  Thus,  for  each 
incident  plon  energy,  spectra  were  accumulated  for  w+  +  p  scattering  with 
CH2  targets  of  the  same  dimensions  as  the  6L1  targets.  Using  the 
Coulomb-corrected  phase-shift  predictions  from  the  computer  code  CROSS 
[Bu-n.d.],  which  uses  the  phase  shifts  of  Rowe,  Salomon,  and  Landau 
[Ro-78],  to  determine  da/dilH,  from  Eq.  ( II I— 7 ) 


5dfl 


(^•(Peak  Area)H*NORMH*AH 
'Jq* (areal  density )CH^ *do/dftH 


14 

T 


(III-8) 


where  the  14/2  corrects  for  one  mole  of  CH2  containing  two  grams  of 
hydrogen  and  twelve  grams  of  carbon.  Inserting  Eq.  (III-8)  into  Eq. 
(III-7)  yields  the  expression  for  the  absolute  differential  scattering 
cross  section 


rdo)c  _  G6U  .  <Peak  Ar«a>6L1  ^  (N0RM)6Li 
dlT  6LI  Gh  (Peak  Area)H  *  (N0RM)H 


w-. . 


.V.-. 


*  *  •  ' 


x 


X 


(areal  densicy)^^ 
(areal  densicy)6Li 
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(HI-9) 


Eq.  (III-9)  reveals  che  possible  sources  of  error  in  the 
differencial  scaccering  cross  seccions.  The  error  bars  presenced  in  che 
figures  and  Cable  are  errors  due  only  Co  che  scaciscical  and  fleeing 
errors  in  che  excracced  peak  areas  escimaced  by  che  compuCer  program 
LOAF.  Addidonal  Cocal  uncerCalnCies  are  escimaced  Co  be  3  t  72  due  eo 
uncercalncies  of  ±  3%  in  chamber  efficiencies,  1  32  in  pion  survival 
fraccion  (from  che  pach  lengeh  correccions) ,  £  22  in  che  speccromecer's 
solid  angle  variacion  wich  momenCum  in  che  speccromecer  (probably  over- 
esdmaCed  since  che  pion  evencs  of  Che  elasCic  and  inelascic  peaks  in  Che 
6L1  Q-value  specCra  had  6gp's  corresponding  co  Che  flac  region  of  Che 
speccromecer's  accepcance  curve),  £  32  in  channel  beam  monicoring,  and 
±  32  in  che  normalizaCion  Co  ir+  elasCic  scaccering  from  hydrogen. 
Furchermore,  che  daca  for  che  3.563-  and  4.25-MeV  sCaCes  concain 
syscemacic  errors  of  £  152  and  t  102,  respeccively ,  due  Co  che 
uncercaincy  i.  che  fleeing  of  che  large  concinuum  background  and  che 
uncercaincy  in  che  posicion  and  widch  of  Che  2+  scace.  These  syscemacic 
errors  were  inferred  from  Che  fleeing  program  LOAF  by  varying  che  order 
of  che  polynomial  fic  co  che  background,  and  by  varying  che  posicion  and 
widch  of  che  2+  scace  from  4.23  Co  4.27  MeV  and  660  Co  700  keV. 
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0.  Data 


The  data  are  presented  in  Tables  A-l  through  A-4  In  Appendix  A 
and  in  Figs.  III-8  through  III-15.  For  all  figures,  the  results  of 

multiple  measurements  at  a  given  pion  energy  and  scattering  angle  are 
plotted  as  a  single  data  point  which  is  the  weighted  average  of  the 

multiple  measurements.  The  incident  pion  energies  and  scattering  angles 
at  which  the  data  were  acquired  are  the  result  of  the  main  goal  of  the 
experiment — the  measurement  of  a  detailed  constant-q  excitation  function 
for  x+  inelastic  scattering  to  the  3.563-MeV  state  of  6Li.  The  incident 
pion  energies  range  from  100  MeV  to  260  MeV,  thus  covering  the  lower 
energy  region  through  the  upper  energy  region  of  the  A(1232)  resonance. 
Tn  *  100  MeV  is  also  the  lowest  pion  energy  normally  used  for  most 

experiments  at  EPICS  due  to  the  decreased  incident  pion  flux,  decreased 
pion  survival  fraction  in  the  spectrometer  and  experimental  time 
limitations.  T^  •  260  MeV  is  also  not  the  highest  usable  pion  energy, 
but  higher  energies  were  not  feasible  for  this  experiment  because  of  the 
rapidly  decreasing  scattering  cross  section  for  the  3.563-MeV  state  at 
T^  >  230  MeV.  Except  for  T^  *  120  and  180  MeV,  there  are  data  at  only 
one  scattering  angle  for  each  incident  pion  energy.  The  scattering 

angles  correspond  to  a  constant  momentum  transfer  of  q  =  109  MeV/c,  with 
an  error  In  q  of  <  7%.  This  momentum  transfer  was  determined  during  on¬ 
line  data  acquisition  from  the  apparent  maximum  of  the  angular 
distribution  at  Tff  ■  120  MeV  for  ir+  inelastic  scattering  to  the  3.563-MeV 
level.  However,  after  subsequent  off-line  data  analysis,  the  maximum  of 
the  120  MeV  angular  distribution  occurs  at  8lafa  »  27°  (the  midpoint  of 
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Fig.  III-9:  Differential  cross  sections  for  n+  elastic  scattering  for 
6Li  at  a  constant  q  *  109  MeV/c. 
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Che  three  highest  data  points  in  Fig.  III-14),  yielding  a  q  »  101  MeV/c. 
No  corrections  have  been  made  to  the  data  for  the  variations  in  q.  Such 
corrections  result  in  small  changes  in  the  differential  scattering  cross 
sections  for  the  3.563-MeV  state,  but  not  for  the  elastic  and  2.185-  and 
4.25-MeV  states  (discussed  below).  T  »  120  and  180  MeV  were  chosen  as 

II 

the  incident  pion  energies  for  the  partial  angular  distributions  in  order 
to  have  angular  distributions  off  and  on  the  peak  of  the  4(1232) 
resonance. 

Fig.  III-8  shows  the  differential  cross  sections  for  w+  elastic 
scattering  from  6Li  at  •  120  and  180  MeV.  Fig.  1II-9  shows 
differencial  cross  sections  for  100  <  <  260  MeV  and  q  *  109  MeV/c. 
The  excitation  function  continually  increases  with  increasing  Incident 
pion  energy,  a  characteristic  feature  for  elastic  scattering  of  a 
strongly  absorbed  probe  [Bl-66],  The  dip  in  the  data  curve  at  T^  •  180 
MeV  is  partially  a  result  of  the  momentum  transfer  being  approximately  5? 
too  high  (q  *  114  MeV/c).  Extrapolation  to  q  =  109  MeV/c  moves  the 
Tw  -  180  MeV  data  point  from  80  mb/sr  to  90  mb/sr. 

The  differential  cross  sections  for  w+  inelastic  scattering  to 
the  3+,  T-0,  2.185-MeV  level  at  *  120  and  180  MeV  are  shown  in  Fig. 
111-10.  Within  the  angular  range  measured,  the  angular  distributions 
have  the  shape  expected  for  a  strong  transition  in  the  lp  shell  in  which 
the  dominant  transition  amplitude  is  for  J-2,  L-2,  S-0  (J  is  the  total 
angular  momentum  transfer,  L  is  the  orbital  angular  momentum  transfer, 
and  S  is  the  spin  angular  momentum  transfer)  [Le-80],  Fig.  I II— 1 1 
presents  the  differential  cross  sections  for  100  <  T.  <  260  MeV  and 

II 

q  *  109  MeV/c.  The  excitation  function  continually  increases  with 


0  10  20  30  40  50  60 
Oem  (deg) 


Fig.  III-10:  Angular  distributions  for  n+  inelastic  scattering  to  the 
3+,  T-0,  2.185-MeV  state  of  6L1  for  T  -  120  and  180  MeV. 
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Fig.  X II— l l :  Differential  cross  sections  for  w  +  Inelastic  scattering  to 
the  3+,  T-0,  2.185-WeV  state  of  6Li  at  a  constant  q  =  109  MeV/c. 
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Increasing  incident  plon  energy  as  expected  for  this  natural-parity, 
AS*0,  transition  (see  Chap.  1). 

Fig.  III-12  shows  the  differential  cross  sections  for  x+ 
inelastic  scattering  to  the  2*,  T-0 ,  4.25-MeV  level  at  Tt  «  120  and  180 
MeV.  The  shape  of  the  angular  distributions  are  the  saae  as  for  the 
inelastic  transition  to  the  2.185-MeV  state,  consistent  with  the 
predictions  of  Lee  and  Kurath  [Le-80].  This  is  also  consistent  with  the 
electron  scattering  fora  factors  for  both  natural-parity  transitions  to 
the  2.185-  and  4.25-MeV  levels.  The  fora  factors  have  the  sane 
dependences  on  the  momentum  transfer  [Be-63],  Fig.  III-13  presents  the 
excitation  function  taken  at  q  *  109  MeV/c.  As  for  the  2.185-MeV  state, 
the  energy  dependence  is  an  increasing  function  of  the  incident  pion 
energy.  The  fluctuations  in  the  data  curve  are  partially  the  result  of 
variations  in  q.  For  example,  the  Tv  -  180  MeV  data  point  corresponds  to 
q  3  113  MeV/c.  Extrapolation  to  q  3  109  MeV/c  moves  the  data  point  from 
0.7  mb/sr  to  0.6  mb/sr.  Also,  the  fluctuations  are  not  as  great  as  they 
appear  in  Fig.  III-13  when  one  accounts  for  the  uncertainties  discussed 
in  Sec.  III-C. 

The  differential  cross  sections  for  x+  inelastic  scattering  to 
the  0+,  T»l,  3.563-MeV  state  are  presented  in  Figs.  III-14  and  111-15. 
The  120-MeV  angular  distribution  peaks  at  an  angle  (6  =  28°)  smaller 
than  the  peaks  *  45°)  of  the  120-MeV  angular  distributions  for  the 
2.185-  and  4.25-MeV  states.  Such  a  result  is  expected  for  this  J«l,  L*0, 
S“1  transition  [Le-80].  The  excitation  function  increases  from  T^  -  100 
to  190  MeV  and  then  decreases.  A  direct  comparison  to  the  excitation 
function  for  w*  inelastic  scattering  to  the  15.11-MeV  state  of  12C  or  to 
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Fig,  XII— 1 3 ;  Differential  cross  sections  for  n+  inelastic  scattering  to 
the  2+,  T-0,  4.25-MeV  state  of  6Li  at  q  *  109  MeV/c. 
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Fig.  III-15:  Excitation  function  at  a  constant  q  »  109  MeV/c  for  w 
inelastic  scattering  to  the  0+,  T-l,  3.563-MeV  state  of  6Li. 
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sin28  (see  Chap.  1)  is  not  meaningful  because  of  distortion  effects  of 
the-  incident  plan  wave  function.  Further  discussion  of  the  excitation 
function  is  given  in  Chap.  V. 


IV.  THEORY 


A  proper  understanding  of  pion-nucleus  (rrA)  scattering  data  would 
require  theoretical  methods  incorporating  relativistic  scattering 
techniques,  since  nA  interactions  Involve  both  the  absorption  and 
production  of  plons.  However,  at  present,  there  are  almost  no  systematic 
relativistic  theories  available.  The  best  theoretical  tool  is  the  non- 
relativistic  multiple-scattering  formalism  develcyed  by  Watson  and  others 
in  the  1950's  [Fr-53,Fr-56,Wa-57 ,Wa-58] ,  along  with  corrections  for 
relativistic  kinematics.  The  nonrelativistic  multiple-scattering 
formalism  leads  quite  clearly  to  the  lowest-order  result  for  the  optical 
potential  for  irA  elastic  scattering  and  the  transition  amplitude  for  ttA 
inelastic  scattering.  (An  optical  potential  is  a  potential  Involving  the 
projectile  (pion)  variables  only  and,  when  used  in  a  one-particle 
Schrtfdlnger  equation  for  the  projectile  (Klein-Gordon  equation  for  the 
pion),  yields  the  solution  to  the  multiple-scattering  problem  for 
projectile-nucleus  (irA)  elastic  scattering.)  The  lowest-order  results  for 
the  optical  potential  and  transition  amplitude  involve  the  basic  pion- 
nucleon  (irN)  T  matrix  and  the  nuclear  ground  state  and  transition 
densities,  respectively,  and  have  been  successful  as  lowest-order 
approximations  to  the  physics  of  wA  elastic  and  inelastic  scattering  in 
the  medium-energy  range  (100  <  Tff  <  300  MeV). 
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A.  Nonrelatlvlstic  multiple-scattering  formalism 


is 


For  a  pion  incident  on  a  nucleus  of  A  nucleons ,  the  Hamiltonian 


H  -  H0(rltr2,..«,rA;r)  +  V(r1,t2,  —  ,rA;?)  , 


(IV- la) 


with 


Ho  “  HN<rl*r2- 


K,(?) 


V  “  * 


(IV-lb) 


Hn  Is  the  nuclear  Hamiltonian  for  A  nucleons  including  the  nuclear  recoil 
(spin  and  isospln  degrees  of  freedom  are  not  shown  explicitly),  1C,,  is  the 
kinetic  energy  operator  for  the  pion,  and  vt  is  the  potential  interaction 
operator  between  the  pion  and  the  ith  nucleon.  The  variables  fi****,^ 
refer  to  the  nucleon  coordinates  and  r  refers  to  the  pion  coordinate. 
The  Llppmann-Schwinger  equation  for  the  irA  transition  operator  ,T,  is 

T  -  V  +  VG.T  -  t  v.  +  Z  v.G_T  ,  (IV-2a) 

°  l-l  1  i-1  1  0 

where 


* 


E  -  HN  - 


+ 


in 


(IV-2b) 


GQ  is  the  free  Green's  operator  of  the  Hamiltonian  HQ  and  E  is  the 
collision  energy  in  the  itA  center-of-mass  system, 


E  -  (k2  +  m2)1/2  +  +  M2)1/2  +  E”ucl 


(IV-3) 


|  kQ  Is  Che  pion-nucleus  relative  momentum,  M  is  the  mass  of  Che  nucleus 

(M  *  Am) ,  mw  is  the  pion  mass,  and  EqUc1  is  the  nuclear  ground  state 
energy  (taken  to  be  zero,  see  Eq.  (IV-21)).  T  as  defined  in  Eq.  (IV-2) 
is  a  many-body  operator  involving  both  nucleon  coordinates  and  pion 
coordinates.  Since  the  free  xN  T  matrix  is  well  parameterized  from 
experimental  data  (see  Chap.  I),  the  objective  of  multiple-scattering 
theory  is  to  express  the  xA  T  matrix  in  terms  of  the  free  xN  T  matrix. 
This  objective  is  accomplished  by  first  rewriting  Eq.  (IV-2)  in  a  form 
which  allows  the  exact  summing  of  all  the  pion  rescatterings  from  a  given 
nucleon  within  the  nucleus  [Ei-80] — 

I 


T 


A 

■  Z  T. 
i-l  1 


(IV-4a) 


with  Tj  defined  as 


Ti  =  vi  +  viGoT  •  ( IV-4b) 

Using  Eq.  (IV-4a)  in  Eq.  (IV-4b), 

n  •  •  <IV-5) 

Subtract  from  both  sides  of  Eq.  ( IV— 5 )  the  quantity  VjG0T^,  thus 


obtaining 
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O  "  vi.Go)Ti  -  vj  +  ViG^S^Tj  . 


(IV-6) 


Dividing  both  sides  of  Eq.  (IV-6)  by  ( 1  -  VjGg)  and  defining 


Ti  5(1-  v^)"^  , 


Eq.  (IV-6)  becomes 


T,  -  T.  +  t.G„  £  T. 

1  1  1  °j*l  J 


(IV-7) 


(IV-8) 


Substituting  this  result  into  Eq.  (IV-4a),  the  Lippman-Schwinger  equation 
for  the  irA  transition  operator  is 


A 

l 

1-1 


T1  + 


A 

Z  T.G  Z  T. 
1-1  1  °j*i  J 


(IV-9) 


Iterating  Eq.  (IV-9)  yields  the  Watson  series  for  the  irA  (more  generally, 
projectile-nucleus)  transition  operator — 


A 

m  Z  Xj 
i-1  1 


A 

Z  t,G„  Z 
l-l  1  °j*l 


A 

T,  +  Z 
J  i-1 


TtGQ  E  T.G.  Z  Tk  + 
1  °j*i  J  °k*j  * 


(IV-10) 


Eq.  (IV-10)  Is  a  general  result  for  any  form  of  the  Green's  operator  ,GQ, 
and  can  be  applied  to  Inelastic  as  well  as  to  elastic  scattering 
processes.  The  physics  of  Eq.  (IV-10)  is  clear  after  interpretation  of 
the  T's.  Multiplying  by  (1  -  VjG0)  and  adding  v^Gqt ^  to  both  sides  of 
Eq.  (IV-7)  gives 


viVi 


(IV-11) 
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Eq.  (IV-11)  Is  the  Llppmann-Schwlnger  equation  for  the  scattering  of  a 
plon  from  a  nucleon  which  Is  In  Interaction  with  the  other  A-l  nucleons 
of  the  nucleus.  t  la  a  complicated  many-body  operator  Involving  the 
nuclear  Hamiltonian  Returning  to  Eq.  (IV-10),  one  sees  that  irA 
scattering  Is  expressed  as  an  Infinite  series  of  wN  scattering  terms  (In 
the  presence  of  other  nucleons)  for  single,  double,  triple,  etc. 
scatterings,  with  successive  scatterings  on  the  same  nucleon  excluded  to 
avoid  double  counting. 

To  solve  Eq.  (IV-10)  for  the  irA  T  matrix,  one  must  be  able  to 
obtain  an  expression  for  the  t's.  This  is  a  very  difficult  problem  since 
the  t's  involve  the  nuclear  Hamiltonian  H^.  One  can  relate  the  x^  to  the 
free  irN  transition  operator,  t., 

tt  -  vt  +  ,  (IV-12a) 


where 


-  Kn  - 


K«  +  in 


(IV-12b) 


Kn  Is  the  kinetic  energy  operator  for  the  nucleon  and  u  1$  a  suitable 
choice  for  the  irN  collision  energy  in  the  irA  center-of-mass  system.  The 
relation  between  and  t^  Is 

Tl  *  h  +  Cl(G0  -  gQ)ri  .  (IV-13) 


With  this  expression  for  x  in  Eq.  (IV-10),  the  irA  transition  operator  Is 
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expressed  In  terms  of  the  free  xN  transition  operators  plus  nuclear 
medium  corrections.  If  the  difference  between  GQ  and  g0  is  small,  one 
can  approximate  »  t^  (impulse  approximation)  and  the  Watson  series  for 
the  irA  T  matrix  is  greatly  simplified,  since  the  itA  T  matrix  then 
Involves  the  experimentally  well-determined  nN  phase  shifts.  The  Impulse 
approximation  is  based  on  the  assumption  that  the  energy  of  the  plon  is 
sufficiently  high  compared  to  the  binding  energy  of  a  single  nucleon,  and 
thus  one  may  ignore  the  nuclear  medium.  However,  the  Impulse  approxi¬ 
mation  does  not  Include  ignoring  the  momentum  of  the  bound  nucleon.  With 
*  tj,  Eq.  (IV-9)  becomes 

A  A 

T  *  It,  +  Z  t . G_  Z  T.  .  (IV-14) 

1-1  i  i-l  l  Oj*i  j 


Eq.  (IV-14)  Is  the  starting  point  for  deriving  the  lowest-order  optical 
potential  for  irA  elastic  scattering  and  also  the  transition  amplitude  for 
irA  Inelastic  scattering. 


B.  Elastic  scattering 

Eq.  (IV-14)  restricts  the  summation  of  scattering  terms  such  that 
there  are  no  successive  scatterings  on  the  same  nucleon.  If  this 
restriction  is  dropped  (large  nucleus  approximation),  then  one  extra  term 
is  included,  introducing  an  error  of  order  l/A,  which  for  large  nuclei  is 
small  [St-74],  With  this  approximation,  Eq.  (IV-14)  is  replaced  by 

A  A 

T  ■  I  t,  +  Z  t.G-ET, 
i-i  i  i-i  i  Qj  i 


(IV-15) 
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Rewriting, 

T  •  T  +  TG0T  ,  (IV-16a) 

where 

T  -  It.  .  (IV- 16b) 

i-1  1 

For  elastic  scattering,  in  which  the  nucleus  remains  in  its  ground  state, 
one  must  take  nuclear  ground  state  matrix  elements  of  the  *A  transition 
operator.  For  this  purpose,  introduce 

|0>  -  •  l«>  *  l*o(^l»^2****’^A>>  (IV-17) 

for  the  nuclear  ground  state  and  excited  state  wave  functions, 
respectively.  The  nuclear  ground  state  matrix  element  of  T  is 

<0|T|0>  -  <0 | T | 0>  +  <0|TGoT|0> 

-  <0 1 T (0>  +  EE,<0|T|e'Xe'|Go|eXe|T|0>  .  (IV-18) 

Noticing  that 

<e'|G0|e>  -  - _i - <e'|e>  -  G0(e)6£.e  ,  (IV-19) 

E  -  E£uci  -  K„  +  in 


Eq.  (IV-18)  can  be  written 
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<0|T|0>  -  <0 | T | 0>  +  I<0|T|e>Go(e)<e(T|0>  .  (IV-20) 

In  Che  second  term  of  Eq.  (IV-20),  there  are  both  diagonal,  e-0,  and 
nondiagonal,  e*0,  matrix  elements.  The  diagonal  matrix  elements  should 
be  larger  than  the  nondiagonal  matrix  elements  due  to  the  expected  good 
overlap  of  the  nuclear  ground  state  wave  function  with  itself  (especially 
for  small  momentum  transfers)  (Ei-80).  The  diagonal  terms  are  usually 
called  the  coherent  terms  since  they  correspond  to  the  elastic  case  where 
the  pion  wavelets  from  different  nucleon  scattering  centers  add 
coherently.  Exhibiting  the  coherent  terms  explicitly,  Eq.  (IV-20)  is 

<0 1  T  |  0>  -  <0 1  T 1 0>  +  <0 1  T 1 0> - l- - <0 1 T I  0>  + 

E  -  k„  +  in  1  1 

[  E  <0 1 T | e>G  ( e)<e 1 T J 0> ]  ,  (IV-21) 

e*0  ° 

where  the  term  in  brackets  is  small  and  the  arbitrary  zero  of  energy  has 
been  determined  by  letting  the  nuclear  ground  state  energy  vanish, 
E-Cl  -0.  Taking  only  the  coherent  terms  (the  coherent  approximation), 

<0|T|0>  -  <0 1 T 1 0>  +  <0| T|0>  ■_  1 - - — <0 1 T 1  0>  .  (IV-22) 

E  "  K.  ^  in 


( IV-23a) 


where  pI*****Pa»  Pi****»Pa  ar*  the  momenta  of  the  individual  nucleons 
Eq.  (IV-24)  can  be  written  as 


<^'|vop|ic> 


A  fd3p1  d3p2 

1-1  (2w!)3(2ir)3 


d3pj  d3p£  d3 


PA 


(2wi3(2ir)i  TZrj* 


*o<PiP2,*,Pi*,*Pa>  * 


| t1(u) | iTp1>*0( pj^* •  'Pi* •  *pA)  . 


(IV-26) 


<k'p^| t^(u)|kpj>  contains  a  momentum  conserving  delta  function.  Thus, 
define  a  reduced  irN  T  matrix  as  [Er-80l 


I tl(<u) jkp±>  -  I ti<0»)|  |^p1>(2ir)36Cp£-pi+q)  ,  (IV-27) 


where  q  *  k'  -  k  is  the  three-momentum  transferred  to  the  pion. 
Substituting  Eq.  (IV-27)  Into  Eq.  (IV-26),  the  plon  momentum  space 
representation  of  the  optical  potential  is 


<k* I vop I ^  ”^nr7y<^~^  lt(u)l  lkP>P(P-q.P)  » 


( IV-28a) 


where  p(p-q,p)  is  the  single-nucleon  density  and  is  given  by 


^  .  £  rd  pi  d3p2  ...d3Pi-id3Pi+i  ..  u  ha 
(p-q,p  JlJ(2w)3(2x)1*  T2X)3  (2w)3  '  ’u^O7 


d3 


Pa 


.  a-i+  ?  ».*,+  +  ♦  ♦  +♦ 

6(p  +  E  Pj  +  k0)®0(P|P2***Pi-i  P_<1  Pi+I  *  *  *  PA) 
*0(PlP2***Pi-lPPi+l*"PA>  • 


( IV-28b) 
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The  single-nucleon  density  contains  a  momentum  conserving  delta  function 
which  ensures  that  the  nucleus  has  a  total  momentum  of  -lcQ  in  the  irA 
center-of-mass  system  (see  Eq.  (IV-3)).  <£'p-q [ | t(w) | | kp>  possesses  both 

spin  and  isospln  degrees  of  freedom  and  may  be  written  in  terms  of  spin 
and  isospin  operators  [  La-73,  Ei-80]~ 

<ic'p-q  1 1 1 ( u>)  |  | £p>  -  <2'p-q|  |  t8S(w)|  |itp>  + 

<£'p-q  1 1  tv8(  til)  j  |&p>o*ft  + 

<k'p-q | | t8V(u)| |kp>I«t  + 

<k'p-q  j  |  tvv(<u)|  |kp>o*fiI»r  ■ 

E<k'p-q||t£(u)| |kp>05  ,  (IV-29) 

*  *  * 

where  a,  d,  I,  and  t  are  defined  in  Chap.  I.  Eq.  (IV-28)  may  then  be 
rewritten  as 

<k^|V0p|£>  *  AZ/^— yj.<S'p-q  |  1 1  |  | kp>p^(p-q  ,p)  ,  (IV-30a) 


with 


♦  P 2  ^  pa  ♦  a  ♦  «► 

Pg<P*^'P)  -  *  ji/j  *  ko>°E 


♦0( P-9P2  *  * 'Pa5  V PP2 ‘ *  PaJ 


(IV-30b) 


The  p^'s  are  the  nuclear  ground  state  spin/isospin  scalar  and  vector 
densities.  Compared  to  <k'p-q| | t^(w)| |£p>,  p^(p-q.p)  is  a  sharply  peaked 
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function  of  the  nucleon  momentum  p.  This  Is  a  result  of  the  nuclear  size 
being  ouch  larger  than  the  range  of  the  *N  Interaction.  Thus,  one  may 
use  Che  factorization  approximation  and  remove  <£'p-q | | t^(w)| |£p>  from 
the  Integral  In  Eq.  (IV-30a),  after  evaluating  at  some  average  nucleon 
momentum  p0.  Eq.  (IV-30)  becomes 

<*'|Vop|£>  -  Al<2'pQ-q| | t5(w0) | |fcp0>p5Cq)  ,  (IV-31a) 


where 


P5(q)  - 


(IV-31b) 


The  value  of  pQ  and  should  be  chosen  so  as  to  make  the  Impulse 
approximation  valid  (see  Eq.  (IV-13))  and  so  as  to  optimally  factorize 
the  plon  momentum  space  representation  of  the  lowest-order  optical 
potential.  The  final  expression  obtained  for  the  optical  potential  is 
simply  a  product  of  the  xN  T  matrix  and  the  nuclear  form  factor. 

The  xN  T  matrix  Involved  In  Eq.  (IV-31a)  is  not  the  free  xN  T 
matrix  In  the  xN  center-of-mass  system.  Denoting  <<' | £(fi>0) | <>  as  the 
free  xN  T  matrix  In  Che  xN  center-of-mass  system,  <<' | C(£>0) | <>  Is  defined 
only  for  the  on-energy-shell  values  | <' |  -  |<|  and  u>Q  -  “0(<).  where  x  Is 
the  plon-nucleon  relative  momentum.  However,  <ic'pQ-q  1 1  c(ui0)  |  j  kpQ> 
Involves  both  on-energy-shell  (|i?'|  -  j£|  -  |ic0|,  (<*»0(k0))  and  off- 

energy-shell  (|k'|  *  |ic|  *  |£Q|)  matrix  elements.  Furthermore,  this  T 
matrix  describes  xN  scattering  in  the  xA  center-of-mass  system.  One 
needs  to  relate  <£'pQ-q  1 1 1 ( u»Q)  |  j  itpQ>  to  <<' | C(«0) | *>,  since  <x' |  C(ffl0)  |  k> 
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is  the  T  matrix  given  by  the  experimentally  determined  xN  scattering 
phase  shifts.  The  transformation  is  ambiguous  and  one  method  is 
explained  in  detail  in  [Co-84],  Briefly,  one  employs  a  frozen-target 
approximation  (pQ  -  -kQ/A)  and  assumes  effective  on-energy-shell 
kinematics.  <£'p0-q| | t(«Q) j |£p0>  is  approximated  as 

<£'p0-q 1 1 1<  w0) | | £p0>  *  <!'  I  1 1< Wg)  |  |f> 

-  Y<K'|C(ffi0)|K>  ,  (IV-32a) 


where  k  is  an  effective  on-energy-shell  momentum  [Co-84]  associated  with 
the  xN  collision  energy  in  the  xA  center-of-mass  system  and 


C  *>EN< k  )  1 1/2 

M  )EM(je/ A) En(V7"A)  J 


(lV-32b) 


Ew ( k)  -  (x2  +  m2)1^2,  En(<)  -  (it2  +  m2)1/2,  etc.  is  given  by 

%.  "  M<)  +  en<JS/a)  •  ( IV— 33  > 

Since  one  is  assuming  effective  on-ener  ’,y-shell  kinematics,  all  kinematic 
variables  (Jc,  w^,  and  u>Q)  are  related  through  use  of  the 

invariance  of  the  four-vector  product  s  -  (P*|  +  pjJ)(Pnu  +  PNu).  The 
effective  on-energy-shell  kinematics  are  phenomenologically  determined  by 
Incorporating  an  energy  shift  parameter  in  EwOO  [Co-84],  i.e., 

Ex<£>  "  Mk>  +  A  •  CIV-34) 
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A  is  a  free  parameter  which  adjusts  the  kinematics  so  that  the  approxi¬ 
mations  (impulse,  factorization,  frozen-target)  made  in  deriving  the 
lowest-order  optical  potential  are  most  correct  [Co-84],  The  final 
result  for  the  pion  momentum  space  representation  of  the  lowest-order 
optical  potential  is 

<*'|Vop|fc>  -  AY<<'|e(uS0)|x>p(q)  .  (IV-35) 

The  connection  between  the  optical  potential  and  the  vN  phase 
shifts  is  made  in  the  usual  way  for  on-energy- shell  Kinematics  [Ei-80] — 

<K'|e<fi„)|£>  -  ,  (iv-36) 

O  *  P 

where  u  is  the  reduced  energy  for  the  pion  and  nucleon,  u  *  [£„(<)  1  + 
EN(«r‘]'‘,  and  F( tc'+x)  is  the  scattering  amplitude  matrix  discussed  in 
Chap.  1.  With  Eq.  (1-7)  and  neglecting  the  spin-dependent  term,  the 
optical  potential  is 

^'|Vop|k>  -  ^Ay[a(S0)  +  b(5>0)x',K]p(q)  ,  (IV-37) 

which  is  known  as  the  Kisslinger  optical  potential  [Ki-55],  Transforming 
to  the  coordinate  space  representation. 


VQp(r)  -  ~Ay[-a(u0)p(r)  +  b(ui0)7*p(r)V]  . 


(IV-38) 
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In  Eq.  (IV-38),  p(r)  is  the  ground  state  nuclear  density  normalized  to 
unity,  /p(r)d3r  -  1.  The  ground  state  nuclear  density  may  be  expressed 
in  terms  of  the  neutron  and  proton  densities  as 


p(r )  -  p«(r)  +  Pz(r)  ,  PN(r)  *  ^p(r)  »  Pz<f)  *  70(f)  *  (IV-39) 

““’A  A 


where  N  is  the  number  of  neutrons  in  the  nucleus  and  Z  the  number  of 
protons.  Substituting  Eq.  (IV-39)  into  Eq.  (IV-38),  the  Kissllnger 
optical  potential  becomes 


Vop(0  »  ^AY[-a(<SQ)PN(?)  -  a(ffi0)pz(r)  + 


b(ao)V«pN(r)V  +  b(S0)V*pz(r)v]  . 


(IV-40) 


A  similar  expression  holds  for  elastic  scattering  from  a  nucleus, 
which  is  obtained  from  Eq.  (IV-41)  by  replacing  Pjj(r)  with  pz(r)  and  vice 
versa.  For  self-conjugate  nuclei  (N  *  Z)  such  as  6Li,  the  Kissllnger 
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STATES  OF  LITHIUH-6CU)  AIR  FORCE  INST  OF  TECH 
WRIGHT-PATTERSON  AFB  OH  R  R  KIZIRH  DEC  84 
UNCLASSIFIED  AFITXCIXNR-84-85D  FXG  20X8  NL 
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optical  potentials  for  *+  and  w"  elastic  scattering  are  Identical  as  seen 
from  Eqs.  (IV-39)  and  (IV-41). 

Differential  cross  sections  for  r  A  elastic  scattering  are 
obtained  by  numerically  solving  a  Klein-Gordon  equation  containing  the 
Kisslinger  optical  potential  and  the  Coulomb  potential,  Vc.  The  radial 
piece  of  the  Klein-Gordon  equation  is 

+  ~J^K<r)  -  (j£  +  ^  -  25Uop  -  25Uc)ut(r)  ,  (IV-42) 

where  Uc  -  Vc(r)/fic,  UQp  •  vop(r)/fi c,  5  *  E/fic,  k  is  the  incident 
momentum  of  the  pion  in  the  wA  center-of-mass  system,  and  E  is  the  total 
energy  of  the  Incident  pion  in  the  *A  center-of-mass  system.  For  the 
elastic  calculations  presented  in  this  dissertation,  a  modified  version 
[Co-80]  of  the  computer  code  PIRK  [El— 74]  was  used  to  solve  Eq.  (IV-42) 
for  the  phase  shifts  needed  to  compute  the  elastic  differential  cross 
sections.  The  modifications  to  the  code  Include  the  transformation  of 
the  i»N  T  matrix  in  the  irA  center-of-mass  system  to  the  uN  center-of-mass 
system  according  to  Eq.  (IV-32)  and  linking  PIRK  co  the  MINUIT  optimizer 
package  [Ja-75a]  in  order  to  search  on  the  energy  shift  parameter 
incorporated  in  the  effective  on-energy-shell  kinematics,  Eq.  (IV-34) 
[Co-80].  The  MINUIT  optimizer  package  also  allows  searching  on  the 
shapes  and  magnitudes  of  the  neutron  and  proton  density  distributions 
which  enter  into  the  Kisslinger  optical  potential.  However,  for  the  6Li 
elastic  calculations,  the  proton  density  distribution  was  taken  directly 
from  electron  scattering  measurements  with  no  variation  in  the  shape  and 
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magnitude,  but  corrected  for  the  finite  size  of  the  proton  charge.  The 
sane  density  distribution  was  used  for  the  neutrons.  Lastly,  the  *N 
partial-wave  amplitudes  in  Eq.  (IV-41)  were  constructed  from  the  *N  phase 
shifts  of  Rowe,  Salomon,  and  Landau  [Ro-78]. 

C.  Inelastic  scattering 


The  theory  for  *A  inelastic  scattering  is  more  complicated  than 
for  irA  elastic  scattering.  As  for  elastic  scattering,  one  uses  the  non- 
relativistic  multiple-scattering  formalism.  However,  »A  Inelastic 
scattering  involves  a  "hard”  interaction  which  induces  the  nuclear 
transition,  an  interaction  not  present  in  irA  elastic  scattering  [Ei-80]. 
To  derive  the  transition  amplitude  for  irA  inelastic  scattering,  the  main 
assumption  is  that  the  reaction  mechanism  consists  of  two  parts:  (1) 
many  scattering  steps  which  serve  only  to  distort  the  pion  wave  function 
and  in  which  the  nucleus  remains  in  the  ground  state  and  (2)  a  one-step 
"hard"  interaction  which  causes  the  nuclear  excitation  [Ei-80], 
Derivation  of  the  transition  amplitude  involves  the  same  approximations 
as  the  derivation  of  the  lowest-order  elastic  optical  potential  in  Sec. 
IV-B,  although  the  approximations  may  be  less  well  founded. 

Beginning  with  Eq.  (IV-I4)  and  using  the  large  nucleus  approxi¬ 
mation. 


A 

T  -  T  +  TG„T,  T  -  E  t, 
0  1-1  1 


(IV-43) 


Since  one  is  interested  in  inelastic  scattering,  matrix  elements  of  Eq 
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(IV-43)  must  be  caken  between  a  nuclear  excited  state,  |0>,  and  the 
nuclear  ground  state,  |0> — 


< 0 1 T 1 0>  -  <S|T|0>  +  <8|TGoT|0>  . 


(IV-44) 


Separating  out  successively  higher  orders  of  diagonal  matrix  elements, 
Eq.  (IV-44)  Is 


<0|T|O>  -  < 0 1 T 1 0>  +  < 0 1 T 1 0>GQ ( 0 ) < 0 1 T ) 0> 
+  <0|T| 0>GO<  0><0|TjO>  + 


l  <0 | T | e>G0( e)<e|T| 0>  , 
e*0 ,  p 


(IV-45a) 


where 


Go<°>  -  E  -  k;  V~ln  ’  Go(S) 


E  -  Egucl  -  K,  +  in 


n-0+  .  (IV-45b) 


Just  as  in  the  elastic  scattering  case,  the  last  term  of  Eq.  (IV-45)  is 
ignored  since  it  contains  more  nondiagonal  nuclear  matrix  elements  than 
the  first  three  terms.  Therefore, 


< 0 1 T| 0>  »  < 0 J  T 1 0>  +  < 0 1 T j 0>GQ ( 0 ) < 0 1 T j 0> 
+  <0\T| 0>GoC  0)<0jT|O>  . 


(IV-46) 


-  - 

•  J 


Rewriting , 
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<8|T|0>  -  {l  -  <8|Tj  0X?o( 6)]“1<8|T|0>  * 

tl  +  Go(0)<0|T|0>}  .  (IV-47) 

In  order  Co  obtain  a  final  expression  for  Che  transition  amplitude,  one 

needs  to  calculate  the  matrix  elements  of  <0 1 Tf 0>  in  the  momentum  space 

of  the  pion.  Letting  j  $£  >  and  J  $£»>  denote  the  plon  initial  and  final 

o  o 

plane  wave  states,  (E  -  >  -  0, 

o 

<*jf-8|T|«  0>  *  <^J[1  -  <S|T|0>Go(0)]‘l<8|T|O>  * 

0  0  0 

[1  +  Go(0)<0|t|0>]|«  >  .  (IV-48) 

o 

|  ♦J  >  +  Go(0)<0]T|0>|^  >  is  the  equation  for  the  plon  scattering  state 
o  o 

in  which  the  complete  optical  potential,  <0|T|0>,  is  involved,  i.e., 

I  4  >  -  >  +  Go(0)<0|T|0>|«  >  ,  (IV-49a) 

0  0  "’O 


where 


(E  -  K,  -  <0|T|0»|4  >  -  0  .  (IV-49b) 

o 

Likewise,  <$j£*|[l  -  < 0 j T |  0>GO(  0) ] ~ 1  is  the  equation  for  the  pion 
scattering  state  in  which  the  lowest-order  optical  potential  is  involved, 
i.e. , 


<XM  -  <#M  +  <xjJ-|<3|T|8>Go(0)  , 
0  0  0 


(IV- 50a) 


■  /air 
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where 

(E  -  Kt  -  <8|T|8»jx|^>  -  0  .  ( IV-50b) 

Using  Eqs.  (IV-49)  and  (IV-50)  in  Eq.  (IV-48),  the  final  result  for  the 
transition  amplitude,  which  describes  *A  inelastic  scattering,  la 

0>  *  <xj'S|T|^  0>  .  (IV-51) 

o  o  o  o 

This  expression  is  referred  to  as  the  distorted-wave  Impulse  approxi- 

f 

matlon  (OWIA)  description  of  *A  inelastic  scattering.  In  practice, 

♦ 

is  replaced  by  xiT  and  I  8>  in  Eq.  (IV-50)  by  |0>,  so  that 
o 

<^-8|T|^o0>  *  <  X^- 1  <  e  |  T 1 0>  I  x^o>  •  ( IV— 52) 

Thus,  the  DWIA  expression  for  the  transition  amplitude  involves  the 
initial  and  final  plon  distorted  waves,  in  which  the  distorting  potential 
is  the  lowest-order  elastic  optical  potential,  <0|T|0>,  and  an  inelastic 
interaction,  <ftjT|0>,  which  induces  the  nuclear  transition. 

Inserting 

1  "  ffa *  <™-53) 

where  |k>  represents  the  plon  state  with  initial  momentum  k  in  the  irA 
center-of-mass  system,  Eq.  (IV-52)  becomes 


;mwJ«;<£')Ji<£'f'|t‘<“>|£o>x£V£)  •  (IV-54> 


+  ^  +  ♦ 

X?  (k)  -  <k|x£  >  and  w  is  an  appropriate  *N  collision  energy  in  the  *A 
center-of-mass  system.  The  spin  and  isospin  dependences  for  both  the 
nucleus  and  pion  are  not  shown  explicitly  in  Eq.  (IV-54).  Denoting 
Upo(k'.£)  as  <It'6|tl(w)|itO>,  Ugo(ii',ie)  la  analogous  to  the  pion 
momentum  space  representation  of  VQp  (see  Eq.  (IV-24)).  Therefore,  after 
the  same  manipulations  as  performed  for  <ic'|VQp|ic>, 


«»<£'.£>  ■  A|<k'p0-^||t5(ui0)||kp0>p5(q)  , 


(IV-55a) 


where 


p^(q)  -  pP^(P-q.P)  . 

4  4  >  .  d  P2  d  +,  A  ♦  ♦ 

p5(p-q.p)  -  /THT5-*'*nwy8’6(p  +  i-2pl  +  k°)05 


(IV-55b) 


*B(p-qP2***PA^o^PP2,"PA) 


(IV-55c) 


The  p^'s  are  the  nuclear  spin/isospin  scalar  and  vector  transition 
densities  and  measure  the  ability  of  a  nucleon  in  the  state  |p>  to  be 
scattered  into  the  state  |p-q>,  while  the  nucleus  is  excited  from  the 
ground  state,  |0>,  to  the  excited  state,  |f3>,  during  a  particular 
spin/lsospin  transition  (La-73, Le-74] . 
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As  for  che  elastic  scattering  case,  the  *N  T  matrix  involved  in 
the  inelastic  interaction  involves  both  on-energy-shell  and  off-energy- 
shell  matrix  elements  and  is  not  the  free  irN  T  matrix  in  the  *N  center- 
of-mass  system.  The  transformation  is  explained  in  detail  in  [Co-84] 
and,  as  for  the  elastic  case,  employs  a  frozen-target  approximation  and 
effective  on-energy-shell  kinematics.  Following  the  development  of  Eqs. 
(IV-32)  through  (IV-33),  the  effective  on-energy-shell  kinematics  are 
phenomenologically  determined  by 

E*(<)  *  Mk>  ~Ex/2  +  A  •  (IV-56) 

where  Ex  is  the  excitation  energy  for  the  nuclear  state  |0>,  and  A  is 
taken  from  elastic  scattering  [Co-84).  Also,  the  collision  energies  for 
the  initial  and  final  pion  distorted  waves  are  E*(£)  *  Eff(k)  +  A  and 
E$(jO  ■  E„(k)  -  Ex  +  A,  respectively. 

Following  the  development  for  the  elastic  scattering  case  but 
including  both  the  spin-dependent  and  spin-independent  terms  of  Eq. 
(1-7),  the  coordinate  space  representation  of  the  DWIA  expression  for  the 
transition  amplitude  is 

^80  *  Jd^rX~*(E$(jc* )  ,r )UgQ( r )x+(E$(_k)  ,r )  ,  (IV-57a) 


with 


U0O(^  “  •7A^t-aSS^o>pss(^)  "  •SV(&0)psv(*)  + 


b8S(0o)V«psg(?)V  +  bsv(ffl0)7.psv(?)V  + 

icvs(ffl0)Vxpvs(?)7.ft  +  ic^fiJVxp^rJV.ftl  . 


o'  ^w' 


( IV-57b) 


The  nuclear  translclon  densities  are  given  by 


p£<r)  *  J1£1/d3rl,,*d3rA*0(^r**fA>5(r-r1)O51*o(?1“*?A) 

•if  <B|  6(r-r .  )0- J0>  . 

Ai-l  1 


(IV-58) 


The  nuclear  transicion  densities  contain  the  nuclear  structure  input  to 
irA  inelastic  calculations,  and  when  possible,  are  taken  from  other 
projectile-nucleus  scattering  measurements  such  as  electron  scattering. 

As  mentioned  in  Chap.  I,  the  above  DWIA  formalism  is  Che  most 
prevalent  theoretical  description  of  *A  inelastic  scattering,  and  cal¬ 
culations  using  the  DWIA  and  well-known  nuclear  transition  densities  are 
adequate  for  describing  the  pion-induced  excitation  of  nuclear  states 
which  are  strongly  excited.  A  few  examples  of  the  success  of  the  DWIA 
were  given  in  Chap.  I.  In  many  of  the  successful  applications  of  the 
DWIA  to  the  description  of  xA  inelastic  scattering,  the  nuclear  structure 
input  is  fixed  from  electron-nucleus  inelastic  scattering  (e,e')  data. 
Since  the  essential  parts  of  the  Inelastic  interaction  are  the  xN  T 
matrix  and  the  nuclear  transition  densities  (see  Eq.  (IV-55)),  un¬ 
certainties  in  the  inelastic  Interaction  are  a  result  of  uncertainties  in 
the  xN  T  matrix.  In  other  words,  within  the  framework  of  the  DWIA,  if 
the  transition  densities  determined  from  fits  to  (e,e'>  data  do  not  give 
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agreement  between  theory  and  data  for  xa  inelastic  scattering  data,  the 
aost  reasonable  part  of  the  DWIA  calculation  to  question  is  the  xN 
interaction,  especially  if  the  expressions  for  (e,e')  and  xA  inelastic 
scattering  Involve  the  transition  densities  in  the  same  manner.  The 
examples  of  Chap.  I,  which  were  selective  of  the  spin-dependent  piece  of 
the  xN  interaction,  suggest  the  correctness  of  the  parameterization  of 
the  spin-dependent  piece  of  the  xN  interaction  [Ca-83].  Furthermore, 
Morris,  et  al.  [Mo-81]  and  Boyer,  et  al.  [Bo-81]  obtain  good  agreement 
between  angular  distributions  and  DWIA  calculations,  which  use  transition 
densities  determined  from  (e,e')  data,  for  inelastic  scattering  to  low- 
lying  collective  states  in  12C,  **°Ca,  “*2Ca,  4‘*Ca,  and  l*8Ca.  The  good 
agreement  suggests  the  correctness  of  the  parameterization  of  the  spin- 
independent  piece  of  the  xN  Interaction.  Since  the  DWIA  calculations  for 
pion  Inelastic  scattering  to  the  2.185-  and  3.563-MeV  states  of  6Li 
presented  in  Chap.  V  use  transition  densities  determined  from  fits  to 
(e,e')  data,  the  following  paragraphs,  taken  from  [Pe-79,Pe-81] ,  briefly 
outline  the  relationship  between  electron-nucleus  and  xA  inelastic 
scattering  differential  cross  sections. 

In  order  to  obtain  clear  and  simple  relations,  the  plane-wave 
Born  approximation  and  plane-wave  impulse  approximation  (PWIA)  are  used 
in  deriving  the  differential  cross  sections  for  electron-nucleus  and  xA 
inelastic  scattering,  respectively.  For  xA  inelastic  scattering,  the 
PWIA  expression  for  the  transition  amplitude  is  given  by  Eq.  (IV-57)  with 
the  initial  and  final  pion  distorted  waves  replaced  by  plane  waves.  The 
expressions  for  the  differential  cross  sections  will  be  given  for 
transitions  in  which  only  one  value  of  J,  the  total  angular  momentum 
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transfer.  Is  involved.  (Although  the  ground  state  of  6Li  has  spin  and 
parity  J*  •  1+,  there  is  only  one  dominant  J  transfer  for  the  transitions 
to  the  3+,  2.185-MeV  and  2+,  4.25-MeV  states.)  For  natural-parity 
transitions  (A*  -  (-1)J,  where  A*  is  the  parity  change),  the  differential 
cross  sections  for  electron-nucleus  and  irA  inelastic  scattering  are 


lj{tS<q)Pa(<l>  +  eircaS(q)pJo(q>>l2 


+  ^f(9)|^S(q)p^(q)|2 


( IV-59a) 


( IV-59b) 


In  Eq.  (IV-59a),  |FL(q)|2  Is  the  longitudinal  form  factor  and  Che 

equation  is  valid  only  for  transitions  which  are  dominantly  longitudinal. 
Pp( q )  is  the  Fourier-Bessel  transform  of  the  proton  transition  density. 
This  transition  density  includes  the  finite  size  of  the  proton  charge  in 
the  electron-nucleus  result  but  not  in  the  vA  result.  In  Eq.  (IV-59b), 
is  a  scaling  factor,  f(8)  is  an  angle  factor,  a  is  an  isospin  index, 
tc(q)  and  t^S(q)  are  the  Fourier-Bessel  transforms  of  the  spin- 
independent  and  spin-dependent  pieces  of  the  tin  interaction,  and  p2s(q) 
and  Pjj(q)  are  Fourier-Bessel  transforms  of  the  spin-orbit  and  spin 
transition  densities.  If  p®j  *  PjS  *  0  and  pp  -  pjj  as  in  N-Z  nuclei,  the 
differential  cross  sections  for  electron-nucleus  and  wA  inelastic 
scattering  involve  only  a  single  density  and  are  directly  related.  This 


Is  Che  basis  for  fixing  the  transition  densities  needed  in  tA  inelastic 
scattering  calculations  from  (e,e')  measurements. 


For  unnatural-parity  transitions  (At  *  (-1)J+*),  the  differential 
cross  sections  for  electron-nucleus  and  tA  inelastic  scattering  are 


do* 

*3IT 


(IV-60a) 


^«^f(0)|EtJ;s(q)p»a(q)| 


LS /  _  \  _ 8Ci/ 


(IV-60b) 


I^T^^I2  19  the  transverse  magnetic  form  factor  and  the  transition 
densities  for  the  electron-nucleus  expression  include  the  finite  size  of 
the  proton  charge.  gg  and  are  the  spin  and  orbital  g-factors,  Pjj(q) 
is  the  Fourier-Bessel  transform  of  the  orbital  current  transition 
density,  and  pj(q)  is  the  transverse  linear  combination  of  spin 
transition  densities. 


■>J<q>  •  -<2ST>l/2“Wq)  -  • 


2J+i  ‘ 


(IV-61) 


Since  Eq.  (IV-60a)  involves  both  the  orbital  current  and  transverse  spin 
transition  densities,  there  is  not  a  direct  connection  between  the 
electron-nucleus  and  vA  inelastic  scattering  differential  cross  sections 
as  for  the  natural-parity  transitions.  Use  of  (e,e')  measurements 
requires  the  separation  of  the  orbital  current  and  transverse  spin 
transition  densities  from  the  transverse  magnetic  form  factor  and  depends 
upon  the  ground  state  and  excited  state  wave  functions.  However,  as  for 
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the  above  natural-parity  case,  the  theoretical  expression  for  pion- 
induced  unnatural-parity  transitions  is  particularly  simple,  involving 
only  a  single  piece  of  the  xN  interaction  and  a  single  transition 
density. 


V.  COMPARISON  OF  THEORY  AND  DATA 


A.  Elastic  scattering 

As  explained  In  Sec.  IV-B,  the  first-order,  zero-range,  Impulse- 
approximation  elastic  calculations  were  performed  with  a  modified  version 
[Co-80]  of  the  coordinate-space  computer  program  PIRK  [El-74].  The  same 
density  distribution  was  used  for  both  the  point  proton  and  point  neutron 
density  distributions  which  enter  Into  the  Klssllnger  optical  potential 
(see  Eq.  (IV-41)).  This  density  distribution  was  taken  from  elastic 
electron  scattering  measurements  [Li-71]  and  is  a  three-parameter 
phenomenological  distribution  of  the  form 


P(r) 


8* 


572[7*”p(i£ 


c2(6b2-r2) 

-exP(7r7JJ 


4b' 


4b*1 


(V-l) 


with  a  -  0.928  fm,  b  -  1.26  fm,  and  c  »  0.48  fm  [Li-71],  The  elastic 
electron  scattering  parameters  in  Eq.  (V-l)  were  corrected  for  Che  finite 
size  of  the  proton  charge  according  to 

<r2>pp  -  <r2>ch  -  (0.8)2  ,  (V-2) 

where  <r2>pp  is  the  mean-square  radius  of  the  point  proton  density, 
<r2>c^  is  the  mean-square  radius  of  the  charge  density  as  determined  from 
electron  scattering  measurements  "  2.56  -  0.05  fm)  [Li-71],  and 


0.8  fm  is  the  root-mean-square  radius  of  the  charge  distribution  for  a 


single  proton  [El-61]. 

The  measured  differential  cross  sections  for  it+  elastic 
scattering  from  ®Li  at  T,  -  120  and  180  MeV  and  the  differential  cross 
sections  for  »“  elastic  scattering  from  6Li  at  -  164  MeV  taken  from 
Zichy's  work  at  SIN  [Zi-80],  along  with  the  elastic  calculations,  are 
presented  in  Fig.  V-l.  Fig.  V-2  shows  the  elastic  calculations  and  *+ 
elastic  scattering  differential  cross  sections  for  100  <  T,  <  260  MeV  and 
q  *  109  MeV/c.  For  the  elastic  calculations,  the  partial-wave  amplitudes 
in  the  Kisslinger  optical  potential  (see  Eq.  (IV-41))  were  evaluated  at 
an  energy  of  30.0  MeV  below  the  Incident  pion  beam  energy,  i.e.,  in  Eq. 
(IV-34),  A  ■  -30.0  MeV.  This  value  of  A  was  chosen  from  the  work  of 
Cottlngame  and  Holtkamp  [Co-80].  They  found  that  a  value  of  A  »  -30  MeV 
gave  better  agreement  between  pion  elastic  scattering  data  and  cal¬ 
culations  for  nuclei  ranging  from  9Be  to  208Pb.  Elastic  calculations 
were  also  performed  for  the  120-  and  180-MeV  and  n-  164-MeV  data  sets 
with  a  one-parameter  search  on  the  shift  in  the  collision  energy.  The 
result  of  the  search  was  A  =  -20  MeV.  However,  the  elastic  calculations 
with  A  =■  -20  MeV  differ  only  slightly  from  those  with  A  =  -30  MeV. 

As  seen  from  Figs.  V-l  and  V-2,  there  is  good  agreement  between 
the  elastic  calculations  and  the  measured  elastic  differential  cross 
sections.  Since  the  lowest-order  elastic  optical  potential  is  the 
distorting  potential  used  for  the  generation  of  the  initial  and  final 
pion  distorted  waves  for  inelastic  calculations  (see  Eq.  (IV-52)),  the 
good  agreement  indicates  adequate  handling  of  the  distortions. 


dq/dn  (rnb/si) 
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Fig.  V-l:  Differential  cross  sections  for  and  it-  elastic  scattering 
for  ^Li  for  Tfl  »  120,  164,  and  180  MeV.  The  calculations  include  a  -30 
MeV  shift  in  the  energy  at  which  the  optical  model  parameters  are 
calculated.  The  164-MeV  data  are  from  [ Zi— 80 J . 


dc/df)  (mb/sr) 


B.  Inelastic  scattering 


The  Inelastic  calculations  presented  In  the  following  sections 
use  the  DWIA  expression  for  the  transition  amplitude  (see  Eq.  (IV-57)). 
The  initial  and  final  pion  distorted  waves  are  calculated  from  the 
lowest-order  elastic  optical  potential  employed  in  the  elastic  calcu¬ 
lations  using  the  coordinate-space  computer  program  UTDWPI  [Bo-n.d.]. 
Calculation  of  the  Inelastic  interaction,  Ugg,  uses  the  frozen-target 
approximation  and  effective  on-energy-shell  kinematics  and  a  collision 
energy  obtained  from  the  elastic  calculations  (see  Eq.  (IV-56)).  For  the 
inelastic  calculations  using  Cohen-Kurath  Intermediate  coupling  p-shell 
wave  functions  [Co-65]  and  pure  LS-coupling  p-shell  wave  functions,  a 
modification  of  the  generalized  inelastic  scattering  potential  code 
ALLWRLD  [Ca-84]  is  used  to  generate  the  nuclear  transition  densities  from 
harmonic  oscillator  wave  functions  and  subsequently  calculate  the 
inelastic  interaction.  The  harmonic  oscillator  parameters  and  renormali¬ 
zation  constants  (shape  and  strength  parameters  for  the  transition 
densities)  include  the  standard  center-of-mass  correction  needed  when 
using  shell-model  wave  functions.  The  differential  cross  sections  are 
then  generated  from  the  code  UTDWPI.  For  all  other  inelastic  calcu¬ 
lations,  only  the  code  UTDWPI  is  used  and  no  center-of-mass  correction  is 
included  in  the  harmonic  oscillator  parameters  and  renormalization 
constants.  However,  these  inelastic  calculations  use  center-of-mass 
corrected  transition  densities  with  the  correction  being  applied  in 
momentum  space  in  the  usual  way. 
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1.  3+,  T-0,  2.185-HeV  state 

Electron-nucleus  Inelastic  scattering  measurements  Indicate  that 
the  natural-parity  transition  to  the  2.183-MeV  level  of  6Li  Is  almost 
completely  longitudinal,  with  measurements  In  the  region  of  q  -  0.7  to 
1.8  fm”1  yielding  a  transverse  form  factor  which  Is  less  than  22  of  the 
longitudinal  form  factor  [Ne-71).  Since  the  spln-orblt  and  spin 
transition  densities  are  approximately  zero  [Ca-84a],  the  wA  inelastic 
Interaction,  U^q,  involves  only  the  spin-independent  central  component  of 
the  itN  Interaction  and  the  isoscalar  matter  transition  density  (see  Eq. 
(IV-59b)  and  discussion).  Inelastic  calculations  for  this  natural-parity 
transition  were  performed  using  two  different  transition  densities 
derived  from  pure  LS-coupling  and  Cohen-Kurath  Intermediate  coupling  p- 
shell  wave  functions.  Each  calculation  employed  equal  transition 
densities  for  the  protons  and  neutrons.  The  harmonic  oscillator 
parameter,  a  (a  *  (Mw/K)  ,  where  M  is  the  nucleon  mass  and  lieu  is  Che 
energy  quantum  of  the  harmonic  oscillator),  and  renormalization  constant 
for  the  transition  density  derived  from  the  pure  LS-coupling  p-shell  wave 
functions  are  0.534  fm"1  and  2.03.  These  values  were  taken  from  [Pe-82] 
and  were  deduced  in  [Br-72]  from  a  fit  to  the  (e,e')  data  of 
[Be-63,Ne-69] .  Using  these  values,  Petrovich,  et  al.  obtained  good 
agreement  between  theory  and  experiment  for  the  6Li(p,p')6Li*(2.185  MeV) 
differential  cross  sections  at  Ep  *  25  and  45  MeV  [Pe-82].  The 
transition  density  derived  from  the  intermediate  coupling  p-shell  wave 
functions  uses  a  *  0.558  fm_l  and  a  renormalization  constant  of  i.93, 
again  determined  from  (e,e')  data  [Ca-84a]. 


Electron  scattering  longitudinal  form  factors  (Be-76,Be-79)  are 
shown  In  Fig.  V-3.  The  solid  (dashed)  curves  correspond  to  the 
transition  densities  computed  from  the  pure  LS-  (intermediate)  coupling 
p-shell  wave  functions.  Both  theoretical  form  factors  are  similar,  with 
the  Intermediate  coupling  form  factor  in  better  agreement  with  the  (e,e') 
measurements.  The  differential  cross  sections  for  x+  inelastic 
scattering  to  the  2.185-MeV  state  and  OWIA  calculations  for  Tv  »  120  and 
180  MeV  are  presented  in  Fig.  V-4.  Fig.  V-3  shows  the  data  and  calcu¬ 
lations  for  Tv  *  100  to  260  MeV  with  the  differential  cross  sections 
corresponding  to  a  constant  momentum  transfer  q  *  109  MeV/c.  From  the 
data  and  calculations,  the  first  maxima  of  the  angular  distributions  for 
the  2.185-MeV  state  are  expected  to  be  at  q  *  164  MeV/c.  The 
Intermediate  coupling  calculation  is  in  good  agreement  with  both  the  120- 
and  180-MeV  experimental  angular  distributions.  Furthermore,  this  calcu¬ 
lation  reproduces  well  the  constant-q  experimental  differential  cross 
sections,  considering  that  the  theoretical  values  plotted  in  Fig.  V-5  are 
taken  from  the  steep  forward  slope  of  the  various  angular  distributions 
where  errors  would  produce  the  greatest  variations.  The  pure  LS-coupling 
calculation  yields  similar  shapes  for  the  120-  and  180-MeV  angular 
distributions  and  the  constant-q  differential  cross  sections  as  does  the 
intermediate  coupling  calculation  but  overestimates  the  magnitudes.  Such 
disagreement  suggests  that  the  renormalization  constant  is  too  large. 
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Fig.  V-4:  Angular  distributions  for  inelastic  scattering  to  the  3  , 
T-0,  2.185-MeV  state  of  6Li  for  Tw  -  120  and  180  MeV.  The  calculations 
used  isoscalar  matter  transition  densities  derived  from  pure  LS-coupling 
(solid  curve)  and  Cohen-Kurath  intermediate  coupling  (dashed  curve)  p- 
shell  wave  functions  with  a  »  0.534  fm"1  and  a  renormalization  constant 
of  2.03  and  a  -  0.558  fm~l  and  a  renormalization  constant  of  1.93, 
respectively. 


,*•»>>> 


.*  >  V 


6Ii(TT+,7r+y)8Li* 

3+  (ai85  MeV) 


i  i  i  i  i 


120  150  180  210  240  270 
T„  (MeV) 


Fig.  V-5:  Differential  cross  sections  for  inelastic  scattering  to  the 
3+,  T»0,  2.185-MeV  state  of  6Li  at  a  constant  q  *  109  MeV/c.  The 
calculations  used  isoscalar  matter  transition  densities  derived  from  pure 
LS-coupling  (solid  curve)  and  Cohen-Kurath  intermediate  coupling  (dashed 
curve)  p-shell  wave  functions  with  a  -  0.534  fm~l  and  a  renormalization 
constant  of  2.03  and  o  ■  0.558  fm”1  and  a  renormalization  constant  of 
1.93,  respectively. 


The  electron  scattering  fora  factors  for  both  natural-parity 
transitions  to  the  2.185-  and  4.25-MeV  levels  have  the  same  dependences 
on  the  momentum  transfer  [Be-63].  Therefore,  the  transition  to  the 
4.25-MeV  state  is  principally  longitudinal,  and  this  transition  was 
treated  as  completely  analogous  to  the  excitation  of  the  2.185-MeV  level. 
However,  even  though  the  4.25-MeV  state  is  observed  in  (e,e")  spectra,  a 
detailed  form  factor  has  not  been  measured  because  of  this  state's  large 
natural  width,  the  overlapping  5.37-MeV  state,  and  the  large  continuum 
background  [Be-82].  As  a  result,  a  transition  density  cannot  be  deduced 
from  (e,e')  data  as  was  the  case  for  the  transition  to  the  2.185-MeV 
level.  Thus,  a  transition  density  was  derived  from  Cohen-Kurath 
intermediate  coupling  p-shell  wave  functions  with  a  -  0.52  fm-1  and  a 
renormalization  constant  of  0.82  determined  from  fitting  the  (*+,w+') 
data  at  T„  ■  120  and  180  MeV.  This  transition  density  yields  a  radiative 
width,  F^O^  ’  of  3.02  eV  in  agreement  with  the  experimental  value  of 
5.4  ±  2.8  eV  [Ei-69].  The  4.25-MeV  experimental  and  theoretical  120-  and 
180-MeV  angular  distributions  and  the  constant-q  differential  cross 
sections  are  shown  in  Figs.  V-6  and  V-7,  respectively.  The  agreement  is 
not  as  good  as  for  the  2.185-MeV  state  but  is  very  reasonable  considering 
the  difficulties  in  extracting  the  cross  sections  and  possible  un¬ 
certainties  in  the  transition  density. 
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Fig.  V-6:  Angular  distributions  for  *'  inelastic  scattering  to  the  2T, 
T«0,  4.25-MeV  state  of  6Li  for  *  120  and  180  MeV.  The  calculation 
used  an  lsoscalar  matter  transition  density  derived  from  Cohen-Kurath 
Intermediate  coupling  p-shell  wave  functions  with  a  «  0.52  fm“l  and  a 
renormalization  constant  of  0.82. 
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Fig.  V-7 :  Differential  cross  sections  for  w+  inelastic  scattering  to  the 
2+,  T-0,  4.25-MeV  state  of  6L1  at  q  *  109  MeV/c.  The  calculation  used  an 
isoscalar  matter  transition  density  derived  from  Cohen-Kurath 
intermediate  coupling  p-shell  wave  functions  with  a  *  0.52  fm”1  and  a 
renormalization  constant  of  0.82. 
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3.  0+,  T-l,  3.563-MeV  state 

As  is  the  case  for  the  natural-parity  transitions  to  the  2.185- 
and  4.25-MeV  levels  of  8Li,  the  DWIA  description  of  Che  pion-induced 
unnatural-parity  transition  to  the  3.563-MeV  state  is  straightforward, 
involving  only  a  single  component  of  the  irN  interaction  and  a  single 
nuclear  transition  density  (see  Eq.  (IV-60b)).  The  spin-dependent  piece 
of  the  i»N  interaction  is  represented  by  the  zero-range  spin-orbit 
operator  [Ca-83]  which  Carr,  et  al.  used  for  their  successful  analysis 
of  the  stretched  4”  and  6-  states  in  l60  and  28Si,  and  the  transition 
density  is  the  transverse  spin  transition  density.  For  the  following 
DWIA  calculations,  three  different  transverse  spin  transition  densities 
were  used,  each  transition  density  determined  from  fits  to  various  (e,e') 
measurements  but  compared  to  the  most  recent  sets  of  (e,e')  data 
[Be-75,Be-79l .  The  proton  and  neutron  transition  densities  were  taken  to 
be  equal.  The  reasons  for  using  three  transition  densities  are:  (1)  The 
p-shell  harmonic  oscillator  basis  is  known  to  result  in  a  poor  descrip¬ 
tion  of  the  3.563-MeV  (e,e')  transverse  magnetic  form  factor  over  the 
entire  second  lobe,  q  >  1.4  fm-1  ( Be— 7 5 ] .  Therefore,  one  transition 
density  which  reproduces  both  the  first  and  second  lobes  of  the  form 
factor  is  used.  (2)  Since  there  is  no  antianalog  to  the  3.563-MeV  state 
as  is  the  situation  for  the  12.71-/15.11-MeV  antianalog-analog  pair  of 
12C,  any  anomalies  in  the  3.563-MeV  excitation  function  depend  upon 
comparison  of  theory  with  data  and  not  data  for  analogs  as  for  the 
15.11-MeV  level.  Thus,  accurate  spin  transition  densities  are 


imperative 


The  first  analysis  for  the  spin-flip  transition  to  the  3.563-MeV 
state  used  a  transverse  spin  transition  density  derived  from  Cohen-Kurath 
Intermediate  coupling  p-shell  wave  functions.  The  harmonic  oscillator 
parameter  was  chosen  to  be  a  -  0.518  fm-1  from  the  work  of  Petrovich, 
et  al.  [Pe-82],  who  fitted  the  (e,e')  data  of  [Ne-71 ,Be-75,Be-79]  using 
pure  LS-coupling  p-shell  wave  functions.  In  the  fitting  procedure,  the 
lower  q  data  were  given  the  most  weight.  The  renormalization  constant 
was  determined  to  be  0.97  from  a  fit  to  the  most  recent  sets  of  (e,e') 
data  for  the  form  factor.  Petrovich,  et  al.  conclude  that  both  sets  of 
wave  functions  provide  an  adequate  description  of  the  experimental  static 
moments  and  transition  probabilities  of  6Li  with  the  exception  of  the 
quadrupole  moment  [Pe-82].  However,  the  Cohen-Kurath  intermediate 
coupling  p-shell  wave  functions  were  chosen  because  they  give  a  slightly 
better  fit  to  the  most  recent  sets  of  (e,e')  data  [Be-75,Be-79]  for  the 
inelastic  Ml  form  factor.  Fig.  V-8  shows  the  calculated  transverse 
magnetic  form  factor  (solid  curve)  and  the  (e,e')  data  for  the  transition 
to  the  3.563-MeV  state. 

The  second  DWIA  calculation  employed  a  transverse  spin  transition 
density  obtained  from  the  empirical  shell-model  wave  functions  of 
Donnelly  and  Walecka  [Do-73].  These  wave  functions  are  almost  identical 
to  the  pure  LS-coupling  p-shell  wave  functions.  Using  p-shell  harmonic 
oscillator  radial  wave  functions  for  the  valence  nucleons,  Donnelly  and 
Walecka  determined  Che  one  body  density  matrix  elements  from  normali¬ 
zation  conditions  and  from  fits  to  the  ground  state  magnetic  dipole  and 
electric  quadrupole  moments  and  the  Ml  form  factors  for  electron  elastic 
and  inelastic  scattering  for  q  <  1.01  fm-1  [Do-73].  A  harmonic 
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Fig.  V-8:  |FT(q)|2  for  the  0+,  T-l,  3.563-MeV  state  of  6Li.  The 
calculations  used  transition  densities  derived  from  Cohen-Kurath 
intermediate  coupling  p-shell  wave  functions  (solid  curve)  with  a  -  0.518 
fm-1  and  a  renormalization  constant  of  0.97,  transition  densities  derived 
from  the  empirical  shell-model  wave  functions  of  Donnelly  and  Walecka 
(dashed  curve)  with  a  -  0.493  fm”1  and  a  renormalization  constant  of 
0.96,  and  phenomenological  transition  densities  based  on  the  work  of 
Bergstrom,  et  al.  (chain-dot  curve).  The  data  are  from  [Be-75]  (open 
circles)  and  [ Be— 79 ]  (solid  circles). 


oscillator  parameter  [Do-73]  a  *  0.493  fas'-'  and  a  renormalization 
constant  of  0.96  (determined  from  a  fit  to  the  most  recent  sets  of  (e,e') 
data)  were  used  in  this  second  set  of  calculations.  The  resultant 
inelastic  Ml  form  factor  (dashed  curve)  is  compared  to  the  (e,e')  data  in 
Fig.  V-8. 

The  third  analysis  used  a  phenomenological  transverse  spin 
transiton  density  based  on  the  work  of  Bergstrom,  et  al.  [Be-79].  This 
transition  density  was  derived  in  the  same  manner  as  the  transition 
density  of  Donnelly  and  Walecka  except  that  the  (e,e")  data  of 
[Be-75, Be-79] ,  which  extend  to  q  »  2.96  fm-1 ,  were  fitted  with  a 
polynomial  form  for  the  p-shell  radial  transiton  density.  The  wave 
functions  for  the  6Li  ground  state  and  3.563-MeV  state  were  taken  to  be 
described  by  the  SASK-A  amplitudes  of  Bergstrom,  et  al.  [Be-79],  The 
radial  transition  density  has  a  phenomenological  form 

R2(r)  -  exp(-r2/b2)(a2r2  +  a^r4  +  a6r6)  ,  (V-3) 

with  b  -  2.02  fm,  a2  *  6.625  *  10“2  fm-5,  *  -5.036  x  10“3  fm-7 , 
a^  -  1.967  x  10-4  fm-9,  yielding  a  reduced  x2  of  1.32  from  a  fit  to  the 
inelastic  Ml  form  factor,  which  is  presented  in  Fig.  V-8  (chain-dot 
curve) . 

The  it*  inelastic  scattering  measurements  and  DWIA  calculations  at 
Tff  -  120  and  180  MeV  for  the  spin-flip  transition  to  the  3.563-MeV  state 
are  shown  in  Fig.  V-9.  At  120  MeV  the  three  transverse  spin  transition 
densities  give  equivalent  shapes  for  the  angular  distribution  in  the 
range  of  the  (n+,it+')  data,  15.7°  <  9r  _  <  48.9°. 
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Fig.  V-9:  Angular  distributions  for  ir  inelastic  scattering  to  the  0+, 
T-l,  3.563-MeV  state  of  6Li  for  T*  -  120  and  180  MeV.  The  calculations 
used  a  transverse  spin  transition  density  derived  from  Cohen-Kurath 
intermediate  coupling  p-shell  wave  funcitons  (solid  curve)  with  a  -  0.518 
fm  and  a  renormalization  constant  of  0.97,  a  transverse  spin  transition 
density  derived  from  the  empirical  shell-model  wave  functions  of  Donnelly 
and  Walecka  (dashed  curve)  with  a  -  0.493  fm”1  and  a  renormalization 
constant  of  0.96,  and  a  phenomenological  transverse  spin  transition 
density  based  on  the  work  of  Bergstrom,  et  al.  (chain-dot  curve). 
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Fig.  V-10:  Excitation  function  at  a  constant  q  *  109  MeV/c  for  x+ 
Inelastic  scattering  to  the  0+,  T-l,  3.563-MeV  state  of  6Li.  The  calcu¬ 
lations  used  a  transverse  spin  transition  density  derived  from  Cohen- 
Kurath  Intermediate  coupling  p-shell  wave  funcitons  (solid  curve)  with 
am  0.518  fm~*  and  a  renormalization  constant  of  0.97,  a  transverse  spin 
transition  density  derived  from  the  empirical  shell-model  wave  functions 
of  Donnelly  and  Walecka  (dashed  curve)  with  a  -  0.493  fm-1  and  a 
renormalization  constant  of  0.96,  and  a  phenomenological  transverse  spin 
transition  density  based  on  the  work  of  Bergstrom,  etal.  (chain-dot 
curve) . 


predict  the  correct  locatioa  for  the  first  maximum  of  0c.m.  “  28°* 
However,  the  Cohen-Kurath  intermediate  coupling  (solid  curve),  Donnelly 
and  Walecka  (dashed  curve),  and  phenomenological  (chain-dot  curve)  calcu¬ 
lations  overestimate  the  magnitude  of  the  first  maximum  by  18%,  29%,  and 
38%,  respectively.  The  three  transverse  spin  transition  densities  yield 
similar  results  for  the  shape  of  the  180-MeV  angular  distribution  through 
the  first  minimum  but  differ  through  the  second  maximum  and  minimum. 
Each  calculation  underestimates  the  magnitude  of  the  first  maximum  by  at 
least  26%.  Fig.  V-10  shows  the  *+  3.563-MeV  excitation  function  data  at 
q  *  109  MeV/c  and  DW1A  analyses.  None  of  the  transverse  spin  transition 
densities  predict  the  measured  shape  and  magnitude  of  the  excitation 
function.  The  three  DWIA  calculations  disagree  with  the  data  by  at  least 
a  20%  at  both  the  low  and  high  incident  pion  energies. 

C.  Discussion  and  results 

The  disagreement  between  theory  and  data  for  the  120-  and  180-MeV 
angular  distributions  and  excitation  function  for  the  unnatural-parity 
transition  to  the  3.563-MeV  state  is  difficult  to  understand,  as  is  the 
disagreement  observed  for  the  unnatural-parity  transition  to  the 
15.11-MeV  state  of  12C.  For  comparison,  the  12.71-  and  15.11-MeV 
excitation  functions  along  with  DWIA  calculations  as  described  in  [Mo-821 
are  shown  in  Fig.  V— 11.  As  noted  in  Chap.  I,  a  simple  DWIA  description 
of  vA  inelastic  scattering,  which  uses  a  single  piece  of  the  *N 
interaction  and  a  single  nuclear  transition  density,  has  been  successful 
in  describing  unnatural-parity  transitions  in  many  nuclei.  Using  the 


zero-range  spin-orbit  operator  and  a  transverse  spin  transition  density 
derived  from  Cohen-Kurath  p-shell  wave  functions,  Morris,  et  al.  [Mo-82 J 
adequately  reproduced  the  measured  12.71-MeV  excitation  function. 
Furthermore,  Cottlngame,  et  al.  [Co-84]  adequately  describe  the 
12.71-MeV  angular  distributions  for  -  100  to  260  MeV.  Even  though 
there  is  an  energy-dependent  enhancement  near  -  180  MeV  in  the 
15.11-MeV  excitation  function,  the  low-energy  (T„  •  100  and  116  MeV) 
angular  distributions  for  this  state  are  reproduced  by  simple  DWIA  calcu¬ 
lations  [Co-84].  Using  the  same  zero-range  spin-orbit  operator  for  the 
representation  of  the  spin-dependent  piece  of  the  wN  interaction  and 
three  different  transverse  spin  transition  densities,  DWIA  calculations 
do  not  adequately  reproduce  either  the  120-MeV  and  180-MeV  angular 
distributions  or  the  excitation  function  for  the  transition  to  the 
3.563-MeV  state.  However,  the  measured  angular  distributions  and 
constant-q  differential  cross  sections  at  energies  from  *  100  to  260 
MeV  for  Che  natural-parity  transitions  to  the  2.185-  and  4.25-MeV  levels 
of  6Li  are  reproduced  using  only  the  spin-independent  central  piece  of 
the  *N  interaction  and  transition  densities  derived  from  Cohen-Kurath 
intermediate  coupling  p-shell  wave  functions. 

The  configuration-space,  zero-range  form  of  the  spin-dependent 
component  of  the  «N  interaction  of  Carr,  et  al.  [Ca-83]  is  most  likely 
not  in  error.  Furthermore,  the  transverse  spin  transition  densities  used 
for  the  transition  to  the  3.563-MeV  level  yield  inelastic  MI  form  factors 
which  agree  with  the  (e,e')  data  for  q  <  1.4  fm_1  (see  Fig.  V-8),  a  range 
of  q  that  sufficiently  covers  the  q  3  0.55  fm~l  at  which  the  excitation 
function  was  measured.  However,  this  does  not  eliminate  uncertainties  in 


che  transverse  spin  transition  densities  since  the  electron  scattering 
Inelastic  Ml  form  factor  depends  upon  both  the  orbital  current  and  spin 
transition  densities  (see  Eq.  (IV-60a)).  Petrovich,  et  al.  [Pe-82] 
noted  that  the  monopole  spin  transition  densities  derived  from  both  the 
pure  LS-coupllng  and  the  Cohen-Kurath  intermediate  coupling  p-shell  wave 
functions  produced  Gamow-Teller  matrix  elements  »  182;  and  *  8Z  larger 
than  the  experimental  values.  Also,  Petrovich,  et  al.  [Pe-82]  did  not 
obtain  agreement  between  theory  and  their  6Li(p,p')6Li*(3.563  MeV)  data 
at  Ep  *  25  MeV,  with  che  theory  failing  to  reproduce  either  the  shape  or 
magnitude  of  the  angular  distribution.  This  disagreement,  however,  was 
not  attributed  to  uncertainties  in  che  spin  transition  density,  but 
Petrovich,  et  al.  suggested  that  other  reaction  processes  in  addition  to 
the  direct,  one-step  reaction  process  were  contributing  to  the 
transition.  Furthermore,  Cammarata  and  Donnelly  [Ca-76]  ,  in  their  study 
of  che  reaction  6Li(ir ,ir+)6He  near  threshold,  conclude  that  the  ratio  of 
the  orbital  current  and  spin  transition  densities  derived  from  the 
Donnelly  and  Walecka  wave  functions  is  probably  correct.  Thus,  as  is  che 
case  for  the  15.11-MeV  state  of  12C  [Mo-82],  che  disagreement  between  che 
simple  DWIA  analyses  and  the  3.563-MeV  angular  distributions  and 
excitation  function  is  difficult  to  explain  as  due  to  uncertainties  in 
the  transverse  spin  transition  density  or  the  spin-dependent  piece  of  the 
*N  Interaction. 

As  explained  above  and  in  Chap.  I,  simple  DWIA  calculations  using 
a  well-tested  spin-dependent  component  of  the  sN  Interaction  and 
transverse  spin  transition  densities  obtained  from  fits  to  (e,e'>  data 
fail  to  reproduce  both  the  shapes  and  magnitudes  of  the  3.563-  and 
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15.11-MeV  excitation  functions.  Furthermore,  the  energy-dependent 
anomalous  shapes  of  both  excitation  functions  are  similar.  Fig.  V-12 
shows  the  ratio  of  the  experimental  differential  cross  sections  to  the 
calculated  differential  cross  sections  for  the  excitation  functions  of 
the  two  unnatural-parity  transitions.  (The  calculation  employing  the 
phenomenological  transverse  spin  transition  density  is  used  for  the  6Li 
ratio).  As  mentioned  in  Chap.  Ill,  one  cannot  directly  compare  the 
3.563-MeV  excitation  function  to  the  15.11-MeV  excitation  function,  which 
is  clearly  evident  from  Figs.  V-10  and  V— 11.  From  Eq.  (IV-57a),  the  DWIA 
expression  for  *A  inelastic  scattering  differential  cross  sections  for 
the  excitation  functions  for  unnatural-parity  transitions  can  be  written 
as  the  product  of  two  energy-dependent  factors,  the  spin-dependent  piece 
of  the  irN  coupling  and  the  distorted  nuclear  form  factor  at  a  fixed 
momentum  transfer,  and  sin2(0)  (Si-81],  For  the  transition  to  the 
15.11-MeV  state  of  l2C,  the  strong  energy  dependence  of  the  spin- 
dependent  piece  of  the  *N  coupling  is  cancelled  by  the  energy  dependence 
of  the  distorted  nuclear  form  factor  [Si-81].  Thus,  within  the  framework 
of  the  DWIA,  the  energy  dependence  of  the  15.11-Mev  excitation  function 
follows  sin2(0).  However,  for  the  transition  to  the  3.563-MeV  state  of 
6Li,  the  distortions  are  not  as  great,  and  the  energy  dependence  of  the 
spin-dependent  piece  of  the  *N  coupling  dominates.  Therefore,  using  the 
DWIA,  Instead  of  having  a  gradually  decreasing  energy  dependence,  the 
3.563-MeV  excitation  function  is  rounded  and  peaked.  Taking  the  ratio  of 
the  experimental  to  theoretical  excitation  functions  eliminates  the 
difference  of  the  effects  of  the  distortions  and  allows  direct  comparison 
of  the  shapes  of  the  3.563-  and  15.11-MeV  excitation  functions.  From 
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Fig.  V-12:  Ratio  of  Che  experimental  differential  cross  sections  to  the 
calculated  differential  cross  sections  for  the  excitation  functions  of 
the  two  AS-&T-1  transitions  to  the  3.563-MeV  state  of  6Li  (solid  circles) 
and  the  15.11-MeV  state  of  1ZC  (solid  squares). 
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Fig.  V-12,  the  ratios  for  the  3.563-MeV  state  of  ®Li  and  the  15.11-MeV 
state  of  12C  are  very  similar  in  shape  with  the  3.563-MeV  ratio  approxi¬ 
mately  one-half  the  15.11-MeV  ratio  at  each  incident  pion  energy.  Both 
ratios  deviate  significantly  from  one  and  possess  an  energy  dependence 
resembling  the  basic  *N  interaction. 

Since  the  measured  excitation  functions  for  the  unnatural-parity 
transitions  to  the  3.563-MeV  state  of  6Li  and  the  15.11-MeV  state  of  l2C 
disagree  in  the  same  manner  with  the  excitation  functions  predicted  using 
the  DWIA,  and  uncertainties  in  the  transverse  spin  transition  densities 
or  the  spin-dependent  piece  of  the  xN  interaction  are  not  likely 
explanations  for  the  disagreement,  perhaps  another  process,  in  addition 
to  a  one-step,  impulse-approximation  mechanism,  contributes  to  both 
isovector  transitions.  A  possible  mechanism  is  the  direct  excitation  of 
A-h  components  of  the  excited  state  wave  functions,  as  proposed  in 
[Mo-82],  and  a  schematic  diagram  of  the  process  is  shown  in  Fig.  V-13 
[Mo-84a],  The  diagram  represents  the  process  in  which  the  incident  pion 
excites  a  nucleon  to  a  A  particle;  the  A  particle  decays  into  a  A 
particle  plus  a  pion,  leaving  the  nucleus  in  an  excited  A-h 
configuration.  Since  this  mechanism  involves  the  intermediate  formation 


of  a  A  particle,  its  contribution  to  the  xA  inelastic  scattering 
differential  cross  sections  should  be  resonant  [Mo-84a],  Furthermore, 
the  isospin  of  the  A  particle  (T»3/2)  and  the  isospin  of  the  nucleon  hole 
(T-l/2)  couple  only  to  a  total  T-l  or  T-2;  thus,  resonant  A-h  production 
cannot  contribute  to  excitation  of  the  T-0,  12.71-MeV  state  of  12C. 


Within  the  A-h  model  interpretation,  one  can  roughly  estimate  the 
amount  of  A-h  admixture  in  the  3.563-MeV  wave  function.  Using  the  peak 
of  the  3.563-MeV  excitation  function  (T^  -  190  MeV),  the  DWIA  calculation 
which  used  the  phenomenological  transverse  spin  transition  density,  and 
the  procedure  of  [Mo-82 ,Mo-84a]  (see  Appendix  B),  the  estimated  range  of 
B  (assumed  to  be  real),  the  probability  amplitude  for  the  A-h  component 
of  the  3.563-MeV  state,  is  0.01  <  0  <  0.13.  This  range  of  values  of  0  is 
similar  to  that  estimated  for  the  15.11-MeV  state  of  12C, 
0.026  <  8  <  0.096  [Mo-84a].  However,  if  the  direct  excitation  of  A-h 
components  is  applicable  to  both  AS*AT«1  transitions,  comparison  of  the 
two  excitation  functions  and  DWIA  calculations  (see  Fig.  V-12)  indicates 
that  the  resonant  A-h  scattering  amplitude  m*y  Interfere  differently  with 
the  p-h  scattering  amplitude  for  the  two  transitions. 


VI.  SUMMARY 


Using  Che  EPICS  system  at  Che  Clinton  P.  Anderson  Meson  Physics 
FacillCy,  it*  Inelastic  scattering  data  have  been  obtained  for  the 
transitions  to  the  3+,  T-0,  2.185-MeV,  0+,  T-l,  3.563-MeV,  and  2*,  T-0, 
4.25-MeV  states  of  6Li.  Simple  DWIA  calculations,  involving  only  the 
spin-independent  central  component  of  the  xN  interaction  and  an  isoscalar 
matter  transition  density  derived  from  Cohen-Kurath  intermediate  coupling 
p-shell  wave  functions  with  the  shape  and  strength  parameters  fixed  from 
(e,e')  data  (transition  to  the  2.185-MeV  state  only),  reproduce  well  both 
the  120-  and  180-MeV  (x+,it+')  angular  distributions  and  the  constant-q 
(q  *  109  MeV/c)  differential  cross  sections  at  -  100  to  260  MeV  for 
the  natural-parity  transitions  to  the  2.185-  and  4.25-MeV  levels.  The 
agreement  between  the  DWIA  predictions  and  the  (tt+,it+')  data  for  these 
two  transitions  is  continuing  confirmation  of  the  adequacy  of  the  DWIA 
description  for  pion-induced  transitions  to  nuclear  states  which  are 
strongly  excited. 

DWIA  calculations  using  a  well-tested  spin-independent  component 
of  the  xN  interaction  and  three  different  transverse  spin  transition 
densities  obtained  from  fits  to  (e,e')  data  fail  to  reproduce  the 
(*+,*+')  120-  and  180-MeV  angular  distributions  and  the  (it+,ir*') 
excitation  function  at  q  =  109  MeV/c  for  the  unnatural-parity  transition 
to  the  weakly  excited  3.563-MeV  state.  The  measured  excitation  function 
exhibits  an  energy-dependent  enhancement  near  incident  pion  energies  of 
T_  3  190  MeV,  with  the  difference  between  the  data  and  DWIA  calculations 
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very  similar  Co  Che  difference  observed  for  Che  excicacion  function  of 
Che  unnatural-parity  CransiCion  co  Che  15.11-MeV  state  of  l2C. 

As  for  the  15.11-MeV  excitation  function  [Mo-82],  uncertainties 
in  the  transverse  spin  transition  density  or  inadequacies  in  the  spin- 
dependenc  piece  of  Che  irN  interaction  are  an  unlikely  explanation  for  Che 
disagreement  between  Che  simple  DWIA  analyses  and  (x  +  ,it+')  data  for  Che 
transition  to  the  3.563-MeV  level.  Rather,  a  more  likely  explanation  is 
chat  an  additional  process  other  chan  a  one-step,  impulse-approximation 
mechanism  contributes  to  the  AS“AT”1  transition.  A  possible  mechanism  is 
resonant  A-h  production  as  was  proposed  in  [Mo-82]  to  explain  the 
anomalous  excitation  function  for  the  transition  to  the  15.11-MeV  state 
of  l2C.  Within  the  A-h  model  interpretation,  the  peak  of  the  3.563-MeV 
excitation  function  is  reproduced  with  an  estimated  probability  amplitude 
for  the  A-h  component  of  the  3.563-MeV  state  of  0.01  <  8  <  0.13,  a  range 
of  values  of  6  consistent  with  the  estimation  of  0  needed  to  account  for 
the  peak  of  the  15.11-MeV  excitation  function,  0.026  c  B  <  0.096 
[Mo-84a],  However,  the  A-h  and  p-h  scattering  amplitudes  may  interfere 
differently  for  the  two  transitions.  Whether  or  not  the  A-h  model  is  the 
correct  explanation  for  the  anomalous  3.563-MeV  excitation  function,  the 
(ir+,s+')  measurements  for  the  transition  to  the  3.563-MeV  state  of  6Li 
provide  another  example  of  the  failure  of  a  simple  DWIA  description  of 
transitions  to  weakly  excited  nuclear  levels. 


APPENDIX  A 


Data  Tabulation 


TABLE  III— 1 :  Differential  cross  sections  for  *+  elastic  scattering, 
6Li(»+,*+). 


T„ 

9lab 

®c.m. 

q 

d°/d«c.«. 

(MeV) 

(deg) 

(deg) 

(fm“l) 

(mb/sr) 

100 

33.0 

34.3 

0.559 

27.2 

t  0.1 

120 

15.0 

15.7 

0.289 

80.5 

t  0.4 

120 

21.0 

21.9 

0.404 

74.1 

t  0.4 

120 

24.0 

25.1 

0.460 

64.0 

i  0.3 

120 

27.0 

28.2 

0.517 

49.1 

i  0.2 

120 

27.0 

28.2 

0.517 

48.8 

t  0.2 

120 

29.0 

30.3 

0.554 

43.0 

i  0.3 

120 

32.0 

33.4 

0.609 

34.9 

t  0.2 

120 

32.0 

33.4 

0.609 

35.1 

±  0.2 

120 

42.0 

43.7 

0.790 

15.0 

i  0.1 

120 

47.0 

48.9 

0.878 

8.2 

i  0.1 

140 

26.0 

27.2 

0.551 

64.9 

i  0.3 

160 

24.0 

25.2 

0.557 

80.5 

t  0.3 

170 

23.5 

24.7 

0.569 

85.4 

*  0.5 

180 

20.0 

21.1 

n  .505 

103.4 

i  0.4 

180 

23.0 

24.3 

0.580 

79.4 

±  0.3 

180 

23.0 

24.3 

0.580 

79.5 

±  0.5 

180 

23.0 

24.3 

0.580 

79.3 

t  0.4 

180 

23.0 

24.3 

0.580 

80.9 

i  0.5 

180 

26.0 

27.4 

0.654 

61.3 

t  0.3 

180 

29.0 

0.727 

44.7 

±  0.2 

20.5 

21.7 

0.557 

100.7 

*  0.7 

215 

19.5 

20.7 

0.558 

107.1 

t  0.5 

18.5 

19.7 

0.557 

105.7 

±  0.7 

230 

18.5 

19.7 

0.557 

103.4 

±  1.5 

260 

17.0 

18.2 

0.560 

109.8 

t  0.8 
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TABLE  III-2:  Differential  cross  sections  for  *+  inelastic  scattering  to 
the  3+,  T-0,  2.185-MeV  state  of  6Li. 


00 

33.0 

34.3 

0.555 

0.40  t  0.02 

“  -- 

20 

15.0 

15.7 

0.288 

0.27  ±  0.03 

20 

21.0 

21.9 

0.401 

0.31  i  0.04 

20 

24.0 

25.1 

0.458 

0.57  ±  0.03 

20 

27.0 

28.2 

0.514 

0.63  ±  0.03 

20 

27.0 

28.2 

0.514 

0.59  ±  0.03 

20 

29.0 

30.3 

0.551 

0.74  t  0.04 

20 

32.0 

33.4 

0.606 

0.80  t  0.03 

20 

32.0 

33.4 

0.606 

0.85  ±  0.04 

20 

42.0 

43.7 

0.786 

0.97  t  0.03 

20 

47.0 

48.9 

0.873 

0.98  i  0.03 

40 

26.0 

27.2 

0.548 

1.11  t  0.05 

60 

24.0 

25.2 

0.555 

1.36  i  0.06 

70 

23.5 

24,7 

0.567 

1.78  ±  0.09 

•s.  ■  .-« 

80 

20.0 

21.1 

0.503 

1.25  t  0.07 

80 

23.0 

24.2 

0.577 

1.74  ±  0.06 

80 

23.0 

24.2 

0.577 

1.62  i  0.08 

v:-; 

80 

23.0 

24.2 

0.577 

1.59  t  0.06 

80 

23.0 

24.2 

0.577 

1.68  t  0.09 

80 

26.0 

27.4 

0.651 

1.97  t  0.06 

h*- 

80 

29.0 

30.5 

0.724 

2.32  t  0.06 

00 

20.5 

21.7 

0.555 

1.96  ±  0.12 

15 

19.5 

20.7 

0.557 

1.88  ±  0.09 

30 

18.5 

19.7 

0.555 

2.03  t  0.12 

* 

18.5 

19.7 

0.555 

2.09  t  0.25 

60 

17.0 

18.2 

0.558 

2.06  ±  0.13 
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TABLE  III-3:  Differencial  cross  sections  for  *+  inelastic  scattering  to 
the  0+*  T-l,  3.563-MeV  state  of  6Li. 


i 

®lab 

®c.m. 

9 

t 

(MeV) 

(deg) 

(deg) 

(fnT1) 

(ub/sr) 

33.0 

34.3 

0.553 

39.8  t  2.2 

-•  - 

15.0 

15.7 

0.287 

33.9  ±  3.0 

1  . 

21.0 

21.9 

0.400 

43.6  t  4.6 

24.0 

25.1 

0.456 

61.5  i  5.6 

'  ■ 

27.0 

28.2 

0.512 

41.3  t  4.7 

27.0 

28.2 

0.512 

48.4  ±  4.1 

29.0 

30.3 

0.549 

57.0  t  4.3 

• 

32.0 

33.4 

0.604 

38.0  t  3.3 

L 

37.0 

38.6 

0.694 

32.8  ±  3.2 

42.0 

43.8 

0.783 

22.3  t  3.2 

47.0 

48.9 

0.870 

12.3  ±  2.3 

26.0 

27.2 

0.547 

60.7  t  4.4 

,■ 

24.0 

25.2 

0.553 

68.5  ±  5.5 

23.5 

24.7 

0.565 

70.2  ±  7.0 

< 

23.5 

24.7 

0.565 

80.7  *  9.4 

T  '  " 

20.0 

21.1 

0.502 

117.0  *  7.6 

23.0 

24.3 

0.576 

100.0  *  7.7 

23.0 

24.3 

0.576 

95.0  t  9.4 

23.0 

24.3 

0.576 

62.0  ±  9.1 

26.0 

27.4 

0.649 

55.0  ±  5.6 

29.0 

30.6 

0.722 

43.1  ±  6.3 

-  .  , 

22.0 

23.2 

0.573 

97.0  ±  12.0 

22.0 

23.2 

0.573 

89.5  t  9.5 

20.5 

21.7 

0.554 

81.7  ±  6.4 

20.5 

21.7 

0.554 

97.0  t  5.7 

19.5 

20.7 

0.556 

73.3  ±  5.6 

.  .. 

18.5 

19.7 

0.554 

55.9  t  7.3 

.. 

18.5 

19.7 

0.554 

55.6  t  6.0 

17.0 

18.2 

0.558 

33.1  ±  10.5 

.**  *** 

•  **  < 
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TABLE  III-4:  Differential  crocs  section*  for  *+  Inelastic  scattering  to 
the  2+,  T-0,  4.25-M*V  state  of  *Li. 


T* 

9l.b 

®c.a. 

q 

do/dQ 

c.a. 

CMeV) 

(deg) 

(deg) 

(fa~l) 

(ab/sr) 

100 

33.0 

34.3 

0.552 

0.15 

i 

0.01 

120 

15.0 

15.7 

0.287 

0.08 

t 

0.01 

120 

21.0 

21.9 

0.400 

0.09 

t 

0.01 

120 

24.0 

25.1 

0.456 

0.14 

i 

0.01 

120 

27.0 

28.2 

0.511 

0.26 

t 

0.01 

120 

27.0 

28.2 

0.511 

0.20 

t 

0.01 

120 

29.0 

30.3 

0.548 

0.28 

± 

0.01 

120 

32.0 

33.4 

0.603 

0.27 

t 

0.01 

120 

37.0 

38.6 

0.693 

0.  O 

± 

0.01 

120 

42.0 

43.8 

0.782 

0.34 

± 

0.01 

120 

47.0 

48.9 

0.868 

0.27 

± 

0.01 

140 

26.0 

27.2 

0.546 

0.32 

t 

0.01 

160 

24.0 

25.2 

0.553 

0.51 

t 

0.01 

170 

23.5 

24.7 

0.565 

0.61 

i 

0.02 

170 

23.5 

24.7 

0.565 

0.64 

± 

0.02 

180 

20.0 

21.1 

0.502 

0.86 

± 

0.02 

180 

23.0 

24.3 

0.575 

0.69 

t 

0.02 

180 

23.0 

24.3 

0.575 

0.80 

+ 

0.02 

180 

23.0 

24.3 

0.575 

0.76 

+ 

0.02 

180 

26.0 

27.4 

0.649 

0.83 

+ 

0.01 

180 

29.0 

30.6 

0.721 

0.88 

0.01 

190 

22.0 

23.2 

0.572 

0.87 

± 

0.03 

190 

22.0 

23.2 

0.572 

0.98 

+ 

0.03 

200 

20.5 

21.7 

0.554 

0.61 

£ 

0.01 

200 

20.5 

21.7 

0.554 

0.74 

£ 

0.01 

215 

19.5 

20.7 

0.555 

0.89 

£ 

0.01 

230 

18.5 

19.7 

0.554 

0.71 

t 

0.02 

230 

18.5 

19.7 

0.554 

0.70 

t 

0.01 

260 

17.0 

18.2 

0.557 

0.71 

£ 

0.02 

APPENDIX  B 


Estimation  of  B 

Wichln  the  A-h  model  interpretation,  the  peak  of  the  excitation 
function  measured  for  the  AS*AT“1  transition  to  the  3.563-MeV  state  of 
6Li  Is  reproduced  with  an  estimated  probability  amplitude  of  the  A-h 
component  of  the  3.563-MeV  state  of  0.01  <  $  <  0.13.  6  is  estimated  in 

the  following  manner.  The  inelastic  scattering  differential  cross 
section  is  written  as  [Mo-82 ,Mo-84a] 

1£<E,0)  -  | A(E , 6)  +  0B(E,0)|2  ,  (B-l) 

dw 

where  |A(E,0)j  is  the  p-h  scattering  amplitude,  j B( E , 0 ) |  is  the  A-h 
scattering  amplitude,  and  0  is  the  probability  amplitude  for  the  A-h 
component  of  the  3.563-MeV  state  with  respect  to  the  ground  state.  The 
3.563-MeV  excitation  function  peaks  at  T^  *  190  MeV  and  0^aj,  “  22°. 
Using  this  energy  and  angle  and  the  DWIA  analysis  employing  the 
phenomenological  transverse  spin  transition  density, 

|A(E, 0)|  -  8.2  •'ub/sr  .  (B-2) 

|B(E,0)|  may  be  roughly  estimated  by  examining  the  isospin 
dependence  of  cross  sections  for  (T^T^)  *  ^Tf,Tzf^  transitions 
Involving  only  p-h  configurations  or  A-h  configurations  of  the  final 
state  with  respect  to  the  initial  state,  and  by  assuming  a  simple  A-h 


model  for  Che  3.563-MeV  level.  The  relative  isospin  dependence  is 
determined  using  Lee  and  Lawson's  [Le-80a]  extraction  of  the  isospln 
dependence  of  the  »A  scattering  amplitude  in  the  DWIA  formalism,  and 
modifying  their  result  to  Include  a  A  in  the  final  state.  Using  P3  3 
dominance,  which  simplifies  Lee  and  Lawson's  result,  the  ratio  of  cross 
sections  for  a  (Tj-0,Tzi«0)  *  (Tf-l,Tzj-0)  transition  for  the  final  state 
being  a  A-h  configuration  to  the  final  state  being  a  p-h  configuration  is 
0.623/0.25.  To  determine  the  possible  A-h  configurations  of  the 
3.563-MeV  state  with  respect  to  the  ground  state,  an  extremely  simple 
model  is  used.  In  this  model,  the  3.563-MeV  state  is  a 
(lpt/2  *  lP3/2_1 )j«0,T»l  configuration  with  respect  to  the 
ground  state.  (This  configuration  was  chosen  from  the  dominant  one  body 
density  matrix  element  predicted  from  Cohen-Kurath  Intermediate  coupling 
p-shell  wave  functions.)  One  could  expect  the  following  A-h  configura¬ 
tions  to  mix  with  the  above  p-h  configuration:  ^ 1 p 1 / 2  x  lp3/2_1^J“0  T*l' 

( lp3/2  x  lp3/2  l\j-0,T-l»  and  ^lp5/2  X  lp3/2  l^J»0,T«*l*  Wich  Cohen" 
Kurath  Intermediate  coupling  p-shell  wave  functions,  the  inelastic 

scattering  differential  cross  section  is  quenched  by  a  factor  of  2.25 
from  the  differential  cross  section  computed  with  the  above  simple  p-h 
configuration.  Assuming  that  this  quenching  is  not  present  for  the  A-h 
configurations  and  that  the  three  transitions  Involving  a  A,  P3/2  *  p5/2’ 
p3/2  *  p3/2*  p3/2  *  P 1 / 2 »  contribute  coherently  to  the  differential  cross 
section  and  are  weighted  as  / 2 j  ^  +  1 , 


The  only  ocher  factor  which  needs  Co  be  considered  is  Che  ratio  of 


coupling  conscancs  f^As^AN**  This  raCio  can  '>e  obtained  from  the  quark 
model  and  is  *  ^25/8  [Mo-84a],  Thus, 


|  B(E, 9) | 


/6  +  /Z  +  /2  /OUT  /25 
- 72 - xTOTx7Tx 


/2.25)A(E,0)|2  . 


(B-3b) 


Substituting  in  the  value  of  )A(E,9)|  from  Eq.  (B-2), 

|B(E,9)|  -  143  /ub/sr  .  (B-3c) 

With  Eqs.  (B-l)  through  (B-3)  and  assuming  that  the  two  scattering 
amplitudes,  p-h  and  A-h,  are  completely  in  phase  or  are  completely  out  of 
phase,  0.01  <  8  <  0.13,  where  8  is  assumed  to  be  real. 
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ABSTRACT 


The  excitation  function  for  t+  inelastic  scattering  to 
the  0+,  T"l,  3.563-MeV  level  of  6ll  has  been  measured  at  a 
constant  momentum  transfer  q  *  109  MeV/c  for  Incident  plon 
energies  from  100  co  260  MeV,  Although  the  differential 
cross  sections  extracted  for  the  natural-parity  transitions 
to  the  3+,  T-0,  2. 185-MeV  and  2+,  T-0,  4.25-MeV  levels  are 
well  reproduced  within  the  framework  of  the  distorted-wave 
Impulse  approximation  (DWIA),  DWIA  calculations  fall  to 
reproduce  the  anomalous  excitation  function  observed  for 
the  transition  to  the  3.563-MeV  level.  The  shape  of  the 
3.563-MeV  excitation  function  Is  similar  co  that  previously 
observed  for  x4  Inelastic  scattering  to  the  the  1  +  ,  T-l, 
15.11-MeV  state  of  12C  (C.  L.  Morris,  et  al. ,  Phys.  Lett. 

108B,  172  (1982)].  The  same  mechanism  may  be  responsible 
for  the  observed  excitation  funcdons  of  both  AS-AT-l 
transitions,  A  possible  mechanism  is  the  previously 
proposed  direct  excitation  of  &  particle-nucleon  hole  (A-h) 
components  In  the  wave  functions. 

PACS  numbers:  25.80.0j,  25.80.Fm 

(NUCLEAR  REACTIONS:  6L1(w+,w'k  );  100  <  T„  <  260  MeV;  measured  o(8); 
DWIA  analysis;  deduced  <,  the  matrix  element  between  A-h  and  p-h 
basts  states  for  the  3.563-MeV  state] 


I.  INTRODUCTION 


Comparison  of  plon-lnelastlc-scattertng  data  In  the  energy  range  of  the 
(3,3)  reaonance  with  calcuiatlone  using  the  dlstorted-^rave  lapulae 
approximation  (DWIA)  and  welldtnown  transition  densities,  Indicates  that  a 
alngle-etep,  lapulae  approximation  la  an  adequate  description  for  plon-lnduced 
transitions  to  nuclear  states  which  are  strongly  excited.  For  exaaple,  Lee  and 
Kurath1  and  Lee  and  Lawson,2  using  the  DWIA  and  transition  densities  derived 
from  shell-model  wave  functions,  qualitatively  reproduce  angular  distributions 
for  plon  Inelastic  scattering  co  excited  states  of  various  p-  and  ad-shell 
nuclei.  Morris,  et  al.3  and  Boyer,  et_al.1*  obtain  good  agreement  between 
angular  distributions  and  DWIA  calculations,  which  use  empirical  transition 
densities  determined  from  electron-scattering  data,  for  Inelastic  scattering  to 
low-lying  collective  states  In  l2C,  l*1ICa,  42Ca,  44Ca,  and  48Ca.  Furthermore, 
the  excitation  functions  measured  at  a  constant  momentum  transfer  for  the 
unnatural-parity  transitions  to  the  4”,  19.25-MeV  scate  and  2”  structure  at 
approximately  18.4  MeV  In  12C  (Ref.  5)  and  the  9/2*,  9.5-MeV  state  In  13C  (Ref. 
6)  decrease  with  Increasing  Incident  plon  energy  as  predicted  by  the  DWIA.7 
Within  the  framework  of  the  DWIA,  these  unnatural-parley  transitions  are 
particularly  simple.  Involving  only  the  spin-dependent  part  of  the  effective 
plon-mucleus  Inelastic  Interaction  and  the  transverse  spin  transition  density.8 
In  particular,  for  the  unnatural-parity  transitions  to  states  of  stretched 
configuration  In  l80  (4-,  17. 79-,  18. 98-,  and  !9.80-MeV  levels)  and  29St  (6", 
It. 58-  and  14.3A-MeV  levels),  Carr,  et  al.8  satisfactorily  reproduce  the 
measured  angular  distributions  with  DWIA  calculations  using  spin  transition 
densities  fixed  from  (e,e‘)  and  (p,p')  data. 


For  plon-lnduced  •scitatioa  of  nuclaar  atataa  thac  are  weakly  axel cad  by  a 
ona-atap,  lapulsa-approxlaatlon  aachanisa,  nuclaar  aadlua  affacts  and  aultlstep 
procaaaaa  aay  ba  Important  and  thara  aay  not  bo  good  agreaaent  between  the 
plon-lMlasclc-scatterlng  data  and  DWIA  calculatlona ,  aa  for  the  above  atrongly 
excited  exaaplas.  Thia  la  lndaad  the  caae  for  the  T«1  aaabar  of  the  weakly 
excited  1*  doublet,  12.71  MeV  (T-0)  and  15.11  MeV  (T-l),  of  l2C.10  Cohen-Kurath 
wave  functions11  daacrlbe  the  15.11-  and  12.71-HeV  atataa  aa  near  analogs,  and 
thle  description  for  the  apln  transition  densities  la  supported  by  (e,e'> 
data.12  Calculations  using  a  ona-atap,  DWIA  aachanisa  predict  the  ratio  of 
cross  sections  for  unnatural-parity  transitions  to  members  of  an  analog- 
antianalog  pair  to  ba  four  to  one.8’10  Both  the  constant-aomentua-transfer- 
axcltatlon  function  and  the  angular  distributions  for  the  12.71-MeV  state  are 
wall  reproduced  by  OUIA  calculations.  However,  the  measured  ratio  of  four 
tlaas  the  averaged  **  and  »”  differential  cross  sacclons  for  the  15.11-MeV 
state  to  the  averaged  **  and  *“  differential  cross  sections  for  the  12.71-HeV 
state  deviates  significantly  from  one,  especially  at  energies  near  180  MeV 
(ratio  is  approximately  three),  and  displays  a  rapidly  varying  energy 
dependence.  (Averaging  the  **  and  differential  cross  sections  removes  the 
effect  of  lsospin  mixing  between  the  two  states  on  the  ratio  to  better  than 
12. l0)  Also,  DWIA  calculations  do  not  agree  with  the  15.11-MeV  angular 
distributions  at  energies  near  the  (3,3)  resonance,12 

Uncertainties  In  the  spin  transition  density  or  Inadequacies  In  the  spin- 
dependent  piece  of  the  effective  pion-nucleus  Interaction  are  an  unlikely 
explanation10  for  the  anomalous  excitation  function  for  the  AS»iT*l  transition 
to  the  15.11-MeV  level  of  12C.  Rather,  a  more  likely  explanation  Is  that  an 
additional  process  other  chan  a  one-step,  Impulse-approximation  mechanism  Is 
contributing  to  Che  lsoveccor  transition.  Therefore,  to  Investigate  further 
plon-lnduced  excitation  of  weakly  excited  nuclear  levels,  we  consider  the 
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AS“AT"1  transition  from  the  1+,  T«0,  ground  state  of  6 11  to  the  O'4-,  T*l, 
3.563-MeV  level.  We  have  measured  an  excitation  function  for  t  +  Inelastic 
scattering  to  the  3.563-MeV  state  of  6U.  at  Incident  pton  energies  from  100  to 
260  MeV  and  partial  angular  distributions  at  120  and  180  HeV.  Ue  also  present 
In  Table  I  the  differential  cross  sections  for  elastic  scattering  and  for 
Inelastic  scattering  to  the  2. 185—  and  4.25^leV  states.  The  2.185-  and 
4.25-HeV  differential  cross  sections  are  coapared  to  alcroscoplc  0W1A 
calculations  using  transition  densities  derived  froa  Cohen-Kurath  wave 
functions.11  The  3.563-MeV  differential  cross  sections  are  coapared  to 
microscopic  DWIA  calculations  using  spin  transition  densities  derived  froa 
Cohen-Kurath  wave  functions,11  the  empirical  shell-model  wave  functions  of 
Donnelly  and  Walecka,12  and  the  phenomenological  wave  functions  of  bergstrom, 
et  at. 


tl.  DATA  ACQUISITION  AND  REDUCTION 

The  data  were  collected  using  the  Energetic  Plon  Channel  and  Spectrometer 
(EPICS)  system1*  at  the  Clinton  P.  Anderson  Meson  Physics  Faclllcy  (LAMPF). 
Four  separate  targets  consisting  of  sheets  of  enriched  lithium,  >18?  6ll, 
fabricated  by  the  Oak  Ridge  National  laboratory  were  used  during  the 
experiment.  Two  of  the  targets,  with  dimensions  of  22.9  cm  x  15.2  cm  and  areal 
densities  of  202  mg/cm2  and  100  rag/cm2 ,  were  used  for  approximately  one-half  of 
Che  total  data  acquisition  and  contained  no  discernible  contamination.  The 
remaining  two  targets,  with  dimensions  of  20  cm  «  10  cm  and  areal  densities  of 
205  mg/cm2  and  95  mg/cm2 ,  were  contaminated  by  exposure  to  the  air  during 
shipping.  Hydrogen  contamination  was  negligible  In  the  thick  20  cm  «  10  cm 
target  (Alt)  but  appreciable  In  the  thin  20  cm  x  10  cm  target  (*5T).  Stnce  we 


clearly  observed  the  6. 1 3—  and  6.92-MeV  states  of  l60  In  sooe  *+  Q-value 
spectra  but  never  observed  any  Inelastic  peaks  of  l1*?),  we  assumed  the  heavy 
Impurity  to  be  160.  With  this  assumption,  the  total  amounts  of  Impurities  for 
the  thick  and  thin  20  cm  x  10  cm  targets  were  2.3Z  and  18Z,  respectively.  A 
comparison  of  differential  cross  sections  for  the  elastic  scattering  from  5Li, 
extracted  from  the  contaminated  and  uneontamlnated  targets  at  the  same  Incident 
plon  energy  and  scattering  angle.  Indicate  that  the  estimated  amounts  of 
Impurities  are  accurate  within  <11.  ; LI  was  not  detected  In  any  of  the 
targets. 

Fig.  I  shows  the  Q-value  spectrum  for  6U(«'t‘,t+’ )  using  the  205  mg/cm2 ,  20 
cm  x  10  cm  target  for  ■  140  MeV  and  S^g  «  26°.  This  spectrum  Is 
representative  of  Che  spectra  used  for  determination  of  cross  sections  for 
excitation  of  the  3.563-  and  4.25-+feV  levels.  Although  the  experimental 
resolution  was  not  Che  best  for  the  205  mg/cm2  target  (240  keV  full  width  at 
half  maximum  (FWHN)),  the  3.  563-MeV  level  Is  clearly  resolved  from  the  large 
conctnuum  background  and  the  4.25-MeV  level.  The  discontinuity  In  the  Q-value 
spectrum  between  the  ground  state  and  2.185-MeV  stace  Is  an  artifact  of  the 
hardware  veto  of  elastic  counts  used  Co  limit  data  rates  during  data 
acquisition  by  rejecting  the  majority  of  events  with  a  0  value  below  *1.2  MeV. 
The  elastic  and  2.185-MeV  cross  sections  were  determined  without  the  use  of 
this  hardware  veto  (see  the  Insert  In  Fig,  t). 

Peak  areas  were  extracted  from  the  Q-value  spectra  using  the  computer 
program  LOAF.1'  The  backgrounds  under  the  elastic  and  3+  peaks  were  fitted  with 
a  first-order  polynomial,  while  the  backgrounds  under  the  0+  and  2*  peaks  were 
fitted  with  a  third-order  polynomial.  Line  shapes  for  the  elastic  and  3+  peaks 
were  extracted  from  the  spectra  for  each  target  at  each  Incident  plon  energy 
because  the  resolution  varied  with  energy  from  200  keV  (220  keV)  FVHM  at 
Tf  «  100  MeV  to  260  keV  (260  keV)  FWKM  at  T„  -  260  MeV  for  the  thin  (thick) 
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targeta.  Hit  lint  shape  of  Cht  0+  ptak  was  caktn  Co  bt  eht  Mat  at  Cht  3+, 
while  cht  lint  thapt  of  cht  2+  wat  conecrucced  by  folding  cht  lint  ahapt  uttd 
for  Cht  3*  wlch  a  Lorentzlan  of  680  keV  FWHM.  The  peak  poeldona  were 
contcralned  co  bt  0.0,  2,  185, 19  3. 563, 19  and  4, 25  MeV.  BoCh  Che  poaldon  and 
nacural  vldch  of  Che  2'*'  level,  4,25  t  0.02  MeV  and  680  t  20  keV,  were 
decernlned  from  flcs  co  several  spacers  where  Chit  scace  was  predominant. 
Previous  determinations  from  ocher  experiments  are  4.27  t  0.04  MeV  and 
690  1  120  keV,19  4.29  t  0.02  MeV  and  850  t  SO  keV,29  and  4.30  t  0.01  MeV  and 
480  *  80  keV.20  The  5.37-  and  5.65-MeV  seaces  were  not  observed  In  the  sptccra, 
and  thus  no  attempt  was  made  to  Include  them  In  the  fits. 

The  consistency  of  the  fits  to  the  various  Q-value  spectra  was  checked  by 
extracting  areas  for  100  keV  wide  segments  of  background  centered  about  3.563 
and  4.25  MeV.  Plots  of  che  background  yields  for  both  the  3.563-  and  4.25-MeV 
segments  resulted  In  smooch  and  continuous  angular  distributions  at  120  and  180 
MeV.  Also  the  background  yields  versus  Incident  pton  energy  are  smooth  and 
continually  Increasing.  Ue  therefore  have  indication  that  our  extraction  of 
the  peak  areas  from  the  large  continuum  background  is  consistent  at  different 
incident  pton  energies  and  scattering  angles. 

Experimental  yields  were  measured  for  >  25°  by  monitoring  the  EPICS 
channel  beam  flux  with  an  Ionization  chamber  located  downstream  from  the 
scattering  target.  For  <  25°  the  ionization  chamber  was  not  used  since  it 
partially  blocked  the  spectrometer  entrance.  For  these  angles,  relative 
normalization  was  accomplished  through  an  Ionization  chamber  located  within  the 
pton  production  target  cell  and  a  charge  Integrating  toroidal  coll  located 
upstream  of  the  plon  production  target.  Monitoring  of  the  ratio  of  the 
ionization  chamber  current  to  the  currents  from  these  two  monitors  of  the 
proton  beam  showed  <2T  fluctuations  for  >  25°,  establishing  them  as 
reliable  beam  flux  monitors.  Absolute  cross  sections  were  calculated  by 


normalizing  to  »+  scattering  on  hydrogen  (CH^  targats  of  dimensions  22.9 
cm  *  15.2  cm  and  20  cm  *  10  cm  corraapondlng  to  tha  two  dlffarant  alzaa  of  *1.1 
targata)  ualng  Coulomb-corractad  phaaa-shlft  pradletlona  from  the  computer  coda 
CROSS21  with  the  phaaa  ahlfta  of  Rowe,  Salomon,  and  landau.22 

The  data  ware  corrected  for  computer  live  time,  multiwire  proportional 
drift  chamber  efficiency,  plon  aurvlval  fraction  through  the  apectrometer,  and 
the  variation  of  the  spectrometer's  solid  angle  with  plon  momentum.  The  quoted 
error  bars  are  statistical  only.  Total  systematic  errors  are  estimated  to  be 
*7Z  due  to  uncertainties  of  ±32  In  chamber  efficiency,  ±32  in  plon  survival 
fraction,  *22  In  the  spectrometer's  solid  angle  variation  with  momentum  in  the 
spectrometer,  * 32  In  channel  beam  monitoring,  and  *32  In  normalization  to  v  + 
scattering  on  hydrogen.  Furthermore,  the  data  for  the  3.563-  and  4.25-MsV 
states  contain  additional  systematic  errors  of  ±152  and  *102,  respectively,  due 
to  the  uncertainly  In  the  fitting  of  the  large  continuum  background  and  the 
uncertainty  In  the  position  and  width  of  the  2*  state.  These  systematic 
errors  were  Inferred  by  varying  the  order  of  the  polynomial  fit  to  the 
background,  and  by  varying  the  position  and  width  of  the  2+  state  from  4.23  to 
4.27  MeV  and  660  to  700  keV. 

IV.  DATA  ANALYSIS 

A.  Elastic  scattering  analysis 

The  first-order,  zero-range.  Impulse-approximation  elastic  calculations 
were  performed  with  a  variation23  of  the  coordinate-space  computer  program 
PIRK,24  which  solves  a  Klaln-Cordon  equation,  using  only  linear  terms  In  the 
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optical-model  potential.  For  all  6L1  alaatlc  calculations,  us  used  tha 
Klasllnger  form2*  of  the  optical-model  potential  given  by 

V(r)  -  -Akjb0p(r)  Ab^pCr)-?  ,  (1) 
where  p(r)  la  the  nucleon  deneley  normalized  to  unity,  kf  Is  the  lab  momentum 
of  the  Incident  plon,  and  A  Is  the  mass  of  the  target.  The  complex  b^  and  bj 
coefficients  are  constructed  from  the  pton^iucleon  phase  shifts  of  Rowe, 
Salomon,  and  Landau22  evaluated  at  an  energy  of  30.0  MeV  below  the  incident 
plon  beam  energy.  This  procedure  has  been  demonstrated  by  Cottlngane  and 
Holtkamp23  to  give  better  agreement  to  plon^laatlc-ecatterlng  data  for  nuclei 
ranging  from  *Be  co  208Pb.  This  phenomenological  shift  In  the  collision  energy 
Is  a  method  for  adjusting  the  plon-nucleus  kinematics  so  as  optimally  to 
factorise  the  optical-model  potential  (see  Ref.  12  for  a  discussion  of  the 
above  procedure).  The  proton  matter  density  distribution  was  characterized  by 
a  three-parameter  phenomenological  distribution  taken  from  elastic  electron 
scattering26  with  the  finite  size  of  the  proton  charge  removed.  This 
distribution  has  the  form 

p(r)  «  Z/(8*^',*)[(l/'a3)exp(-r2/4a2)  - 

(cZfbbi-rilMbMexpf-r2/^2)!  ,  (2) 
with  a  «  U.928  fm,  b  •  1.26  fm,  and  c  »  9.48  fm.26  This  phenomenological 
distribution  was  also  used  for  the  neutron  matcer  den9lcv  distribution. 

The  present  120-  and  180-+leV  x  +  elastic  data,  the  SIN  164-MeV  x  ”  elastic 
data,22  and  elastic  calculations  are  presented  In  Fig.  2.  Since  the  elastic 
optical  potential  Is  used  for  the  generation  of  the  distorted  waves  for 
Inelastic  calculations,  the  good  agreement  Indicates  adequate  handling  of  the 
distortions. 


8.  Inelastic  scattering  analysis 
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The  Inelastic  calculation*  presented  In  this  papar  art  Born  approximation, 
DHIA  calculation*  In  which  tha  plon-nuclaua  eranaltlon  amplitude  la  a 
configuration-space,  folded  product  of  a  distortion  function  and  fora  factor. 
The  distortion  fuactloa,  a  product  of  Initial  and  final  plon  dlstortad  waves. 
Is  coaputed  f roe  tha  elastic  optica l-node 1  potential  employed  la  tha  elastic 
calculations.  Calculation  of  the  fora  factor,  whose  specific  fora  depends  upon 
the  Inelastic  transition  and  Is  the  folded  product  of  the  plon-nucleon 
Interaction  and  the  nuclear  transition  density,  uses  the  frozen-target  approxl¬ 
aatlon  aasualng  on-shell  kinematics  and  a  collision  energy  obtained  froa  the 
elastic  calculations.  For  Inelastic  calculations  using  pure  LS-coupllng  and 
Cohcn-Kurath  intermediate  coupling  wave  functions,11  we  use  a  modification  of 
the  generalized  inelastic  scattering  potential  code  AU.WRLD28  to  calculate  the 
fora  factor.  The  harmonic  oscillator  parameters  and  renormalization  constants 
include  the  standard  center-of-aass  correction  needed  when  using  shell-model 
wave  functions.  The  distortion  function  and  differential  cross  sections  are 
then  generated  from  the  cod*  UTDWPI.29  For  all  other  inelastic  calculations  w* 
use  only  the  code  UTDWPI  and  no  center-of-mass  correction  is  included  in  the 
harmonic  oscillator  parameters  and  renormalization  constants.  However,  these 
inelastic  calculations  use  center-of-mass  corrected  form  factors  with  the 
correction  being  applied  In  momentum  space. 


-:i 


1.  3*.  T-0,  2. 185-MeV  state 


Electron-scattering  daca  indicate  that  the  natural-parity  transition  to 
the  2.185-M*V  level  of  6ll  is  almost  completely  longitudinal ,  with  measurements 
in  the  region  of  q  •  0.7  to  1.8  fm-1  yielding  a  transverse  form  factor  which  is 
less  than  22  of  the  longitudinal  form  factor.30  Since  che  spin-orbit  transition 
density  Is  approximately  zero,  the  form  factor  Involves  only  the  central 
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compooent  of  the  ploa-audooa  interaction  and  tha  lioacaiar  aattar  cranalton 
density.31  Inelastic  calculations  for  this  natural-parity  transition  have  been 
performed  using  two  different  transition  densities  derived  froa  pure  LS- 
coupllng  and  Cohen-Kurath  Intermediate  coupling  p-shell  wave  functions.  Each 
calculation  aaployed  equal  transition  densitlea  for  the  protons  and  neutrons, 
the  haraonic  oscillator  paraaater.  a.  and  renoroalltatlon  constant  for  the 
transition  density  derived  froa  the  pure  LS-coupllng  p-shell  wave  functions  are 
0.534  fa-1  and  2.03.  These  values  were  required  to  fit  the  lnelastle-elactron- 
scatterlng  data  considered  in  Ref.  32  and  give  good  agreeaent  between  theory 
and  experiment  for  tha  6U.(p,p')6Li*(2.185  MeV)  differencial  cross  sections  at 
Ep  »  25  and  45  MeV.32  The  transition  density  derived  froa  tha  lnteraedlate 
coupling  p-shell  wave  functions  uses  a  •  0.558  fa-1  and  a  renormalization 
constant  of  1.93,  again  determined  froa  (e,e')  data.33  Electron-scattaring 
longitudinal  form  factors  from  Refs.  14  and  34  are  shown  In  Fig.  3.  The  solid 
(dashed)  curves  correspond  to  che  transition  densities  computed  froa  the  pure 
IS-  (intermediate)  coupling  p-shell  wave  functions.  both  theoretical  form 
factors  are  similar,  with  the  intermediate  coupling  form  factor  in  better 
agreeaent  with  che  electron-scattering  data.  The  **  2.185-MeV  data  and 
calculations  for  Tff  ■  120  and  180  MeV  are  presented  in  Fig.  4.  Fig.  5  shows 
the  data  and  calculations  for  T„  •  100  to  260  MeV  with  che  differential  cross 
sections  corresponding  to  a  constant  momentum  transfer  q  *  109  MeV/c.  From  the 
data  and  calculations,  che  first  maxima  of  che  angular  distributions  for  the 
2.185-MeV  state  are  expected  to  be  at  q  *  164  MeV/c.  The  Intermediate  coupling 
calculation  is  in  good  agreement  with  both  che  120-  and  the  180-MeV 
experimental  angular  distributions.  Furthermore,  this  calculation  reproduces 
well  che  constanc-q  experimental  differential  cross  secclons,  considering  that 
the  theoretical  values  plotted  in  Fig,  5  are  taken  from  che  steep  forward  slope 
of  che  verlous  angular  distributions  where  errors  would  produce  the  greetest 
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varlatlons.  Tha  pur  a  LS -coup ling  calculation  yields  similar  ahapaa  for  tha 
120-  and  180^!aV  angular  distributions  and  tha  constant-^  diffarantlal  croaa 
sectlona  aa  doaa  tha  Intermediate  coupling  calculation  but  overestimates  tha 
magnitudes.  Such  disagreement  suggests  chat  the  renormalization  constanc  la 
too  large. 

2.  2+,  T«0,  4.25-W«V  state 

The  electron-scattering  form  factors  for  both  natural-parity  transitions 
to  the  2.185-  and  4.25-MeV  levels  have  the  same  dependences  on  the  momentum 
transfer.35  Therefore,  the  transition  to  the  4.25-MeV  state  Is  principally 
longitudinal,  and  we  treat  this  transition  aa  completely  analogous  to  the 
excitation  of  the  2. 185-MaV  level.  However,  even  though  the  4.25-MeV  state  Is 
observed  In  (e,e')  spectra,  a  detailed  form  factor  has  not  been  measured 
because  of  this  state's  large  natural  width,  the  overlapping  5.37-MeV  state, 
and  the  large  continuum  background.36  We  therefore  cannot  deduce  a  transition 
density  from  electron-scattering  data  as  was  the  case  for  the  2. 185-MeV 
transition.  Thus,  we  simply  used  a  transition  density  derived  from  Cohen- 
Kurath  Intermediate  coupling  p-shell  wave  functions  with  a  *  0.52  fm”l  and  a 
renormalization  constant  of  0.82  determined  from  fitting  only  our  (*+,s+'>  data 
at  T^  •  120  and  180  MeV.  This  transition  density  yields  a  radiative  width, 
of  3.02  eV  In  agreement  with  the  experimental  value  of  5.4  t  2.8  eV.19 
The  4.25-MeV  experimental  and  theoretical  120-  and  180-MeV  angular 
distributions  and  the  constant-d  differential  cross  sections  are  shown  In  Figs. 
6  and  7,  respectively.  The  agreement  Is  not  as  good  as  for  the  2.185-HeV  state 
hut  Is  very  reasonable  considering  tha  difficulties  In  extracting  the  cross 
sections  and  possible  uncertainties  In  the  transition  density. 


3.  0*,  T-t,  3.563-MaV  state 


As  Is  the  case  for  the  natural-parity  transitions  to  the  2.185-  and 
4.25-fleV  levels  of  6Li,  the  Born  approximation  DWIA  description  of  the  plon- 
induced  unnatural-parity  transition  to  the  3.563-MeV  state  Is  straightforward. 
Involving  only  a  single  component  of  the  plon-nucleon  Interaction  and  a  single 
nuclear  transition  density.  Using  a  single  scattering  model,  the  spin- 
dependent  piece  of  the  Interaction  Is  represented  by  the  zero-range  spin-orbit 
operator,4  and  the  transition  density  Is  the  transverse  spin  transition 
density.  Ue  employed  three  different  transverse  spin  transition  densities 
determined  from  fits  to  various  (e,e‘)  measurements  but  compared  to  the  most 
recent  sets  of  (e,e')  data  (Refs.  14  and  37).  The  proton  and  neutron 
transition  densities  were  taken  to  be  equal.  The  reasons  for  using  three 
transition  densities  are:  (1)  The  p-shell  harmonic  oscillator  basis  Is  known 
to  result  In  a  poor  description  of  the  3.563-MeV  electron-scacterlng  transverse 
form  factor  over  the  entire  second  lobe,  q  >  1.4  fm-1 .  Ue  therefore  employ  one 
transition  density  which  reproduces  boch  the  first  and  second  lobes  of  the  form 
factor.  (2)  Since  there  Is  no  analog  to  the  3.563-MeV  state  as  Is  the 
situation  for  the  12. 71-/15. 1 1-MeV  analog-antlanalog  pair  of  l2C,  any  anomalies 
In  the  3.563-MeV  excitation  function  depend  upon  comparison  of  theory  with  data 
and  not  data  for  analogs  as  for  the  l5.U-MeV  level.  Thus,  accurate  spin 
transition  densities  are  Imperative. 

The  first  analysis  for  the  spin-flip  transition  to  the  3.563-MeV  state 
used  a  spin  transition  density  derived  from  Cohen-Kurath  intermediate  coupling 
p-shell  wave  functions.  The  harmonic  oscillator  parameter  was  chosen  to  he 
0.518  fm-1  from  the  work  of  Petrovich,  et  al. , who  fitted  low  q  electron¬ 
scattering  data  using  pure  LS-coupllng  p-shell  wave  functions.  We  used  a 
renormalization  constant  of  0.97,  Petrovich,  et  al.  conclude  that  both  secs 


of  chaorecical  wav*  function*  provide  *n  adequate  deecrlption  of  the 
experimental  etaclc  aoaant*  and  cranaltlon  probabilities  of  *11  with  the 
exception  of  the  quadrupole  nonane.32  However,  wa  chose  the  intermediate 
coupling  p-ehall  wav*  functions  becausa  they  glv*  a  slightly  batter  fit  to  the 
•oat  recent  eat*  of  (*,*')  data  (Refs.  14  end  37)  for  the  Inelastic  HI  fora 
factor.  Fig.  8  show*  th*  calculated  transverse  fora  factor  (solid  curve)  and 
(*,*')  data  for  the  transition  to  eha  3.563-M*V  state. 

Th*  second  calculation  eaployad  a  spin  transition  density  obtained  froa 
the  eaplrlcal  shell-model  wave  functions  of  Donnelly  and  Walecka.13  These  wave 
functions  are  almost  Identical  to  the  pure  LS-coupllng  p-shell  wav*  function*. 
Using  p-shell  harmonic  oscillator  radial  wave  functions  for  the  valence 
nucleons,  Donnelly  and  Walecka  determined  the  on*  body  density  matrix  elements 
froa  normalization  conditions  and  froa  fits  to  th*  ground  state  magnetic  dipole 
and  electric  quadrupole  aoaents  and  the  HI  fora  factors  for  electron  elastic 
and  Inelastic  scattering  for  q  <  1.02  fa-1.13  A  harmonic  oscillator  parameter13 
a  •  0.493  fa-1  and  a  renormalization  constant  of  0.96  were  used  in  our 
calculations.  Th*  resultant  Inelastic  HI  form  factor  (dashed  curve)  is 
compared  to  the  (*,*')  data  In  Fig.  8. 

Th*  third  analysis  used  a  phenomenological  spin  transition  density  based 
on  the  work  of  Bargstroa,  et  si.1’*  This  transition  density  was  derived  In  the 
same  manner  as  th*  transition  density  of  Donnelly  and  Ualecka  except  that  we 
fitted  the  (*,*')  data  of  Refs.  14  and  37,  which  extend  to  q  •  2.96  fa"1 ,  and 
assumed  a  polynomial  form  for  the  p-shell  radial  transition  density.  The  wave 
functions  for  the  SL1  ground  state  and  3.563-MeV  state  were  taken  to  be 
described  by  the  SASK-A  amplitudes  of  Bergstrom,  et  al.1*  The  radial  transition 
density  has  a  phenomenological  form 

R2(r)  -  e-r2^2(a2r2  *  a,r-  4-  a&r6)  ,  (3) 
with  b  -  2.02  fa,  a2-  6.625  *  10'2  fa-s ,  a„-  -5.036  *  10"3  fa-7, 


1.967  *  I0~v  fa"^ ,  yielding  a  reduced  x*  of  1.32  from  a  fit  to  the 
tnelaeclc  Ml  fora  factor,  which  la  presented  in  Pig.  8  (chain-dot  curve). 

The  data  and  OUIA  calculations  at  T„  •  120  and  180  MeV  for  the  spin- 
flip  transition  to  the  3. 563-MeV  state  are  shown  in  Fig.  9.  At  120  MeV  the 
three  spin  transition  densities  give  equivalent  shapes  for  the  angular 
distributions  In  the  range  of  our  data,  15.7°  <  8c-m<  <  48.9°.  All 
calculations  predict  the  correct  location  for  the  first  maximum  of  8„  _  •  28°. 

CsIQs 

However,  the  Cohen-Kurath  Intermediate  coupling  (solid  curve),  Donnelly  and 
Ualecka  (dashed  curve),  and  phenomenological  (chain-dot  curve)  calculations 
overestimate  the  magnitude  of  the  first  maximum  by  182,  292,  and  382, 
respectively.  The  three  spin  transition  densities  yield  similar  results  for 
the  shape  of  the  130-MeV  angular  distribution  through  the  first  mlmlmum  but 
differ  through  the  second  maximum  and  minimum.  Each  calculation  underestimates 
the  magnitude  of  the  first  maximum  by  at  least  262. 

The  »+  3.563-MeV  excitation  function  data  and  DWTA  analyses  are  presented 
In  Fig.  10.  The  momentum  transfer  for  the  excitation  function  was  determined 
to  be  q  i  109  MeV/c  from  Che  first  maximum  of  the  120-MeV  angular  distribution. 
All  other  data  points  for  the  various  energies  are  within  +42  of  this  q.  None 
of  che  spin  transition  densities  predict  the  measured  shape  and  magnitude  of 
Che  3.563-MeV  excitation  function.  The  Cohen-Kurath  Intermediate  coupling, 
Donnelly  and  Ualecka,  and  phenomenological  calculations  disagree  with  che  data 
by  at  least  *202  at  both  the  low  and  high  Incident  pton  energies. 


V.  DISCUSSION  AND  RESULTS 


The  dl<«|n«Mnc  between  theory  end  dete  for  the  120-  end  180^1eV  angular 
distribution*  end  exeltetlon  function  for  the  unneturel-perlty  treneltton  to 
the  3.563-HeV  state  la  difficult  to  understand,  as  Is  the  disagreement  observed 
for  the  unnatural-parity  transition  to  the  15. 11-KeV  state  of  12C.  For 
comparison,  the  12.71-  and  15. 11-MeV  excitation  functions  along  with  DUIA 
calculation*  as  described  in  Ref.  10  are  shown  in  Fig.  11.  As  noted  in  the 
Introduction,  a  simple  OWIA  description  of  plon^iucleus  Inelastic  scattering, 
which  uses  a  single  component  of  the  plon-siucleon  Interaction  and  a  single 
transition  density,  has  been  successful  in  describing  unnatural-parity 
transitions  In  many  nuclei.  Using  the  zero-range  spln-orblt  operator  and  a 
spin  transition  density  derived  from  Cohen-Kurath  p-shell  wave  functions, 
Morris,  et  al.10  adequately  reproduced  the  experimental  12.71-MeV  excitation 
function.  Furthermore,  Cottlngame,  et  al.12  adequately  describe  the  12.71-MeV 
angular  distributions  for  «  100  to  260  MeV.  Even  though  there  Is  an  energy- 
dependent  enhancement  near  *  180  MeV  In  the  15.11-MeV  excitation  function,10 
the  low-energy  (T  »  100  and  116  MeV)  angular  distributions  for  this  stace  are 
reproduced  by  simple  DWIA  calculations. 12  Using  the  same  zero-range  spln-orblt 
operator  and  three  different  spin  transition  densities,  we  cannot  adequately 
describe  either  the  120-MeV  and  180-MeV  angular  distributions  or  the  excitation 
function  of  the  transition  to  the  3.  563-MeV  state.  However,  we  are  able  to 
reproduce  the  angular  distributions  and  constant-q  differential  cross  sections 
at  energies  from  •  100  to  260  MeV  for  the  natural-parity  transitions  to  the 
2.185-  and  &.25-aieV  levels  of  6L1  using  only  the  central  component  of  the  plon- 
nucleon  Interaction  and  transition  densities  derived  from  Cohen-Kurath  p-shell 
wave  functions.  The  configuration-space,  zero-range  form  of  the  spin-dependent 
component  of  the  plon-nucleon  Interaction  of  Carr,  et  al.°  Is  most  likely  not 
In  error.  Furthermore,  the  transition  densities  we  used  for  the  transition  to 
the  3.563-MdeV  level  yield  Inelastic  Ml  form  factors  which  agree  with  the  (e,e') 
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data  Cor  q  <  1.4  fa-1  (aaa  Fig.  8),  which  la  a  raoga  of  q  chat  adequately 
covera  the  q  ■  0.55  fa”1  at  which  the  excitation  function  wae  aeaaured. 
However,  thla  doea  not  allainata  uncertainties  In  the  spin  transition 
denaltlaa  alnce  the  electron-scattering  Inelastic  Ml  fora  factor  depends  upon 
both  the  orbital  and  spin  transition  densities.  Petrovich,  et  al.32  noted  that 
the  aonopole  spin  transition  densities  derived  froa  both  the  pure  LS-coupllng 
and  the  Cohen-Kurath  Intermediate  coupling  p-shell  wave  functlona  produced 
Caaov-Teller  aacrlx  eleaents  17. 9Z  and  8.4Z  larger  than  the  experlaental 
values.  Also,  Petrovich,  et  al.32  did  not  obtain  agreeaent  between  theory  and 
their  6U(p,p')*Ll*0.563  MeV)  data  at  Ep  »  25  MeV,  with  the  theory  falling  to 
reproduce  either  the  shape  or  aagnltude  of  the  angular  distribution.  This 
disagreement,  however,  was  not  attributed  to  uncertainties  In  the  spin 
transition  density,  but  Petrovich,  et  al.  suggested  that  ocher  reaction 
processes  In  addition  to  the  direct,  one-step  reaction  process  were 
contributing  to  the  transition.  Furthermore,  Caaaarata  and  Donnelly,38  In 
their  study  of  the  reaction  6L1( Y,'i+)8He  near  threshold,  conclude  that  the 
ratio  of  the  orbital  and  spin  transition  densities  derived  froa  the  Donnelly 
and  Unlocks  wave  functions13  Is  probably  correct.  Thus,  as  Is  the  case  for  the 
15.11-MeV  scare  of  12C,10  the  dlsagreeaent  between  our  staple  DWIA  analyses  and 
the  3.563-MeV  angular  distributions  and  exclcacion  function  is  difficult  to 
explain  as  due  to  uncertainties  In  the  spin  transition  density  or  the  spin- 
dependent  piece  of  the  plon-nucleus  Inelastic  Interaction. 

Since  staple  DWIA  calculations  using  a  well-tested  spin-dependent 
component  of  the  plon-nucleus  Inelastic  Interaction  and  spin  transition 
densities  obealned  froa  fits  to  (e.e*)  data  fall  to  reproduce  both  the  shapes 
and  magnitudes  of  the  3.563-  and  15.11-MeV  excitation  functions,  we  suggest 
that  the  same  aechanlsa  aay  be  responsible  for  the  observed  energy-dependent 
anoaalous  shapes  in  both  sLt  and  12C.  A  possible  aechanlsa  is  che  direct 
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excitadon  of  A-h  components  of  che  excited  state  wave  functlone,  as  proposed 
In  Kef.  10.  Within  this  A-h  nodel  interpretation,  we  can  estimate  the  amount 
of  A-h  admixture  in  che  3.563-HaV  wave  function.  Using  the  peak  of  the 
3.563-HaV  excitation  function  (T,  •  190  HeV),  the  DWIA  calculation  which  used 
the  phenomenologlcsl  spin  transition  density,  and  Che  procedure  of  Refs.  10  and 
39,  we  estimate  a  mixing  matrix  element  of  3  <  <  <  37  HeV.  This  range  of 
values  of  <  Is  similar  to  that  estimated  for  the  lS.U-MeV  state  of  lzC 
(7.4  <  <  <  27.3  HeV).39  However,  if  the  direct  excltadon  of  A-h  components  is 
applicable  to  both  transitions,  comparison  of  the  two  excitation  functions  and 
DWIA  calculations  (see  Figs,  10  and  11)  Indicates  that  the  resonant  A-h 
scattering  amplitude  may  Interfere  differently  with  the  p-h  scattering 
amplitude  for  the  two  transitions. 


VI.  CONCLUSIONS 

Describing  plon-nucleus  Inelastic  scaccerlng  with  che  Born  approximation 
and  DWIA,  we  obtain  predictions  which  agree  wlch  che  n+-lnelasclc-scacterlng 
measurements  for  che  3+,  2.185-HeV  and  2*,  4.25-HeV  states  of  6L1,  Our 

predictions  fall  to  reproduce  both  the  120-  and  ISO-HeV  angular  distributions 
and  che  excitation  function  for  the  unnacuraL-parlty  transition  to  che 
3.563-HeV  state.  The  measured  excltadon  function  exhibits  an  anomalous 

energy-dependent  shape  similar  to  that  observed  for  the  1*,  15.11-HeV  state  of 
1ZC.  Similar  mechanisms  may  be  responsible  for  the  anomalous  excitation 
functions  for  these  Isovector  transitions.  If  this  mechanism  is  the  direct 
excitation  of  A-h  components  in  the  final  state  wave  functions,  the  peak  of  che 
3.363-HeV  excltadon  function  is  reproduced  with  an  esdmated  mixing  matrix 

element  of  3  <  <  «  37  HeV,  a  range  of  values  consistent  wlch  che  esdmated 
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aixlttt  eacrlx  element  of  7.4  <  <  <  27.3  MeV39  for  the  13.1l-MoV  excitation 
function. 
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TABLE  I:  Differencial  croee  sections  for  w  elastic  scattering  and  Inelastic 
scattering  to  the  2.185-MeV,  0+,  3.563-MeV,  and  2*,  4.25-*eV  levels  of  *Li. 
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The  q  values  were  calculated  for  an  excitation  energy  of  2.1*5  MeV 
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FIGURE  CAPTIONS 

Fig.  1.  A  typical  energy-loss  spectrum  taken  at  T,  -  140  MeV  and  -  26° 

Cor  which  tha  hardware  veto  of  elaatic  events  was  used.  The  Insert  is 
the  energy-loss  spectrua  without  the  uaa  of  the  hardware  veto. 

Fig.  2.  Differential  croas  sections  for  n*  and  »“  elaatic  scattering  for  SL1 
for  Tf  •  120,  164,  and  180  MeV.  Tha  calculation  include  a  -30  MeV 
shift  in  tha  energy  at  which  tha  optical  aodel  paraaaters  are 
calculated.  The  164-MeV  data  are  froa  Ref.  27. 

Fig.  3.  F^(q)  for  the  3*,  T-0,  2.185-MeV  scate  of  6L1.  The  calculations  used 

microscopic  transition  densities  derived  from  pure  LS-coupllng  (solid 
curve)  and  Cohen-Kurath  intermediate  coupling  (dashed  curve)  p-shell 
wave  functions  with  a  *  0.534  fm~l  and  a  renormalization  constant  of 
2.03  and  a  «  0.556  fm~l  and  a  renormalization  constant  of  1.93, 
respectively.  The  data  are  from  Refs.  14  (solid  circles)  and  34  (open 
circles) . 

Fig.  4.  Angular  distributions  for  inelastic  scattering  to  the  3*,  T»0, 

2.185-HeV  state  of  5Li  for  T„  •  120  and  180  MeV.  The  calculations 
used  microscopic  transition  densities  derived  from  pure  LS-coupllng 
(solid  curve)  and  Cohen-Kurath  intermediate  coupling  (dashed  curve)  p- 
shell  wave  functions  with  a  •  0.534  fm_1  and  a  renormalization 
constant  of  2.03  and  a  ■  0.558  fm~^  and  a  renormalization  constant  of 
1.93,  respectively. 

Fig.  5.  Differential  cross  sections  for  **  inelastic  scattering  to  the  l*. 

T»0,  2.185-MeV  state  of  6Ul  at  a  constant  q  *  109  MeV/c.  The 
calculations  used  microscopic  transition  densities  derived  froa  pure 
LS-coupllng  (solid  curve)  and  Cohen-Kurath  intermediate  coupling 

(dashed  curve)  p-shell  wave  functions  with  a  •  0.534  fm“l  and  a 
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Fig.  6 
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rtaoraallsadoa  constant  of  2.03  and  a  «  0.538  fa-1  and  a 
renormalization  constant  of  1.93,  respectively. 

.  Angular  distributions  for  v+  lnalaatic  scattaring  to  tha  2+,  T»0, 
4.25-MaV  atata  of  ®L1  for  t(  *  120  and  180  MeV.  Tha  calculation  usad 
a  alcroscoplc  transition  density  derlvad  from  Cohen-Kurath 
lntaroadlata  coupling  p-shell  wavs  functions  with  a  -  0.52  fa-1  and  a 
ranoraallsatlon  constant  of  0.82. 

.  Differential  cross  sactlons  for  w*  inelastic  scattering  to  tha  2*. 
T-0,  4.25-MeV  state  of  ®L1  at  q  *  109  MeV/c.  Tha  calculation  usad  a 
alcroscoplc  transition  density  derived  froa  Cohen-Kurath  lnteraadlate 
coupling  p-shell  wave  functions  with  u  -  0.52  fa"1  and  a 
renormalization  constant  of  0.82. 

.  F^(q)  for  the  0*.  T*l,  3.563-MeV  state  of  SL1.  The  calculations  usad 

a  microscopic  transition  density  derived  from  Cohen-Kurath 
Intermediate  coupling  p-shell  wave  functions  (solid  curve)  with 
o  ■  0.518  fm"1  and  a  renormalization  constant  of  0.97,  a  microscopic 
transition  density  derived  from  Che  empirical  shell-model  wave 
functions  of  Donnelly  and  Walecka  (dashed  curve)  with  a  »  0.493  fm”1 
and  a  renormalization  constant  of  0.96,  and  a  phenomenological 
transition  density  based  on  the  work  of  Bergstrom,  et  al .  (chaln-doc 
curve).  Tha  data  are  from  Refs.  14  (solid  circles)  and  37  (open 
circles). 

.  Angular  distributions  for  «*  inelastic  scattering  to  the  04’,  T»l, 
3.563-MeV  state  of  6L1  for  T„  »  120  and  180  MeV.  The  calculations 
used  a  microscopic  transition  density  derived  from  Cohen-Kurath 
Intermediate  coupling  p-shsll  wave  functions  (solid  curve)  with 
t>  ■  0.518  fm"1  and  a  renormalization  conscanc  of  0.97,  a  alcroscoplc 


transition  density  derived  from  the  empirical  shell-model  wave 


function*  of  Donnelly  and  Walecka  (daahad  curva)  with  a  -  0.493  fa-1 
and  a  renormalization  con* tan t  of  0.96,  and  a  phanooanologlcal 
tranaltlon  danalty  baaed  on  the  work  of  Bergstrom,  at  al.  (chain-dot 
curve). 

Fig.  10.  Excitation  function  at  a  conatant  q  »  109  HaV/c  for  **  lnelaatlc 
acaccertng  to  the  0* ,  T-l,  3.563-MaV  atac*  of  6Li.  The  calculation* 
uaed  a  microscopic  tranaltlon  denalty  derived  froa  Cohen-Kurath 
Intermediate  coupling  p-ahell  wave  functlona  (aolld  curva)  with 
a  -  0.518  fa~L  and  a  renormalization  conacant  of  0.97,  a  microscopic 
transition  density  derived  from  tha  empirical  shell-model  wave 
functions  of  Donnelly  and  Walecka  (dashed  curve)  with  a  -  0.493  frn-1 
and  a  renormalization  constant  of  0.96,  and  a  phenomenological 
transition  density  based  on  the  work  of  Bergstrom,  et  al.  (chain-dot 
curve). 

Fig.  11.  Exclcatlon  functlona  (averaged  »+  and  x~  cross  sections  to  remove 
isospin  mixing)  at  a  constant  q  •  124  HeV/c  for  pion  Inelastic 
scattering  to  the  l+,  T-0,  12.71-MeV  (circles)  and  1+,  T»l,  15.11-MeV 
(squares)  states  of  ^C.  The  calculations  used  microscopic  transition 
densities  derived  from  Cohen-Kurath  p-sheil  wave  functions.  The  solid 
curve  la  for  the  12.71-MeV  state.  The  dashed  curve  Is  for  the  15.11— 
MeV  state  and  has  been  multiplied  by  four  (see  Ref.  L0). 
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